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SECTION  1.0 

INTRODUCTION  AND  SUMMARY 

This  report  describes  results  of  an  Investigation  of  radiation  effects 
on  semiconductor  devices  used  In  electro-optical  sensor  applications.  Emphasis 
Is  placed  on  determining  the  basic  mechanisms  of  the  Interaction  of  radiation 
with  such  devices  with  a  view  toward  gair.  .g  understanding  of  benefit  to  devel¬ 
opers  of  radiation-tolerant  devices.  Sectl*'  2.0  describes  results  of  a  study 
of  ionizing  radiation  effects  on  silicon  MOS  devices  at  cryogenic  temperatures. 
In  Section  3.0,  findings  of  an  investigation  of  the  mechanisms  of  radiation 
effects  on  extrinsic  silicon  detector  material  are  presented.  Section  4.0 
describes  results  of  a  study  of  neutron-  and  ioni zing-radiation  cf  facts  on 
charge -coupled  devices.  In  Section  5.0,  the  topic  of  a  neutron  damage  coeffi¬ 
cient  appropriate  for  silicon  depletion  regions  is  addressed.  The  remainder 
of  the  present  section  sunmarizes  major  findings  reported  In  this  document. 

A  study  of  ionizing  radiation  effects  on  silicon  MOS  devices  was  per¬ 
formed  at  temperatures  ranging  from  4  to  77°K.  Emphasis  was  placed  on  com¬ 
paring  behavior  at  4°K  with  that  at  77°K  for  the  purpose  of  determining  whether 
electron  trapping  is  important  at.  the  lower  temperature.  Various  experimental 
techniques  were  employed  to  obtain  appropriate  pre-  and  post-irradiation  elec¬ 
trical  data  on  the  MOS  capacitors  employed  as  test  vehicles.  Conventional  C-V 
measurement  techniques  are  inapplicable  at  4°K  because  of  carrier  freeze-out. 

It  is  therefore  necessary  to  generate  sufficient  majority  carriers  to  maintain 
charge  neutrality  as  the  silicon  substrate  is  swept  from  accumulation  to  deple¬ 
tion.  This  was  accomplished  here  by  two  approaches:  1)  use  of  optical  excita¬ 
tion;  2)  use  of  the  Co*^  gamma  rays  employed  for  irradiation  as  a  source  of 
carriers,  thereby  permitting  C-V  measurements  during  bombardment.  Comparison 
C-V  measurements  were  made  at  18-20°K  using  solely  thermal  excitation  of  car¬ 
riers.  It  is  the  compelling  conclusion  of  this  work  that  charge  buildup  at 
4°K  in  irradiated  SI 0^  films  can  be  described  exclusively  in  terms  of  hole 
trapping  since  results  observed  at  77°K,  which  previously  have  been  success¬ 
fully  modeled  on  this  basis,  are  identical  to  the  lower^ temperature  findings. 
Strictly  speaking,  this  conclusion  should  be  modified  to  apply  specifically  to 
the  case  of  the  wet  oxides  studied  and  to  the  particular  irradiation  conditions 
employed  (steady-state  ionizing  radiation--20  minute  time  scale--10^  rads(Si)). 
It  remains  a  possibility  that  significant  electron  trapping  may  be  observed  on 


shorter  time  scales,  and  this  Issue  should  be  addressed.  We  also  raise  the 
possibility  that  electron  trapping  and  transport  In  Si 02  may  be  dependent  on 
oxide  processing  conditions. 

An  investigation  of  the  basic  mechanisms  of  radiation  effects  on  extrinsic 
silicon  detectors  was  initiated.  Such  devices  are  becoming  technologically  im¬ 
portant  because  of  their  suitability  for  the  fabrication  of  monolithic  focal - 
plane  arrays  to  be  utilized  in  infrared  imaging  applications.  The  work  reported 
here  constitutes  the  beginning  of  an  effort  that  will  continue  in  a  subsequent 
program  during  1979.  We  first  provide  a  detailed  background  discussion  of  radi¬ 
ation  effects  on  extrinsic  silicon.  Experimental  results  obtained  on  neutron- 
irradiated  S':S,  Si:Zn,  and  Si:Au  are  then  presented.  Hall  measurements  were 
performed  as  a  function  of  temperature  before  and  after  bombardment  in  order  to 
determine  carrier  removal  rates.  It  was  found  that  the  introduction  rate  of 
compensating  defeats  depends  strongly  on  the  depth  within  the  forbidden  gap  of 
the  intentional  dopant  level.  Specifically,  a  donor  level  0.19  eV  from  the 
conduction  band  (sulfur)  appears  to  be  compensated  nearly  as  effectively  as  a 
shallow  donor  (arsenic).  On  the  other  hand,  the  carrier  removal  rate  was  found 
to  be  considerably  smaller  for  an  acceptor  level  0.31  eV  from  the  valence  band 
(zinc).  Interpretations  of  the  data  obtained  to  date  are  given,  as  are  sug¬ 
gestions  for  relevant  future  investigations  in  this  area. 

A  detailed  study  of  radiation  effects  on  charge-couplcJ  devices  has  oeen 
performed,  with  emphasis  placed  on  examining  changes  in  dark  current  density 
produced  by  bombardment  with  neutrons  and  with  ionizing  radiation.  Devices 
from  two  manufacturers  were  examined:  CCDs  and  MOS  capacitors  from  Hughes  Air¬ 
craft  (HAC)  and  CCDs  from  Fairchild  Semiconductors  (CCD  311).  The  MOS  capaci¬ 
tors  were  used  to  obtain  a  neutron  damage  coefficient  appropriate  for  describing 
increases  in  the  generation-center  density  in  silicon  depletion  regions.  Capac¬ 
itors  were  also  employed  to  obtain  gamma-ray  damage  coefficients  associated  with 
increases  in  surface  generation  velocity  at  the  SiO^-Si  interface.  The  rate  of 
increase  in  dark  current  density  for  neutron- irradiated  CCDs  was  measured  and 
found  to  be  very  similar  for  HAC  and  Fairchild  CCDs  (3.5  x  10-*1  nanoamps  per 
fission  neutron  for  the  Fairchild  device).  A  discussion  is  presented  which 
indicates  that  a  nx  derate  reduction  in  the  operating  temperature  of  a  CCD  should 
substantially  reduce  its  sensitivity  to  the  production  of  dark  current  by  neutron 
bombardment.  The  HAC  CCDs  studied  were  found  to  be  quite  sensitive  to  the  effects 


10 


TfnTr 


of  ionizing  radiation  when  irradiated  under  power.  Short-term  annealing 
measurements  were  performed  in  which  dark  current  density  in  the  CCD  311  was 
monitored  as  a  function  of  time  following  pulsed  neutron  bombardment.  As  an 
example,  at  20  sec  after  irradiation  the  dark  current  annealing  factor  was 
~1.7.  This  means  that,  at  a  given  fluence,  the  neutron- induced  increase  in 
dark  current  density  is  a  factor  of  1.7  larger  at  20  sec  after  bombardment 
than  it  is  at  some  much  longer  time  (such  as  30  minutes). 

In  the  MOS  capacitor  studies  mentioned  above,  a  generation-lifetime  dam¬ 
age  coefficient  (Kg)  was  obtained  that  is  appropriate  for  introduction  of  damage 

in  silicon  depletion  regions  by  neutron  bombardment.  The  value  obtained,  7.0 
6  2 

x  10  n-sec/cm,  appears  to  apply  for  depletion  regions  in  both  n-  and  p-type 
silicon  at  room  temperature  and  also  is  a  good  value  for  use  over  a  moderate 
temperature  range  (~240  to  ~370°K,  the  range  studied).  We  explored  the  possi¬ 
bility  that  this  damage  coefficient  might  be  useful  for  predicting  the  effects 
of  neutron  bombardment  on  silicon  devices.  We  have  found  that  leakage  current 
data  obtained  for  various  neutron-irradiated  silicon  devices,  including  CCDs, 
PIN  diodes,  and  JFETs,  can  be  described  reasonably  well  using  a  model  which  in¬ 
corporates  the  unique  value  of  Kg  obtained  in  this  investigation.  Agreement 
between  calculations  and  experimental  results  is  within  a  factor  of  5  2  in 
most  cases,  which  suggests  that  the  calculational  approach  employed  is  suit¬ 
able  for  making  reasonably  accurate  estimates  of  neutron- induced  increases  in 
dark  current  for  silicon  devices. 
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IONIZING  RADIATION  EFFECTS  ON  MOS  DEVICES  AT  CRYOGENIC  TEMPERATURES 
2.1  INTRODUCTION  AND  BACKGROUND 

The  transport  properties  of  holes  in  silicon  dioxide  films  grown  on  sili- 

J  .v? 

con  were  the  subject  of  study  in  our  laboratory  several  year*  aqb,  and  V 

mechanisms  of  transport  are  now  reasonably  well  understood,  at  least  to  the 
point  where  predictions  can  be  made  for  the  performance  of  MOS  devices  operated 
at  cryogenic  temperatures.  *  The  hole  conduction  process  was  found  to  be  gov¬ 
erned  by  traps  distributed  in  energy  and  space.  This  process  was  therefore 
strongly  thermally  activated,  so  that  near-permanent  trapping  of  holes  occurs 
at  low-to-moderate  applied  fields  at  cryogenic  temperatures  (77°K). 

In  the  above  studies,  which  were  performed  at  77°K  or  above,  electrons 
were  found  to  be  transported  through  the  Si 0^  films  on  a  time  scale  short  com¬ 
pared  to  the  experimental  resolution  (10  p,sec).  Consequently,  there  was  no 
evidence  of  trap-limited  conduction  for  electrons.  Moreover,  previous  studies 

of  electron  transport  through  Si0o  films  revealed  that  the  electron  mobility 

^  2 

was  quite  high  for  amorphous  material  (on  the  order  of  20  cm  /V-sec),  and  ex¬ 
hibited  a  decrease  with  increasing  temperature,  which  is  consistent  with  domi- 
nance  by  scattering  processes  rather  than  activation  processes.  It  thus  ap¬ 
peared  that  electron  transport  was  of  a  fundamentally  different  character  than 
nole  transport.  Electrons  are  apparently  governed  by  their  intrinsic,  or 
scattering-limited,  mobility  whereas  holes  drift  with  their  characteristic 
scattering-limited  mobility  only  on  a  time  scale  of  nanoseconds  before  being 
trapped.  It  is  also  possible,  however,  that  trap-limited  conduction  of  elec¬ 
trons  might  manifest  itself  at  temperatures  low  compared  to  the  previous  region 
of  measurement.  If  that  is  in  fact  the  case,  then  operation  of  devices  at  such 
temperatures  at  which  electrons  are  semipermanently  trapped  could  serve  as  a 
partial  remedy  for  the  problem  of  positive  charge  buildup  at  cryogenic  tempera- 
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tures.  The  compensation  of  trapped  holes  by  trapped  electrons  could  have 
practical  merit  even  if  it  were  transient.  In  any  case,  the  possibility  of 
such  trapping  in  SiC^  films  deserved  examination. 

To  make  it  plausible  that  electron  trapping  at  low  temperatures  might  be 
observable,  it  is  worthwhile  to  discuss  the  hole  conduction  process  qualitatively 
in  more  detail.  The  trap-limited  conduction  process  has  been  successfully 
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modeled  by  two  very  different  approaches.  One  is  based  on  a  multiple- trapping 
picture,  in  which  holes  transported  thit>ugh  an  oxide  film  under  the  Influence 
of  an  applied  field  incur  a  probability  of  being  trapped  instates  which  are 
distributed  in  energy  within  the  band  gap.7  The  deeper  states  have  longer  de- 
trapping  times,  so  that  a  given  hole  population  which  Is  suddenly  Introduced 
progressively  distributes  Itself  in  the  deeper  states,  thus  reducing  the  ap¬ 
parent  mobility  with  time,  and  giving  rise  to  exceedingly  "dispersive"  trans¬ 
port.  The  trap  density  of  states  which  fits  the  data  for  the  oxide  studied  in 
detail  in  our  laboratory  is  a  simple  one  in  which  the  density  is  an  exponenti¬ 
ally  declining  function  of  depth  In  energy  from  the  band  edge.7  Such  a  density 
of  states  function  has  the  virtue  of  being  physically  reasonable  tor  a  disord¬ 
ered  solid. 

The  alternative  model  for  the  description  of  hole  transport  is  the  purely 

statistical  stochastic  transport  model  of  Scher  and  Montroll,  as  applied  to  the 

8- 10 

case  of  Sit^  films  by  McLean  et  al .  That  model,  referred  to  as  Continuous 
Time  Random  Walk  (CTRW),  does  not  presuppose  a  physical  mechanism  for  the  ori¬ 
gin  of  the  dispersion  in  carrier  transport.  In  particular,  it  has  been  shown 
to  encompass  the  multiple-trapping  model  as  a  particular  case  under  certain 
conditions.11  In  general,  however,  the  dispersion  in  transport  can  arise  from 
a  number  of  possible  mechanisms.  A  generalized  description  has  the  advantage 
that  all  manner  of  observed  dispersive  transport  processes  can  be  treated 
within  the  same  framework,  but.  one  difficulty  is  that  the  parameters  of  such  a 
model  are  not  easily  associated  with  unique  physical  interpretations. 

R.  C.  Hughes  has  l  ikewise  treated  the  hole  transport,  problem  in  S3 
films  in  terms  of  the  CTRW  model,  and  ascribed  deviations  f^om  that  model  at 
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early  tunes  to  the  manifestation  of  intrinsic  hole  mobility  prior  to  trapping. 
This  intrinsic  mobility  is  also  found  to  be  thermally  activated,  and  to  be  rel¬ 
atively  small  compared  to  electron  mobility.  The  activation  energy  of  about 

q 

0.16  eV  is  taken  to  imply  polaronic  hopping  motion.  McLean  et  al.  observed 
charge  transport  governed  by  activation  energies  of  less  than  0.1  eV  at  low 
temperatures  which  were  likewise  interpreted  in  terms  of  polaron  motion.  The 
higher  temperature  behavior  is  governed  by  activation  energies  similar  to  those 
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observed  in  an  oxide  studied  at  Northrop  (~  0.4  eV). 

If  it  car.  be  established  that  the  hole  transport  data  necessarily  imply  a 
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spread  in  the  density  of  states  at  the  valence  band  edge  of  a  certain  magnitude, 
then  It  is  possible  to  argue  that  the  same  mechanisms  give  rise  to  a  distribution 
of  states  at  the  conduction  band  edge,  although  the  width  in  energy  may  be  differ- 
|  ent.  In  the  multi  pie- trapping  model  of  hole  transport,  the  trap  energy  levels  which 

•  are  found  are  measured  with  respect  to  the  valence  band  edge,  and  the  Implication 

I  Is  that  the  density  of  states  falls  off  exponentially  from  this  edge,  with  a  spread 

I  on  the  order  of  0.4  eV.  On  the  other  hand,  a  number  of  different  mechanisms  of  trap- 

r  limited  transport  could  be  subsumed  under  the  CTRW  model.  It  the  polaron  model  is 

val id--and  the  weight  of  evidence  in  papers  referenced  above  appears  to  be  strongly 
In  favor  of  small -polaron  formation  for  holes--then  the  dominant  activation  ener¬ 
gies  characterizing  transport  (0.4  eV)  would  be  interpreted  in  terms  of  the  polaron 
binding  energy.  A  spread  in  the  observed  activation  energies  could  be  interpreted 
in  terms  of  a  distribution  of  band-tail  states  ascribable  to  disorder,  but  an  in¬ 
terpretation  in  terms  of  a  site-dependence  of  the  polaron  binding  energy  may  be 
more  appropriate,  in  which  case  no  inference  could  be  readily  drawn  for  the  den¬ 
sity  of  states  at  the  valence  band  edge.  However,  even  in  the  event  that  no  ap¬ 
peal  can  be  made  to  our  knowledge  about  the  valence  band  edge,  general  arguments 
can  be  given  for  expecting  at  least  some  tail  in  the  density  of  states  at  the  con¬ 
duction  band,  and  these  are  given  below. 

It  is  helpful  to  recall  at  this  point  that  the  oxides  which  were  originally 
fitted  in  terms  of  the  multi  pie- trapping  model  arid  the  CTRW  model,  respectively, 
were  in  fact  quite  different  in  their  dispersive  transport  behavior.  The  oxide 
studied  at  Northrop  indicated  a  wide  distribution  of  energy  levels,  whereas  the 
oxides  which  yielded  the  best  fit  to  the  cne-parameter  CTRW  model  showed  very 
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little  spread  in  activation  energies.  The  implication  is  that  processing  differ¬ 
ences  between  different  oxides  could  result  in  a  considerably  different  spread  of 
the  density  of  states  at  the  valence  band  edge,  and,  by  implication,  at  the  con- 
auction  band  edge. 

It  has  been  argued  above  that  a  spread  of  energies  at  ♦he  valence  band  edge 
has  similar  implications  for  the  conduction  band  edge,  and  yet  it  is  clear  from 
the  electron  transport  behavior  already  observed  that  the  "scale"  of  this  spread 
Js  very  different  in  the  two  cases.  Molecular  orbital  theory  can  be  used  to  ex¬ 
plain  this  difference.  Electronic  structure  calculations  have  been  performed  on 
quartz  which  showed  that  changes  in  the  Si-O-Si  bond  angle  introduce  significant 
changes  in  the  valence  band  edge.  A  distribution  cf  such  bond  angles,  as  is 
likely  to  be  present  in  amorphous  SiOgj  is  therefore  sufficient  to  explain  the 
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spread  of  energies  Implied  by  the  multiple  trapping  picture  (0.4  etf).  The  effect 
of  this  variation  In  bond  angle  on  the  conduction  band  edge  is  considerably  re¬ 
duced,  according  to  the  molecular  orbital  picture  Similar  conclusions  may  be 

drawn  from  another  study  of  the  electronic  structure  of  amorphous  SiO«j,  quartz 

14  *• 

and  crlstobal ite.  The  effect  of  bond  angle  variations  on  the  valence  band  den¬ 
sity  of  states  is  greater  than  that  of  topological  disorder.  By  contrast,  the 
conduction  band  density  of  states  reveals  a  relatively  larger  influence  of  dis¬ 
order. 

In  addition  to  the  above  very  general  arguments,  we  may  cite  the  observation 

of  a  tail  in  the  conduct ion- band  density  of  states  by  means  of  photon-assisted  tun- 
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neling  measurements  of  electrons  from  aluminum  into  the  Sif^  conduction  band. 

A  Gaussian  tail  of  width  O.leVgave  a  fit  to  the  data  in  most  cases,  except  for  an 
rnarmealed  oxide,  where  a  width  of  0.25  eV  was  indicated.  A  profound  influence  of 
processing  on  the  apparent  band  tail  is  suggested  by  these  findings.  The  magnitude 
of  the  tail,  however,  is  ruch  that  one  would  expect  to  see  trapping-dominated  mo¬ 
bility  even  at  room  temperature,  unless  the  tail  spates  are  likewise  "extended" 
states,  or  characteristic  only  of  the  interface.  It  is  possible  in  the  latter  case 
that  the  applied  field  is  sufficient  to  promote  electrons  (which  have  tunneled  into 
such  tail  states)  into  extended  states  of  the  buli 

The  above  considerations  justify  the  attempt  t(  ascertain  the  existence  of 
shallow  electron  trapping.  That  the  trapping  must  be  shallow  is  clear  from  the  ob¬ 
servation  of  scattering-limited  conduction  at  room  temperature  on  the  one  hand, 
and  the  fact  that  trapping  does  not  manifest  itself  at  times  <10  nsec  near  77°K 
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on  the  other. 

In  the  present  study,  charge  buildup  measurements  nder  gamma  irradiation 
were  employed  to  reveal  whetner  compensation  of  trapped  holes  by  electrons  occurs 
at  4°K.  A  number  of  different  experimental  methods  were  utilized  and  these  will 
be  discussed  in  turn,  along  with  the  results  obtained  with  each. 

2.2  EXPERIMENTAL  METHODS  AND  RESULTS 

2.2.1  Measurements  of  Flatband  Shifts  at  4°K  Ui  ing  Optical  Excitation 

Charge  buildup  measurements  were  performed  on  room- temperature  radiation- 
hard  oxides  grown  pyrogenically  at  925°C  to  various  thicknesses  by  Hughes  Air¬ 
craft.  Company  (Lot  KA-09).  Dry  oxide  films  (1000°C  growth)  had  previously  been 
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characterized  in  our  laboratory  at  77  K.  ’  Standard  capacitance-versus-voltage 
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data  were  obtained  at  1  MHz  In  that  case  to  determine  the  charge  buildup  via 
shifts  In  the  flatband  voltage.  Corresponding  measurements  were  now  required 
at  4°K.  In  the  latter  case,  however,  carrier  freeze-out  renders  conventional 
techniques  Inapplicable.  It  Is  therefore  necessary  to  generate  in  some  manner 
sufficient  majority  carriers  to  maintain  charge  neutrality  conditions  in  the 
bulk  as  the  gate  voltage  is  swept  from  accumulation  to  depletion.  That  Is, 
the  dielectric  relaxation  frequency  (=  pc,  where  p  Is  the  resistivity  and  c 
the  dielectric  constant)  must  be  sufficiently  high  for  the  majority  carriers 

to  follow  the  applied  ramp  waveform.  This  requires  a  silicon  resistivity  of 
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no  more  than  about  10  ohm-cm,  or  a  carrier  concentration  on  the  order  of  10 

cm“^.  A  second  criterion  is  imposed  by  the  requirement  that  the  majority  car¬ 
riers  must  be  able  to  maintain  equilibrium  conditions  of  depletion  width.  That 
is,  the  rate  of  capture  and  emission  of  majority  carriers  for  the  shallow  donor 
sites  must  be  sufficiently  rapid. 

In  the  initial  experiment  of  this  study,  optical  excitation  was  used  to 
furnish  majority  carriers,  and,  to  minimize  the  degree  of  excitation  required, 
C-V  measurements  were  made  at  3  Hz  rather  than  1  MHz.  In  practice,  the  optical 
signal  level  was  increased  until  the  3-Hz  capacitance  signal  in  accumulation 
nearly  saturated,  and  thus  the  dielectric  relaxation  frequency  increased  to 
approximately  3  Hz.  The  optical  excitation  was  furnished  by  an  incandescent 
filament  (removed  from  a  flashlight  bulb)  which  was  placed  in  proximity  to  the 
sample  (about  3  nm).  The  sample  and  filament  were  placed  in  a  heliuri  exchange 
gas  environment  to  permit  variation  of  sample  temperature,  and  this  prevented 
premature  oxidation  of  the  filament. 
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Samples  were  exposed  to  a  total  dose  of  10  rads  (±  10%)  in  the  exposure 
room  of  the  Northrop  Co-60  source.  This  required  a  period  of  20  minutes,  not 
including  a  period  of  about  5  minutes  for  raising,  and  5  minutes  for  lowering, 
the  source.  Irradiation  wac  performed  in  the  dark,  and  tha  optical  signal  was 
employed  only  during  periods  of  measurement.  Detection  of  the  3-Hz  signal  was 
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The  limiting  modulation  angular  frequency  is  given  by  w  =  cn(n  +  ni),  where 
cn  is  tne  (electron)  capture  probability,  n  is  the  equilibrium  electron  con¬ 
centration,  and  n^  is  the  electron  concentration  if  the  Fermi  level  were  at 
the  shallow  donor  energy  level  position.  No  requirements  are  placed  on  equi¬ 
librium  with  respect  to  minority  carriers  (Ref.  16). 
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by  means  of  an  Ithaco  39? A  heterodyne  loci  -In  detector.  Insufficient  phase 
stability  was  observed  for  th1*  Internal  reference  of  that  lock-in  (due  to  pre¬ 
vailing  line  voltage  transients)  and  an  HP-3310A  (or  B)  was  eiuployed  Instead, 

In  conjunction  with  a  voltage  regulator. 

The  sample  holder  provided  for  the  simultaneous  irradiation  of  four  ca¬ 
pacitors.  These  were  measured  sequential 1y  starting  four  minutes  after  the 
source  descent  was  begun,  at  which  time  the  dose  rate  was  negligible.  Due  to 
the  use  of  3-Hz  modulation,  filter  time  constants  larger  than  1  sec  had  to  be 
used,  requiring  low  ramp  rates.  About  one  minute  was  required  t**  obtain  one 
CV-curve.  In  order  to  forestall  any  field- induced  transport  during  the  meas¬ 
urement  time,  sample  biases  were  reduced  on  all  four  devices  at  the  same  4- 
minute  point.  During  the  measurement  itself,  biases  were  kept  sufficiently 
low  tha t  field-induced  hole  transport  did  not  occur.  For  the  typical  sample 
thickness  of  970  A,  the  bias  was  held  to  2 V  (2  x  10JV/cm) ,  just  sufficient  to 
display  the  flatband  shifts  which  were  typically  on  the  O  der  of  one  volt. 

The  circuit  for  the  experiment  described  is  shown  in  Figure  2-1.  Sample  tem¬ 
perature  measurements  were  made  by  means  of  a  carbon  resistor  and  a  platinum 
resistor. 

Since  the  majority  carriers  are  frozen  out  at  4°K,  the  device  equivalent 

circuit  at  this  temperature  is  that  of  a  capacitor  with  a  dual  dielectric  of 

Si02  and  silicon.  It  is  important  to  establish  what  fraction  of  the  applied 

bias  appears  across  the  oxide  in  such  a  device.  It  was  hypothesized  that  under 

irradiation  the  induced  conductivity  of  the  silicon  would  be  ar  greater  than 

that  of  the  SiC^,  and  that  therefore  the  applied  bias  would  be  dropped  almost 

entirely  across  the  oxide.  The  narrower  band  gap  of  silicon  results  in  a  pair 

generation  effectiveness  which  is  five  times  greater  than  that  of  S10?.  Both 
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carriers  are  mobile  in  silicon,  and  the  mobility  ratio  is  in  excess  of  10  with 
respect  to  estimated  electron  mobilities  in  SiOg .  The  electric  field  strength 
is  a  factor  of  five  greater  in  the  SiOr,  than  in  silicon  by  virtue  of  differ¬ 
ences  in  dielectric  constant.  On  the  other  hand,  the  silicon  thickness  is 
considerably  greater  than  that  of  the  oxide  (by  a  factor  of  10  ),  and  little 
can  be  said  about  relative  lifetimes  in  the  two  materials  a  priori.  During  the 
course  of  the  experiments  a  direct  test  of  the  relative  conductivities  under 
irradiation  confirmed  that  the  bias  is  indeed  dropped  entirely  across  the  ox¬ 
ide,  and  this  test  is  discussed  below. 
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Results  of  measurements  of  flatband  shift  after  10*  vaus  (SI)  Irradiation 
at  4°K  are  shown  as  a  function  of  gate  bias  in  Figure  2-2  for  a  sample  of  970  k 
thickness.  Shown  on  the  same  figure  are  comparable  data  taken  at  77°K  using  the 
same  experimental  method  (but  no  optical  excitation).  These  latter  data  have 
been  previously  described  successfully  in  terms  of  a  f1eld-depand*.nt  charge 
yield  at  low  field  Intensities,  rvj  hole  transport  out  of  the  oxide  at  high 
field  Intensities.*  In  that  previous  modeling  of  transport.  It  was  assumed 
that  all  electrons  /fere  swept  out  of  the  oxide  prior  to  measurement.  The  agree¬ 
ment  of  the  sets  of  data  at  the  two  temperatures  suggests  strongly  that  no  sig¬ 
nificant  electron  trapping  occurs  at  4°K  on  the  time  scale  of  the  measurement. 

Several  questions  may  be  raised  about  the  validity  and  generality  of  this 
finding.  There  is  first  of  all  the  possibility  that  the  optical  excitation 
also  resulted  In  local  heating  of  the  sample.  The  temperature-sensing  resistor 
was  located  some  5-6  mm  away  from  the  dual-in-line  package  to  which  the  device 
was  bonded.  Whereas  it  indicated  a  temperature  of  the  sample  mount  of  no  more 
than  6°K  during  measurement,  the  proximity  of  the  optical  source  to  the  device 
makes  it  difficult  to  rule  out  some  local  heating  of  the  oxide  fPm  with  re¬ 
spect  to  the  sample  holder.  An  upper  bound  of  less  than  19°K  for  the  resultant 
sample  temperature  is  implied  by  the  fact  that  carrier  excitation  was  optical 
rather  than  thermal  in  nature.  Thermal  excitation  of  sufficient  carriers  to 
permit  a  CV  measurement  would  have  required  a  silicon  temperature  of  about 
19°R.  The  conclusion  that  the  excitation  was  in  fact  optical  rather  than 
thermal  was  breed  on  the  observation  that  the  power  which  had  to  be  supplied 
to  the  filament  to  obtain  the  required  excitation  in  the  silicon  gradually 
increased  over  time,  suggestive  of  gradual  oxide  cion  of  the  filament.  Since 
the  thermal  resistance  between  the  filament  and  the  device  should  not  be  a 
'unction  of  time,  the  excitation  mechanism  had  to  be  optical  rather  than 
thermal . 

A  second  question  about  this  measurement  concerns  possible  excitation  of 
electrons  out  of  shallow  traps  by  means  of  the  optical  signal  used  to  generate 
the  ir »'ority  carriers.  The  fact  that  we  did  not  observe  any  change  In  flatband 
shift  during  the  course  of  measurement,  however,  suggests  that  excitation  and 
sweep  out  of  trapped  charge  would  l.ave  to  occur  on  a  time  scale  short  compared 
to  a  minute,  implying  inordinately  large  photoionization  cross  sections  for 
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Figure  2-2.  Platband  voltage  shift  in  MOS  capacitors  as  a 
function  of  bias  applied  during  irradiation  to  104 
rads (Si ) .  Data  obtained  at  4°K  using  optical  excita¬ 
tion  and  at  77°K  without  optical  excitation  are  shown. 
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such  traps.  A  third  possibility  Is  that  any  shallow  electron  traps  may  be 
subject  to  field-enhanced  excitation  much  more  readily  than  the  deeper  hole 
V.'  :,ps.  Whereas  this  cannot  be  ruled  out,  it  can  be  stated  simply  that  if  elec¬ 
tion  trapping  does  not  prevail  at  typical  sample  biases,  then  it  holds  no  prac¬ 
tical  interest  in  any  case. 

2.2.2  Measurements  of  Flatband  Shift  at  18-2Q°K 

The  concern  about  possible  optical  excitation  out  of  electron  traps  moti¬ 
vated  an  alternative  approach  in  which  excitation  of  the  requisite  majority 
carrier  concentration  was  achieved  thermal l>.  The  irradiation  was  still  per¬ 
formed  at  4°K,  and  the  sample  temperature  subsequently  raised  just  enough  to 
permit  CV -measurement  at  3  Hz.  Preliminary  to  this  measurement,  the  5-Hz  ca¬ 
pacitance  signal  in  accumulation  was  measured  as  a  function  of  temperature  in 
order  to  ascertain  the  temperature  dependence  of  the  majority  carrier  concen¬ 
tration,  or,  more  directly,  of  the  dielectric  relaxation  frequency.  These  data 
were  obtained  with  an  Ithaco  395  narrowband  voltmeter,  which  is  preferred  here 
by  virtue  of  the  fact  that  it  yields  the  magnitude  of  tiie  signal  independent  of 
signal  phase.  (The  magnitude  is  derived  from  the  vector  of  sum  of  the  outputs 
of  two  phase-sensitive  detectors  in  quadrature.)  5  Hz  was  the  lowest  frequency 
of  operation  at  which  phaselock  to  the  reference  could  be  maintained.  Results 
are  shown  in  Figure  2-3.  Saturation  of  the  capacitance  signal  is  seen  to  re¬ 
quire  about  19-20°K  for  this  sample.  Results  of  measurements  made  using  ther¬ 
mal  excitation  during  measurement  are  shown  in  Figure  2-4.  Comparison  data 
taker  at  77°K  are  shown  in  Figure  2-5,  and  the  composite  of  these  data  sets  is 
illustrated  in  Figure  2-5.  Again,  results  indicate  no  significant  difference 
between  the  charge  buildup  observed  at  77°K  (where  electrons  clearly  do  not 
play  a  role)  and  that  observed  after  4°K  irradiation. 

On  the  basis  of  these  results,  one  may  draw  the  conclusion  that  electron 
trapping  does  net  play  a  role  in  charge  buildup  at  19-20°K  on  the  tinie  scale  cf 
the  present  experiments.  This  temperature  is  just  sufficient  to  excite  a  small 
fraction  of  the  shallow  phosphorus  donor  levels  located  0.045  eV  from  the  con¬ 
duction  band.  The  fraction  ionized  may  be  determined  as  follows:  the  dielec¬ 
tric  relaxation  frequency  of  3  Hz  at  operating  temperature  implies  a  carrier 

3-3  4  2 

concentration  in  the  bulk  of  2  x  10  cm  '  ,  on  the  assumption  of  2  x  10  cm  /V- 

15  -3 

sec  mobility.  The  net  donor  doping  concentration  is  2  x  10  cm  ,  yielding 

- 12 

an  ionized  fraction  of  10  .  We  would  expect  similar  fractional  iomza^en 
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Figure  2-3.  5-Hz  capacitance  signal  in  accumulation  as 
a  function  of  reciprocal  temperature.  The  fall- 
off  at  low  temperature  governed  by  dielectric 
relaxation  is  illustrated. 
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figure  2-4.  Flatband  voltage  shift  in  MOS  capacitors 
as  a  function  of  bias  applied  during  irradiation 
to  l(r  rads  (Si)  at  4°K.  Measurements  were  made4 
using  thermal  excitation  (19-20°K). 
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Figure  2-5.  Hatband  voltage  snift  in  MOS  capacitors 
as  a  function  of  bias  applied  during  irradiation 
to  104  rads (Si)  at  77°K. 


of  traps  of  such  energy  level  position  with  respect  to  the  conduction  band  In 
the  S102  If  density  of  states  factors  were  the  same.  The  latter  affects  the 
argument  quantitatively  to  some  degree,  but  we  can  state  within  the  bounds  of 
this  approximation  that  electron  traps  on  the  order  of  40  meV  or  deeper  are  not 
emptied  significantly  by  raising  the  sample  temperature  to  19°K  for  measure¬ 
ment.  We  have  therefore  ruled  out  the  existence  of  a  significant  population  of 
electron  traps  (i.e.,  significant  relative  to  the  density  of  hole  traps)  deeper 
than  about  40  meV  in  this  oxide  by  the  present  test. 

We  were  required  to  operate  at  19°K  for  the  purpose  of  achieving  a  dielec¬ 
tric  relaxation  frequency  of  3  Hz  to  accommodate  CV  measurements.  The  same  in¬ 
formation  may  be  obtained  by  quasistatic  techniques,  as  first  described  by 
Kuhn.17  In  this  measurement,  the  majority  carrier  current  is  measured  as  the 
MOS  capacitor  is  swept  between  accumulation  and  depletion.  The  condition  for 
the  validity  of  such  data  is  that  the  dielectric  relaxation  frequency  be  high 
enough  for  the  majority  carriers  to  follow  the  ramp.  The  signal  bandwidth  is 
thus  reduced,  permitting  reduction  of  the  sample  temperature  cormiensurately. 

If  the  same  ramp  rate  is  used  as  in  the  above  measurement,  the  signal  band¬ 
width  is  reduced  an  order  of  magnitude  with  respect  to  the  3-Hz  CV  method,  and 
reference  to  Figure  2-3  shows  that  this  permits  operation  at  a  temperature  of 
about  17°K. 

The  experiment  was  therefore  repeated  for  a  number  of  sample  biases  using 
the  quasistatic  method.  Results  of  a  pre-irradiation  run  are  shown  in  Figure 
2-7  for  a  sample  of  970  A  oxide  thickness.  The  capacitance  voltage  curve  shows 
two  dips  between  the  accumulation  bias  condition  and  flatband.  A  dip  can  re¬ 
sult  from  the  depopulation  of  shallow  dopant  levels  as  the  Fermi  level  crosses 
the  dopant  energy  level  position.  The  appearance  of  two  dips  in  this  case  was 
not  explained.  The  dips  were  not  observed  post-irradiation.  Likewise,  they 
were  not  observed  under  3-Hz  modulation,  presumably  due  to  overmodulation  in 
that  case. 


Results  of  flatband  voltage  shifts  using  the  quasistatic  method  were  com' 
pared  with  3-Hz  modulation  data  taken  at  20°K.  The  data  are  summarized  in 
Table  2-1.  Good  agreement  between  the  two  methods  was  found;  the  differences 
that  do  exist  were  unsystematic,  and  indicative  of  the  basic  reproducibility 
of  our  methods. 
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Table  2-1.  Comparison  of  3-Hz  and  quasistatic  measurements  of  flatband  voltage 
shift  at  20°K  after  a  lO^-rad  irradiation  at  4°K. 


Device  # 

Applied  Field 

Flatband  Voltage  Shift  (Volts) 
Quasistatic  3-Hz 

1 

1.06  x  106 

1.53 

1.55 

2 

2.12  x  106 

1.82 

1.79 

3 

3.2  x  106 

1.72 

1.70 

4 

4.26  x  106 

0.74 

0.78 

J 

A  second  measurement  was  therefore  made,  in  which  the  quasi  static  data 
were  taken  at  the  lowest  convenient  temperature.  Under  conditions  of  a  sweep 
rate  no  less  than  100  mV/sec,  hysteresis  in  the  quasistatic  C V  curve  manda  -d 
an  operating  temperature  of  no  less  than  18°K.  Results  obtained  under  these 
conditions  are  shown  in  Table  2-2,  where  they  are  compared  to  3-Hz  data  taken 
at  20°K  before  and  after  irradiation.  After  the  latter  measurement,  the  quasi - 
static  value  was  again  obtained  at  18°K.  The  data  are  once  again  in  good  agree¬ 
ment,  although  in  this  instance  the  3-Hz  value  of  flatband  shift  tends  uniformly 
to  be  somewhat  smaller  than  the  quasistatic  value.  That  this  is  unrelated  to 
charge  rearrangement  between  18  and  20°K  is  confirmed  by  the  reproducibility 
of  the  quasi  static  measurements  nerformed  both  before  and  after  the  3-Hz  meas¬ 
urement  at  the  higher  temperature  (Table  2-1).  We  may  therefore  conclude  that 
the  results  of  measurement  at  18°K  are  identical  to  those  obtained  at  20°K. 

The  range  of  electron  trap  energy  levels  which  could  be  effective  in  charge 
buildup  is  thus  even  further  restricted. 


Table  2-2.  Additional  comparison  of  3-Hz  and  quasistatic  measurements  of 
flatband  voltage  shift  (dose  =  10^  rads  (Si)). 


Device 

# 

Applied 
Field 
x  106V/cm 

Pre- 

Flatband 

Rad 

Voltage, 

Post 

VF3 

-Irradiation 

AV 

FB 

3-Hz 

(20°K) 

Quasi¬ 

static 

(18°K) 

Quasi - 
static 
(18°K) 

3  Hz 
(20°K) 

Quasi¬ 

static 

(18°K) 

Quasi¬ 

static 

(18°K) 

J-Hz 

(20°K) 

1 

1.61 

-.150 

-.151 

+1.51 

+1.47 

+1.50 

1.66 

i.62 

2 

2.12 

-.120 

-.129 

+1.61 

+1.57 

+  1.60 

1.74 

1.69 

3 

2.68 

-.062 

-.153 

+1.69 

+1.68 

+1.69 

1.84 

1.74 

4 

3.23 

- .  134 

-.190 

+1.43 

+1.41 

+1.43 

1.62 

1.59 
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2  2.3  Direct  Measurements  of  Flatband  Shift  at _4^K 

The  first  experimental  attempt  to  perform  a  4.2°K  experiment  utilized  the 
feet  that  majority  carrier  charge  for  sensing  the  transition  between  accumula¬ 
tion  and  depletion  can  be  made  available  from  the  charge  stored  at  the  inter¬ 
face  under  accumulation  bias  conditions.  The  experimental  method  is  as  follows: 
The  capacitor  is  brought  to  equilibrium  in  accumulation  prior  to  irradiation  by 
momentarily  raising  the  sample  temperature  to  19CK  under  accumulation  bias.  The 
sample  temperature  is  then  lowered  to  4.2°K.  As  the  temperature  is  decreased, 
the  majority  carriers  freeze  out  throughout  the  silicon  bulk,  but  the  accumula¬ 
tion  charge  remains  in  the  conduction  band  in  a  thin  region,  on  the  order  of 
100  }.  or  so,  at  the  surface.  To  establish  the  pre-irradiation  flatband  voltage, 
the  bias  is  swept  from  accumulation  to  inversion.  At  flatband,  all  of  the  major¬ 
ity-carrier  charge  in  the  accumulation  layer  will  have  diffused  into  the  silicon 
bulk,  and  the  resulting  current  transient  can  bo  sensed  in  the  external  circuit 
in  the  usual  configuration  for  quasistatic  capacitance  measurement.  Once  the 
charge  has  dissipated,  however,  it  cannot  be  restored  (in  finite  timel)  at  that 
temperature,  so  the  measurement  is  one-shot.  After  restoration  of  the  appropri¬ 
ate  accumulation  bias,  irradiation  of  the  sample  is  initiated.  The  capacitor 
may  either  be  brought  to  equilibrium  charge  distribution  again  at  19°K  prior  to 
irradiation  at  4°K,  or  the  irradiation  may  proceed  directly  at  the  latter  tem¬ 
perature,  in  that  the  carriers  generated  by  the  irradiation  will  restore  the 
capacitance  to  equilibrium  in  accumulation  promptly.*  After  irradiation,  the 
flatband  voltage  may  again  be  sensed  by  sweeping  the  bias  once  from  the  irradi¬ 
ation  bias  to  depletion. 

★ 

The  statement  that  such  equilibration  does  take  place  under  gamna  irradiation  is 
closely  allied  with  the  earlier  assertion  that  the  applied  bias  is  dropped  en¬ 
tirely  across  the  oxide  even  in  the  carrier  freeze-out  temperature  regime.  That 
this  is  in  fact  the  case  was  demonstrated  by  an  experiment  in  which  the  charge 
transport  through  the  oxide  was  measured  under  bias  at  the  prevailing  gamma  dose 
rate.  A  current  of  2.5  x  10"10.4  was  measured  at  41.3V  bias  (4.26  x  106  V/cm). 

This  implies  an  effective  insulator  resistance  of  1.65  x  10*1  ohms.  Since  the 
observed  current  is  far  larger  than  that  which  could  arise  from  carrier  genera¬ 
tion  within  the  oxide,  it  is  in  fact  dominated  by  currents  generated  within  the 
associated  cabling  and  sample  mounts,  which  are  likewise  exposed  to  the  gamma 
source.  Thus,  the  above  value  of  resistance  represents  an  overestimate  of  the 
induced  conductivity  of  the  oxide.  Turning  now  to  the  silicon,  the  observed 
dielectric  relaxation  frequency  of  120  Hz  indicates  a  resistance  of  10'  ohms. 

The  resistance  ratio  of  greater  than  104  implies  that  the  bias  is  applied  ex¬ 
clusively  across  the  oxide.  A  second  argument  can  be  given  from  the  data  them¬ 
selves.  The  agreement  of  the  bias  dependence  of  flatband  shift  at  4°K  with  that 
at  77°K  is  proof,  a  posteriori,  that  identical  bias  conditions  must  have  existed 
across  the  oxide  in  both  cases. 
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In  an  Initial  measurement  of  the  pre- Irradiation  flatband  voTtage  by 
this  method,  results  were  .Tound  t-*  differ  from  our  expectations  In  several  re¬ 
spects.  In  the  first  place,  the  observed  signal  was  found  to  be  considerably 
smaller  than  expected.  The  current  waveform  was  therefore  dominated  Wrongly 
by  capacitor  displacement  current.  The  capacitor  dielectric  1r,  this  case  con¬ 
sists  not  only  of  the  Si02  but  of  the  silicon  bulk,  yielding  a  total  device 
capacitance  somewhat  less  than  1  pF.  To  compensate  this  displacement  cirrent, 
a  bucking  circuit  was  added,  and  the  resulting  experimental  design  is  illus¬ 
trated  in  Figure  2-8.  Data  obtained  with  this  circuit  are  shown  in  Figure  2-9, 
where  curve  (a)  shows  the  current  waveform  after  the  interface  has  been  brought 
to  equilibrium  in  accumulation,  and  curve  (b)  shows  a  subsequent  trace  in  which 
accumulation  charge  is  absent.  The  difference  between  the  two  curves  is  due  to 
the  accumulation  charge,  and  the  bias  at  which  the  charge  disappears  should 
yield  the  flatband  condition.  The  magnitude  of  the  observed  charge  was  found 
to  be  considerably  less  than  that  which  would  be  observed  if  the  majority  car¬ 
riers  were  to  transit  the  whole  silicon  bulk.  One  possible  explanation  for 
this  is  in  terms  of  trapping  of  the  majority  carriers  in  ionized  donor  sites. 

A  certain  concentration  of  compensating  acceptors  inevitably  exists  in  the  ma¬ 
terial,  and  an  equal  concentration  of  shallow  donors  will  therefore  have  yielded 

their  charge  to  these  compensating  sites.  A  good  estimate  of  such  compensation 

13  -3 

in  Czochrilsk'  material  is  3  x  10  cm  .  The  capture  cross  section  at  ionized 

18 

donor  sites  is  quite  large,  so  that  carriers  are  readily  trapped.  Numeri¬ 
cally,  it  is  found  that  the  charge  released  from  the  interface  is  approximately 
equal  to  the  number  of  ionized  donors,  under  the  above  estimate.  The  mean 
transit  distance  of  carriers  is  therefore  predicted  to  be  on  the  order  of  one- 
half  of  the  silicon  thickness.  The  observed  signal  current  was  small  compared 
to  the  value  expected  on  this  basis,  and  other  mechanisms  must  be  operative. 

It  is  possible,  for  example,  that  mobile  carriers  are  trapped  not  only  at  ion¬ 
ized  donor  sites  but  also  at  neutral  sites,  leading  to  a  shorter  transit  dis¬ 
tance.  The  result  in  any  event  is  the  buildup  of  space  charge  in  the  silicon, 
since  the  detrappirsg  times  are  large  compared  to  experimental  time  scales.  The 
measured  flatband  voltage  is  thus  determined  as  much  by  the  fixed  charge  exist¬ 
ing  within  the  silicon  as  that  existing  within  the  oxide.  It  becomes  clear 
that  this  experiment  requires  more  refined  treatment.  Such  analysis  would  have 
been  performed  but  for  the  availability  of  a  more  direct  experimental  option. 
Under  the  circumstances,  this  experimental  approach  was  abandoned,  and  is 
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Figure  2-8.  Schematic  diagram  of  the  experimental  method  for  measuring 
flatband  voltage  utilizing  accumulation  charge.  The  method  is  es¬ 
sentially  a  quasi  static  measurement  with  a  bucking  circuit  to  com¬ 
pensate  for  the  prevailing  displacement  current  transient. 


Figure  2-9.  Quasistatic  measurement  of  flatband  voltage  at  4°K  using 
accumulation  charge.  Shown  are  results  with  accumulation  charge 
present  (Curve  a)  and  without  such  accumulation  charge  (Curve  b). 
The  difference  in  curve  shapes  is  ascribed  to  the  accumulation 
charge,  and  the  break  in  the  curve  taken  to  yield  the  flatband 
voltage. 
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reported  here  for  the  reason  that  it  does  constitute  a  useful  tool  for  the 
characterization  of  semiconductors  in  the  freeneout  regime  when  suitably  an¬ 
alyzed. 

An  unambiguous  method  of  establish"'  >c  wrcther  electron  trapping  exists 
at  4°K  which  is  unobservable  at  18-20°K  comists  in  measuring  any  charge  trans¬ 
port  which  occurs  under  bias  as  the  sample  temperature  is  raised  to  the  higher 
temperature  after  irradiation.  This  relatively  more  straightforward  method 
was  employed  to  ascertain  that  no  significant  electron  trapping  occurred  in 
these  oxides  at  4.2°K.  The  experimental  method  consisted  in  irradiation  of 
capacitors  under  bias,  with  maintenance  of  the  bias  as  the  source  is  lowered 
(i.e.,  dose  rate  -♦  zero),  and  subsequent  measureuent  of  charge  transport  as 
the  sample  temperature  is  raised  from  4  to  20°K.  The  sample  bias  is  held  con¬ 
stant  throughout.  A  Keithley  616  was  used  in  the  charge  integration  mode  to 
yield  charge  directly;  also,  current  was  measured  using  a  Burr-Brown  3523L 

O 

electrometer  operational  amplifier.  The  total  charge  p  (coul/cm  )  generated 

by  an  ionizing  dose  y  (rads  (Si))  in  the  oxide  car  be  calculated  from  the 
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known  pair  creation  energy  of  18  eV.  The  relationship  is  p  =  1.26  x  10  y. 

At  the  applied  bias  of  about  25V  used  in  this  experiment,  nearly  all  of  this 
charge  would  be  collected  in  the  absence  of  trapping;  i.e.,  the  probability  of 
geminate  recombination  is  small. ^ 

In  the  first  experiment,  four  capacitors  were  connected  in  parallel  to 

enhance  the  signal-to-drift  ratio.  The  devices  were  exposed  to  a  dose  of 

5  x  10  rads(Si).  The  maximum  possible  trapped  electron  charge  therefore  amounts 

to  about  9  x  10  ^C.  If  all  of  this  charge  were  to  be  transported  through  the 

oxide,  the  charge  measurement  in  an  external  circuit  would  be  half  of  this,  or 

4.5  x  10  10 C.  The  amount  of  charge  actually  measured  when  the  temperature  was 

-12 

raised  came  to  about  7  x  10  '  C.  If  this  residual  charge  were  indeed  due  to 
electron  trapping,  the  trapped  fraction  is  clearly  too  small  to  be  of  practical 
significance.  It  is  not  in  fact  possible  to  ascribe  the  observed  charge  ex¬ 
clusively  to  electrons  released  from  traps  during  the  temperature  excursion 
since  it  is  known  that  a  hole  trap  distribution  exists,  and  it  is  possible  that 
some  holes  are  released  from  shallow  traps  as  well.  So  the  measured  charge 
places  an  upper  limit  on  the  possible  electron  trapping,  and  this  limit  is  in¬ 
deed  a  very  low  one. 
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Since  this  experiment  was  critical  to  our  conclusions,  it  was  subsequently 

repeated  with  fresh  devices.  Irradiation  was  to  10*  rads  (Si),  and  the  charge 

-12 

transport  ooserved  on  a  single  device  amounted  to  4  x  10  C,  or  1.8%  of  the 

maximum  possible  charge.  A  measurement  of  another  capacitor  on  the  same  chip 

-12 

after  an  additional  irradiation  yielded  a  value  of  9.4  x  10  C.  Although  the 
quantity  of  charge  measured  is  not  highly  reproducible,  the  conclusion  remains 
firm  that  electron  trapping  on  the  time  scale  of  the  experiment  (20  minutes)  Is 
insignificant.  This  experimental  method  would  be  quite  useful  to  examine  the 
possibility  of  radiation-induced  shallow  electron  traps,  by  extension  of  the 
measurement  to  larger  doses,  but  this  was  not  done. 

An  additional  experimental  method  for  CV-measurement  at  4°K  was  pursued 
in  our  desire  to  make  the  experiment  as  direct  as  possible,  and  the  conclusion 
rigorous.  Advantage  was  taken  of  the  fact  that  the  carrier  concentration  re¬ 
quired  to  perform  low-frequency  CV-measurements  can  be  generated  by  the  ioniz¬ 
ing  radiation  itself.  The  generation  rate  of  the  Co-60  gamma  source  was 

14  1 

3  x  10  carriers/cm  sec  for  our  conditions,  and  with  an  estimated  carrier 
lifetime  of  10"*®  sec  (both  majority  and  minority  carrier  lifetimes  need  to  be 
considered  here),  a  carrier  concentration  of  3  x  10*  cm-3  was  estimated,  lead¬ 
ing  to  the  expectation  of  a  rather  low  dielectric  relaxation  frequency  (about 
7  Hz).  Initial  measurements  were  therefore  made  under  gamma  irradiation  with 
an  excitation  frequency  of  1  Hz  using  the  Ithaco  391A  lock-in.  The  conditions 
of  measurement  were  not  auspicious,  since  gamma  irradiation  of  the  long  cable 
runs  to  the  sample  holder  gave  rise  to  considerable  noise.  The  resulting  CV- 
curves  were  not  promising  for  the  recovery  of  flatband  voltage.  Operation  at 
the  low  frequency  of  1  Hz  is  expected  to  yield  the  typical  low-frequency  ca¬ 
pacitance  curve,  in  which  the  minority  carriers  are  able  to  maintain  equilibri¬ 
um  with  the  modulation  in  inversion,  resulting  in  the  observation  of  the  oxide 
capacitance  both  in  accumulation  and  in  inversion.  Measurements  at  higher  fre¬ 
quencies  were  expected  to  yield  the  high-frequency  capacitance,  and  this  was 
confirmed.  Signal-to-noise  ratios  also  improved  at  higher  frequencies.  To 
establish  the  magnitude  of  the  dielectric  relaxation  frequency  under  the  given 
generation  rate,  measurements  were  then  made  of  the  capacitance  signal  in  ac¬ 
cumulation  as  a  function  of  frequency  using  the  Ithaco  395  lock-in,  and  the 
results  arc  shown  in  Figure  2-10  for  the  frequency  regime  of  5  Hz  to  2  kHz. 

The  fit  to  the  data  indicates  a  dielectric  relaxation  frequency  of  about 
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120  Hz.  Measurements  were  therefore  made  of  CV-curves  using  100  Hz  modulation, 
and  the  curves  obtained  were  free  of  significant  noise. 

Values  of  shifts  in  the  flatband  voltage  at  4°K  obtained  using  the  above 
technique  were  found  to  be  lower  than  those  obtained  using  the  conventional 
3-Hz  CV-measurement  at  20°K.  Comparisons  were  made  for  two  oxide  thicknesses 
and  for  a  number  of  gate  biases.  The  observed  shifts  ranaed  from  0.60  to  0.83 
of  the  corresponding  conventionally  measured  values.  The  possibility  existed 
that  the  measurement  method  involved  a  certain  constant  error,  and  one  irradi- 

4 

ation  was  therefore  performed  to  a  dose  of  4  x  10  rads  to  yield  a  larger  abso¬ 
lute  shift.  The  4°K  measurement  again  yielded  a  value  which  was  about  0.83  of 
the  conventionally  measured  one.  These  results  prompted  the  revival  of  the 
hypothesis  that  some  electron  trapping  might  have  been  present  and  responsible 
for  this  reduction  in  net  charge  buildup.  An  additional  experiment  was  there¬ 
fore  undertaken  to  resolve  this  conflict  with  our  previous  findings.  It  was 
supposed  that  since  the  experimental  method  employed  he^e  is  unconventional, 
some  difficulties  relating  to  the  method  itself  might  have  given  rise  to  the 
discrepancy. 

A  sequence  of  measurements  was  made  post-irradiation  in  which  the  flat- 
band  voltage  was  first  measured  under  garrma  excitation  at  4°K,  then  the  sample 
temperature  was  raised  to  20°K  and  the  flatband  voltage  remeasured  convention¬ 
al  ly.  Finally,  the  temperature  was  once  again  dropped  to  4°K,  and  the  source 
raised  to  remeasure  the  flatband  voltage  at  4°K.  Results  of  this  sequence  of 
tests  are  shown  in  Table  2-3  for  two  runs.  The  flatband  shifts  measured  at 
4°K  are  once  again  smaller  than  those  measured  at  20°K,  but  the  significant 
finding  is  that  these  lower  values  are  reproduced  after  the  sample  has  been 
taken  to  20°K,  and  any  (hypothetically)  trapped  electrons  have  been  released. 
The  explanation  for  the  observation  of  a  lesser  shift  at  4°K  therefore  must  lie 
elsewhere  then  with  electron  trapping,  and  our  earlier  conclusion  stands  con¬ 
firmed.  It  may  be  noted  that  whereas  the  flatband  shift  at  4°K  was  reproduced 
very  well  before  and  after  the  "anneal",  the  reproducibility  <rom  run  to  run 
is  quite  poo*'.  The  randomizing  variable  has  not  been  identified  and  the  issue 
was  not  pursued  since  the  principal  purpo:  s  of  the  experiment  had  been  accom¬ 
plished.  Some  additional  comments  on  th^s  technique  are  in  order,  however. 
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Table  2-3.  Sequential  measurement  of  flatband  shifts  at  4°K  and  20°K. 

First  Run: 

1.  AV^  (4°K,  under  gaimia  excitation,  100  Hz):  1.060V 

2.  (20°K,  thermal  excitation,  100  Hz):  1.275V 

3.  AV^  (4°K,  under  gamma  excitation,  100  Hz):  1.060V 

Second  Run: 

1.  AV^  (4°K,  under  gamma  excitation,  100  Hz):  0.762V 

2.  AV^h  (20°K,  thermal  excitation,  100  Hz):  1.077V 

3.  AV^  (4°K,  under  gamma  excitation,  100  Hz):  0.762V 

It  had  been  clear  from  the  outset  that  this  experimental  approach  was 
fraught  with  pitfalls.  The  gamma  dose  rate  is  clearly  sufficient  to  affect 
the  capacitance- voltage  curve  by  charge-up  of  existing  interface  states.  If 
the  difference  in  observed  flatband  shift  may  be  ascribed  to  such  effects, 
then  it  is  not  unreasonable  to  suggest  that  a  dose  rate  does  exist  which  on 
the  one  hand  permits  CV-measurements  to  be  made  at  tractable  modulation  rates, 
and  on  the  other  hand  produces  only  small  perturbations  of  the  capacitance- 
voltage  relationship.  If  a  dose  rate  on  the  order  of  10%  of  the  prevailing 
dose  rate  had  been  used,  for  example,  an  operating  frequency  of  30  Hz  would 
still  have  been  feasible,  and  the  prompt--as  opposed  to  cumulative--effects  of 
the  gamma  dose  rate  on  the  observed  flatband  condition  would  be  reduced  an 
order  of  magnitude,  to  quite  manageable  proportions.  This  is  particularly 
true  when  it  is  noted  that  one  need  not  operate  with  such  small  flatband 
shifts  as  those  to  which  we  have  restricted  ourselves  in  this  experiment. 
Unfortunately  these  speculations  were  not  tested,  but  are  nevertheless  pre¬ 
sented  as  a  possible  avenue  for  the  examination  of  oxides  under  carrier- 
freezeout  conditions. 

In  an  attempt  to  understand  gamna-ray  effects  on  capacitance-voltage  re¬ 
lationships,  measurements  with  and  without  gamma  excitation  were  compared  at 
various  temperatures,  and  the  results  are  shown  in  Figure  2-11.  In  Figure  2-lla, 
comparison  is  made  of  a  C V  curve  obtained  at  4°K  under  gamma  irradiation  with  a 
20° K  CV  curve  obtained  without  gamma  irradiation.  The  flatband  voltages  in  both 
cases  are  approximately  the  same  in  this  pre-irradiation  measurement.  The  sig¬ 
nal  amplitude  in  accumulation  is  not  the  same  in  the  two  cases,  however,  despite 
the  use  of  an  identical  modulation  amplitude  and  the  use  of  a  modulation  frequency 
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Figure  2-11.  Comparison  of  CV  curves  obtained  at  various  temperatures  under  differing  measurement  conditions 


(100  H2 )  less  than  the  dielectric  relaxation  frequency.  This  could  be  due  to 
voltage  division  between  the  silicon  bulk  resistance  and  the  lock-in  amplifier 
input  impedance,  but  that  possibility  was  not  checked.  When  CV  measurements 
made  at  20°K  with  and  without  gamma  excitation  are  compared  (Figure  2- lib),  a 
broadening  of  the  capacitance  step  is  observed  which  is  indicative  of  charge-up 
of  interface  and/or  bulk  states  that  are  susceptible  to  charge  exchange  on  the 
experimental  time  scales.  The  curve  obtained  under  gamma  flux  has  the  character 
of  a  low-frequency  CV  curve,  so  that  the  inversion  capacitance  is  never  exhibited. 
When  1-MHz  CV  measurements  made  at  77°K  with  and  without  gamma  excitation  are 
compared  (Figure  2-llc),  the  effect  of  the  gamma  flux  is  simply  to  maintain  the 
inversion  layer  in  equilibrium. 

2.2.4  Oxide  Thickness  Dependence  of  Charge  Buildup 

Using  the  conventional  low-frequency  capacitance  method  described  above, 
samples  with  identically  grown  oxides  but  different  oxide  thickness  were  char¬ 
acterized  in  terms  of  flatband  voltage  shift  as  a  function  of  gate  bias.  (Wafers 
were  etched  back  to  obtain  the  different  oxide  thicknesses.)  The  results  are 
shown  in  Figure  2-12.  The  scatter  is  indicative  of  the  experimental  resolution 
under  the  given  conditions  of  long  cable  runs,  etc.  The  thickness  dependence 
of  the  maximum  rates  of  charge  buildup  is  shown  in  Figure  2-13,  where  data 
points  corresponding  to  a  field  strength  of  2.8  x  106  V/cm  are  replotted  from 
Figure  2-12.  The  observed  thickness-squared  dependence  is  consistent  with  uni¬ 
form  charge  buildup  within  the  oxide,  as  is  characteristic  of  higher- temperature 
behavior  as  well  (prior  to  significant  hole  transport). 

The  maximum  charge  buildup  which  is  expected  to  occur  if  all  the  holes 
are  trapped  in  the  oxide  near  their  point  of  generation  and  all  the  electrons 
are  swept  out  may  be  determined  as  follows/  For  uniform  charge  distribution, 
the  shift  in  flatband  voltage  AV,^  is  given  by 

<« 

The  maximum  number  of  free  electron-hole  pairs  induced  by  a  given  dose  is  7.9 

12  3 

x  10  pairs/cm  rad  (employing  18  eV  per  pair).  We  therefore  obtain 

4Vf[)  =  1.8  x  106  v  tox2  .  (2) 
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Figure  2-12.  Flatband  voltage  shift  as  a  function  of  field 
applied  during  irradiation  for  MOS  capacitors  with  var 
ious  i -*ide  thicknesses.  Irradiations  were  performed  at 
4°K  and  CV  measurements  then  made  at  20°K. 


AV-  (volts) 


Figure  2-13.  Flatband  voltage  shift  at  2.4  x  106  V/cm 
as  a  function  of  oxide  thickness  (derived  from 
Figure  2-12). 


For  an  oxide  thickness  of  970 A,  this  yields  a  maximum  possible  shift  of  1.71V  for 

a  dose  of  10  rads(Si).  The  observed  maximum  shifts  of  1.25V  In  Figure  2-12  and 

1.45V  in  Figure  2-2  are  somewhat  less  than  the  maximum  possible  value  due  to  field- 

induced  transport  of  charge  on  the  time  scale  of  the  measurement,  which  is  already 

significant  at  field  strengths  of  3  x  10®  V/cm.  Previous  modeling  of  such  hole 
4  n 

transport  data  at  77  K  also  indicates  that  the  charge  yield  is  less  than  100%  at 
this  field  strength.  The  net  affect  is  that  a  charge  of  about  80%  of  the  maximum 
value  is  expected,  or  a  peak  flatband  shift  of  about  1.4V.  The  data  are  in  reason¬ 
able  agreement  with  this  expectation. 

The  data  of  Figure  2-6  indicate  a  much  larger  flatband  voltage  shift  than 
the  two  figures  discussed  above.  The  peak  value  of  about  1.9V  is  clearly  in  ex¬ 
cess  of  the  maximum  theoretical  shift.  A  scale  factor  error  is  indicated;  since 
the  data  are  otherwise  well-behaved.  The  principal  point  here  is  the  correspon¬ 
dence  between  the  4°K  and  the  770|K  data,  the  latter  being  well  understood,  and 
this  purpose  is  still  served  despite  any  uncertainty  in  scale  factor. 

Dosimetry  fo*'  all  measurements  was  performed  by  means  of  an  ionization  gauge 
calibrated  for  Co-60  gammas.  The  gauge  was  placed  at  the  sample  position  in  the 
absence  of  the  dewar,  and  the  absorption  by  the  dewar  needed  to  be  determined. 

This  was  accomplished  in  two  ways.  First,  the  absorption  was  measured  by  means  of 
TLD's ,  and  determined  to  be  11%.  In  view  of  the  uncertainty  in  TLD  dosimetry, 
which  was  taken  to  be  10%,  this  value  is  not  too  meaningful.  A  second  measurement 
was  made  sing  charge  buildup  itself  as  a  dosimetric  technique.  Charge  buildup  was 
measured  in  the  experimental  samples  at  77°K  at  1  MHz  in  the  helium  dewar,  and  the 
measurements  were  repeated  in  a  glass  dewar,  for  which  the  absorption  was  negli¬ 
gible  (^2%).  By  this  latter  method,  the  dewar  absorption  factor  was  found  to  be 
negligible  as  well,  and  this  latter  measurement  was  accepted  as  the  more  accurate. 

A  second  measurement  was  made  to  ascertain  whether  dose  enhancement  prevailed 
at  the  sample  position.  The  sample  had  originally  been  covered  with  a  thin  (30  mil) 
layer  of  aluminum  to  assure  dose  equilibrium,  but  a  small  hole  had  been  cut  into 
the  shield  to  admit  the  optical  signal  during  the  early  runs.  The  helium  exchange 
gas  chamber  was  made  of  stainless,  so  the  possibility  of  some  dose  enhancement 
could  not  be  ''verlooked.  An  additional  run  was  made,  therefore,  in  which  the  alum¬ 
inum  shield  was  replaced  with  one  without  an  aperture,  and  identical  results  in 
terms  of  flatband  shift  were  obtained. 
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2.3  DISCUSSION  AND  CONCLUSION 
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It  is  the  compelling  conclusion  of  the  present  work  that  charge  buildup 
at  4°K  in  irradiated  SiOg  films  can  be  described  exclusively  in  terms  of  the 
hole  trapping  behavior  which  has  already  been  successfully  modeled  at  higher 
temperatures  (notably  77°K).4’5  Strictly  speaking,  this  conclusion  should  be 
modified  to  apply  specifically  to  the  case  of  the  wet  oxides  studiad  and  to 
the  particular  irradiation  conditions  (steady-state  ionizing  radiation--20 
minute  time  scale).  It  remains  a  possibility  that  significant  electron  trap¬ 
ping  may  be  observed  on  shorter  time  scales,  and  this  is  a  matter  which  should 
be  further  investigated.  It  seems  highly  likely  on  general  principles  that 
trap-limited  electronic  conduction  should  be  observed  on  some  time  scale  at 
these  low  temperatures.  If  this  time  scale  is  sufficiently  long,  it  would 
still  be  of  technological  interest  as  a  means  of  compensating  for  trapped 
holes  on  that  time  scale. 


It  has  not  been  mentioned  thus  far  that  shallow  electron  traps  in  oxides 
have  been  the  subject  of  study  by  others.  Most  recently,  shallow  neutral  elec¬ 
tron  craps  have  been  studied  in  SiO,  films  by  thermally  stimulated  current 
19  ^ 

measurements.  The  distribution  in  energy  of  such  traps  was  found  to  be  about 
300  meV  below  the  conduction  band,  with  an  uncertainty  of  about  50  meV,  and  a 
half-width  of  the  distribution  of  about  200  meV.  These  findings  are  not  in 
conflict  with  those  reported  here.  The  electron  traps  in  question  have  very 
low  capture  probabilities,  and  were  found  under  conditions  of  electron  injec¬ 
tion  which  exceeded  by  three  to  four  orders  of  magnitude  the  charge  generated 
in  the  present  penetrating-irradiation  experiments.  Under  conditions  in  which 
both  electrons  and  holes  are  available  for  trapping,  we  have  found  the  residual 
trapped  electron  charge  after  a  20-minute  low- temperature  irradiation  to  be  a 
small  fraction  (no  more  than  about  4%  in  the  present  oxides)  of  the  trapped 
hole  charge.  As  a  final  pcint,  we  have  raised  the  possibility  that  electron 
trapping  and  transport  in  Si may  be  dependent  on  oxide  processing.  This  pos¬ 
sibility  should  be  considered  in  further  studies. 
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SECTION  3.0 


MECHANISMS  OF  RADIATION  EFFECTS  ON  EXTRINSIC  SILICON  DETECTORS 
3.1  INTRODUCTION 

Extrinsic  silicon  detectors  are  becoming  technologically  important 
because  of  their  suitability  for  the  fabrication  of  monolithic  focal  plane 
arrays.  By  appropriate  choice  of  dopant,  the  desired  spectral  region  can  be 
covered,  ranging  from  the  1-2  pm  range  for  the  deep  gold  level  to  about  30  pm 
for  arsenic  or  boron.  Operation  of  the  detectors  is  confined  to  the  tempera¬ 
ture  regime  in  which  the  carriers  are  substantially  frozen  out,  which  ranges 
from  slightly  below  0°C  for  the  gold  acceptor,  to  about  10°K  for  boron  and 
arsenic. 

Extrinsic  silicon  detectors  make  it  possible  to  take  advantage  of  the 
well-established  silicon  integrated  circuit  fabrication  technology  to  incor¬ 
porate  both  detector  elements  and  associated  readout  electronics  in  a  monolith¬ 
ic  structure.  Intrinsic  detectors,  by  contrast,  generally  require  coupling  of 
the  optically  generated  signal  into  silicon  readout  devices,  such  as  CCD's. 
Progress  is  being  made  in  the  development  of  MIS  technologies  for  InSb  and 
other  candidate  intrinsic  detector  materials,  which  would  permit  monolithic 
focal  plane  arrays  using  these  other  materials.  However,  the  criteria  which 
must  be  satisfied  by  a  charge-coupled  device  suitable  for  the  infrared  are  far 
more  stringent  thai.  those  which  apply  in  the  visible  region,  and  the  problems 
in  that  case  have  been  formidable  enough.  The  reason  for  this  is  simply  that 
the  signal  of  interest  is  generally  a  small  perturbation  on  the  total,  or 
background,  optical  signal.  As  an  example,  in  thermal  imaging,  small  tempera¬ 
ture  differences  must  be  made  discernible,  and  thus  comparably  small  differ¬ 
ences  in  optical  flux  must  be  measurable.  A  1°K  temperature  difference  yields 
a  1.8%  difference  in  photon  ftux  at  a  detector  sensitive  out  to  12  pm,  for  ex¬ 
ample.  If  a  100-tlement  CCD  is  envisioned  to  couple  out  the  detected  signal 
sequentially,  then  the  requirements  on  permissible  carri  v  transfer  ineffi¬ 
ciency  are  manifestly  stringent.  Thus,  there  is  :onsiderable  intercut  in 
extrinsic  silicon  detector  arrays. 

This  section  describes  results  of  an  investigation  of  the  mechanisms  of 
radiation  effects  on  extrinsic  silicon  detectors.  The  work  reported  here 
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constitutes  the  beginning  of  an  effort  that  w<ll  continue  in  a  subsequent 
program  during  1979.  We  first  provide  a  background  discussion  of  radiation 
effects  on  extrinsic  silicon,  followed  by  presentation  and  discussion  of  ex¬ 
perimental  results. 

3.2  BACKGROUND 

Radiation  effects  on  extrinsic  silicon  detectors  are  best  discussed  in 
terms  of  the  various  optical  and  transport  parameters  qoverning  operation  of 
such  devices,  and  we  review  these  briefly.  For  specificity,  the  discussion 
will  be  in  terms  of  zinc  as  a  dopant.  Zinc  introduces  an  acceptor  level  at 
about  0.31  eV  above  the  valence  band,  as  established  by  means  of  carrier  con¬ 
centration  versus  temperature  measurements.  At  temperatures  sufficiently  low, 
the  zinc  levels  are  substantially  unionized,  and  susceptible  to  optical  exci¬ 
tation  of  holes  into  the  valence  band.  This  is  shown  schematically  in  Figure 
3-1.  The  probability  of  optical  absorption  and  resultant  carrier  excitation 
is  referred  to  as  the  quantum  efficiency,  T)(\),  and  is  given  in  terms  of  the 
optical  absorption  coefficient,  c*(a),  and  the  reflection  coefficient,  r(\),  by 
the  equation 


mo  - (1  - ,  <n 

l-r(A)  e“'(,l)d 

where  d  is  the  sample  thickness.  This  equation  is  straightforwardly  derived 
under  the  assumption  of  normal  incidence,  uncoated  surfaces,  and  with  inclusion 
of  multiple  internal  reflections.  The  absorption  coefficient,  in  turn,  is  given 
by 


a{X)  =  cr(X)  NZno  ,  (2) 

where  o{\)  is  the  spectral  photoionization  cross  section  and  NZno  is  the  number 
density  of  unionized  zinc  impurities.  Photoionization  cross  sections  typically 
fall  in  the  range  of  lCf1^  to  10"^  cm^,  and  the  solubility  limits  for  typical 
desirable  dopants  typically  fall  in  the  range  of  10^10^  cm"3,  so  that  ab¬ 
sorption  lengths  at  the  peak  response  wavelength  tend  to  fall  in  the  range  of 
0.1-10  mm  for  practical  impurities.  For  simplicity,  we  may  assume  that  the 
sample  thickness  is  less  than  the  absorption  length,  and  that  the  absorption 
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Distance 


VALENCE  BAND 


Figure  3-1.  Band  diagram  for  the  lowest  acceptor 
state  of  zinc  in  silicon.  (A  second,  deeper 
acceptor  state  exists  when  this  level  is 
filled. ) 
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Is  thus  uniform  throughout  the  sample.  The  quantum  efficiency  then  reduces 
simply  to 


^  f 


T)(>)  =  aNd  .  (3) 

The  optically  excited  carriers,  majority  carriers  in  this  case,  undergo 
transport  under  an  ar plied  field  E  by  an  amount  given  in  the  mean  by  the 
Schubweg,  or  transit  distance,  *j,tE,  where  p  is  the  carrier  mobility  and  t  the 
majority  carrier  lifetime.  If  the  contacts  on  the  sample  are  nonblocking  and  in¬ 
jecting,  then  majority  carriers  swept  out  at  one  electrode  may  be  replaced  at  the 
other,  so  that  the  excited  carriers  contribute  up  to  their  whole  Schubweg  to  the 
transport  even  if  it  exceeds  the  sample  thickness.  A  photoconducti ve  gain  G  is 
thus  defined,  given  by  (|j.TE)/i,  where  l  is  the  separation  Between  electrodes. 

The  quantum  yield  is  then  defined  as  the  ratio  of  the  number  of  mobile  carrier 
transits  across  the  photoconductor  to  the  number  of  incident  photons.  Under 
the  assumption  that  the  sample  thickness  is  considerably  less  than  the  absorp¬ 
tion  length,  and  the  absorption  thus  uniform  throughout  the  sample,  the  quan¬ 
tum  yield  is  given  by 

Quantum  Yield  =  T|G  =  UptE/x  .  (4) 

A  related  quantity  in  terms  of  which  detectors  are  typically  character¬ 
ized  is  the  current  responsivity,  Rj,  which  is  the  device  signal  current  nor¬ 
malized  to  incident  optical  power: 

f^U)  =  qT)(\)G/hv  =  qXTl(A)  G/nc  ,  (5) 

where  h  is  the  Planck's  constant,  c  is  the  velocity  of  light,  and  hv  is  the 
energy  of  the  incident  photon  of  wavelength 

Tc  permit  observation  of  the  photoinduced  carrier  density,  detectors  are 
operated  in  the  carrier  freezeout  regime  in  which  the  thermally  generated  car¬ 
rier  concentration  is  less  than  or  equal  to  that  generated  by  the  background 
photon  flux.  The  theoretical  limit  imposed  on  the  detection  of  an  optical  sig¬ 
nal  superimposed  on  the  optical  background  is  given  by  generation- recombination 
(g-r)  noise  in  the  detector.  This  noise  source  derives  from  the  statistical 
fluctuation  in  generation- recombination  events,  and  the  statistics  are 
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independent  of  whether  the  generation  process  is  predominantly  thermal  or  optical. 
Under  conditions  In  which  the  optically  Induced  carrier  density  dominates  the 
thermally  induced  density,  the  operation  of  the  detector  attains  its  limiting, 
i.e.,  background-limited  (BLIP),  value.  A  concrete  treatment  of  g-r  noise 
must  be  based  on  certain  assumptions  with  regard  to  detector  specifications, 
to  which  attention  is  now  turned. 

In  general,  extrinsic  silicon  detectors  for  thermal  Imaging  in  the  8-14 
nm  regime  employ  one  of  the  shallow  III-V  dopants,  such  as  arsenic.  Residual 
acceptors  in  such  material  act  to  compensate,  or  Ionize,  the  intentional  donor 
impurity,  reducing  the  concentration  available  for  photoionization,  and  de¬ 
creasing  the  carrier  lifetime  since  all  the  ionized  donor  sites  are  recombina¬ 
tion  (or,  preferably,  trapping)  sites  for  excited  majority  carriers.  The  ma¬ 
jority  carrier  lifetime  is  thus  given,  to  good  approximation,  by 


°nvth 


for  the  assumed  case.  Here  cn  is  the  capture  probability  at  ionized  arsenic 
sites,  the  concentration  of  which  is  given  by  NAs+-  is  the  capture  cross 
section,  and  vth  the  thermal  velocity.  As  a  practical  matter,  boron  is  the 
principal  residual  acceptor  impurity  in  such  material,  and  the  ionized  donor 
concentration  is  approximately  equal  to  the  residual  boron  concentration. 

For  the  case  of  a  practical  arsenic-doped  detector,  therefore,  the  car¬ 
rier  concentration  temperature  dependence  may  be  described  in  terms  of  a  two- 
level  model  involving  merely  the  intentional  donor  level  of  concentration  Ng 

and  the  residual  acceptor  level,  of  concentration  N..  The  charge  balance  equa- 
20  H 

tion  yields 


no(no+NA) 

PNc(ND  '  NA  e*p(-Ea/kT) 


Here  Ea  is  the  ionization  energy  of  the  donor  level,  Nc  is  the  conduction  band 
density  of  states,  p  is  the  statistical  weight  of  the  donor  level,  and  nQ  is 
the  equilibrium  carrier  concentration  at  temperature  T.  In  the  freeze-out 
regime,  this  equation  reduces  to  the  simpler  form 
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(8) 


"o-  ■^MP(-E.AT) 

The  carrier  concentration  is  seen  to  be  governed  by  the  ratio  of  donor  to  ac¬ 
ceptor  concentrations 

For  deeper  impurities,  which  are  suitable  for  near-infrared  detection, 
the  above  discussion  needs  to  be  modified  only  slightly.  Referring  again  to 
the  case  of  the  zinc  acceptor,  it  is  clear  that  residual  boron  will  dominate 
the  carrier  concentration  temperature  dependence  unless  it  is  completely  com¬ 
pensated  by  a  donor  (typically  shallow).  Any  excess  compensation  will  serve  to 
ionize  the  zinc  level  partially.  If  the  intentional  dopant  concentration  is  suf 
ficiently  large  compared  to  the  residual  boron  dopant  concentration,  then  a  two- 
level  analysis  will  once  again  be  adequate.  We  have,  in  general, 

ep0(p0  +  nd} 

VnA  -  nD  -  p0>  *  exp('Ea/kT)  '  (9> 

which  reduces  at  low  temperature  to  the  familiar  form 

NvNA  , 

po  '  J„-exp(-Ea/kT)  (10) 


and  at  high  temperature  to 

«7Ka° -  *  e*P(-EaAT>  .  (11) 

where  it  is  assumed  that  >  pQ  »  Np.  Here  Np  is  the  net  compensation,  given 
by  Np  =  NAs  -  Np.  In  the  extreme  n igh-temperature  limit,  we  have,  of  course, 

P0  =  Na,  independent  of  statistical  weight  (strictly  speaking,  pQ  =  -  Np). 

Intrinsic  carrier  concentration  has  been  neglected. 

The  above  equations  reveal  that  a  complete  specification  of  the  carrier 
concentration  as  a  function  of  temperature  yields  the  intentional  dopant  con¬ 
centration,  the  concentration  of  compensating  impurity,  and  the  statistical 
weight.  Hall  measurements  as  a  function  of  temperature  are  employed  to  yield 
this  information,  as  will  oe  discussed  further  below. 
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For  extrinsic  detector  material  which  Is  partically  compensated  (this 

21 

covers  all  cases  of  Interest),  van  Vllet  has  shown  that  the  spectral  inten¬ 
sity  of  the  fluctuations  in  device  current  which  arise  from  the  generation- 
recomblnatior.  process  Is  given  per  unit  frequency  Interval  by 

s<fl 1 '  (1 ' E)  ( 77^7 )  <1?) 

where  I  is  the  quiescent  sample  current,  Nc  is  the  total  number  of  carriers 
undergoing  transport  in  the  sample,  g  is  the  fractional  ionization,  and  u>  the 
angular  frequency.  In  the  freeze-out  regime  §  is  negligible.  The  rms  fluctu¬ 
ation  of  current  in  the  time  domain  is  given  by 

co 

J2  2 

S(f)df  =  (1  -  5)  wi-  .  (13) 

o  0  0 

It  is  apparent  from  these  equations  that  an  explicit  lifetime  dependence  is  re¬ 
vealed  only  in  the  spectral  dependence  of  the  fluctuations,  not  in  the  rms  value 
of  the  fluctuations  in  the  time  domain.  The  equations  are  applicable  to  the 
case  of  thermal  imaging  or  near  infrared  detection,  in  which  the  quiescent  car¬ 
rier  concentration  is  a  well-defined  quantity,  i.e.,  the  fluctuations  are  less 
than  the  mean  value.  For  low-background  applications,  a  different  formulation 
is  appropriate. 

As  a  practical  matter,  almost  all  extrinsic  silicon  detectors  currently 
under  development  are  characterized  by  majority  carrier  lifetimes  which  are 
quite  short,  on  the  order  of  microseconds  at  best.  (This  excludes  detectors 
which  employ  a  doubly-changed  acceptor,  such  as  Zn~,  for  pnotoexcitation) .  The 
frequency  domain  of  interest  is  that  given  by  frame  times  in  a  practical  case,  or 
by  circuit  RC  limitations  in  the  present  experimental  context.  In  either  event, 
the  domain  of  interest  is  less  than  5  x  104Hz,  for  which  the  condition  out  «  1 
is  well  satisfied.  In  this  regime,  the  noise  current  in  a  bandwidth  Af  about  fre¬ 
quency  f  is  given  by 

ij!(f)  =  S(f)Af  =  -J— Af  .  (14) 

o 

Substituting  for  the  current  in  terms  of  Nq  and  photoconducti ve  gain  G,  we  have 
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(15) 


1 


n,ms 


The  quiescent  carrier  concentration  is  comprised  of  a  thermally  generated  com¬ 
ponent  nth  and  an  optically  generated  component  n$,  in  terms  of  which  we  have 


i 


n,rms 


(16) 


where  V  is  the  detector  volume.  It  is  seen  to  be  physically  Immaterial  whether 
the  noise  arises  from  thermally  or  optically  generated  carriers. 


The  limits  of  detection  which  are  imposed  by  material  parameters  are  gen¬ 
erally  expressed  in  terms  of  detectivity  (D),  the  reciprocal  of  the  noise  equiv¬ 
alent  power  (NEP).  The  latter  quantity  is  simply 


NEP  = 


^n.rms^I 


» 


where  Rj  is  the  responsivity,  given  by  Equation  (5).  The  specific  detectivity, 
or  D*.  is  defined  as  the  detectivity  normalized  to  unit  area  of  illumination  and 
unit  bandwidth,  giving 

D*  =  (NEP)"1  (A)*5  (Lf)h  =  RI(AAf)J'/in>rms  .  (17) 


Thus, 


0*  = 


(18) 


Uniform  absorption  over  the  sample  volume  is  assumed.  The  absence  of  photo- 
conductive  gain  from  this  equation  should  be  noted.  This  signifies  that  the 
limits  of  detection  are  influenced  by  photoconductive  gain  only  insofar  as 
amplifier  noise  is  important.  Under  conditions  of  ng  »  n^,  and  writing  ns 
as  T,(A.)0(\)T/d,  where  tf>(X)  is  the  photon  flux  density,  we  obtain 


D* 


.  -X.. 
2hc 


(19) 


This  basic  equation  has  been  developed  here  on  the  one  hand  because  its 
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elegant  simplicity  deserves  once  again  to  be  appreciated,  and  on  the  other 
hand,  because  a  number  of  treatments  of  it  in  the  literature  are  not  without 
error,  or  start  with  assumptions  that  are  no  more  self-evident  than  the  result. 
The  detectivity  is  seen  to  depend  only  on  the  wavelength,  the  quantum  effici¬ 
ency,  and  the  background  flux  density.  The  wavelength  dependence  stems  simply 
from  the  normalization  to  incident  power. 

For  various  applications,  figures  of  merit  for  detectors  other  than  the 
detectivity  are  more  appropriate.  It  is  generally  more  relevant  to  define  the 
detection  limits  in  terms  of  minimum  resolvable  temperature  difference  or  the 
noise  equivalent  flux  density.  The  appropriate  figures  of  merit  for  extrinsic 
detectors  have  been  treated  by  Nelson*^  and  by  LoVecchio^  for  the  near- I R 
applications.  However,  from  a  basic  mechanisms  standpoint,  D*  remains  a  use¬ 
ful  figure  of  merit.  Peak  quantum  efficiencies  are  typically  on  the  order  of 
0.1.  Unity  quantum  efficiency  is  rarely  achieved  due  to  limited  solubilities 
and  low  photoionization  cross-sections.  If  quantum  efficiencies  are  much  less 
than  0.1s  the  detector  is  generally  not  of  practical  use. 

The  background  flux  density  is  usually  not  confined  to  the  wavelength  or 
interest,  as  assumed  in  Equation  (19),  but  is  broadband  and  governed  princi¬ 
pally  by  the  cutoff  wavelength  of  the  detector,  Xc,  given  by  the  energy  con¬ 
servation  criterion  as 


-  1-24  tun 
c  "  Ea|eVT 


(20) 


For  the  detectivity  at  the  peak  response  wavelength  (Xp)  under  conditions  of 
broad  background,  for  example,  we  would  have 

D~(XJ  =  [  . — )  y2hc  ,  (21) 


*  /  V 

,(v*U£) 


where  the  effective  background  flux  density  is  given  by 


■  *  \\ 


(X)o(X)  dX 


(22) 
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Additional  insight  may  be  gained  by  considering  the  dependence  of  nr  and 
n^,  respectively,  on  material  parameters.  We  may  write  (for  donor-type  dopants) 


ns  =  ND0B,eff°pT 


nNA 


0B(^)o(X)  dX 


(.24) 


where  it  has  been  assumed  that  t  =  (c^)"*.  Comparison  of  this  equation  with 
Equation  (8)  for  nth  reveals  that  both  quantities  depend  on  the  compensation 
ratio,  Nq/Na.  From  the  standpoint  of  detector  noise  alone,  then,  it  is  found 
that  the  temperature  at  which  background-limited  conditions  of  operation  ere 
reached  under  a  given  background  is  independent  of  the  degree  of  compensation, 
and  this  latter  parameter  is  again  seen  to  be  significant  only  with  respect  to 
amplifier  noise  limitations.  With  respect  to  radiation  effects  on  detectors, 
this  implies  a  considerably  different  sensitivity  to  detector  parameters  for 
the  case  of  medium-background  thermal  imaging  or  near-infrared  detection  com¬ 
pared  to  low-background  conditions  of  operation,  where  amplifier  noise  is  more 
likely  to  be  dominant.  With  this  brief  introduction,  the  effects  of  neutron  ir¬ 
radiation  on  extrinsic  detectors  are  now  briefly  discussed  in  terms  of  their 
permanent  effects. 

is  useful  to  distinguish  between  the  effects  of  neutrons  on  extrinsic 
decectors  which  affect  the  performance  permanently,  or  cumulatively,  and  those 
effects  which  are  transient.  The  relevant  time  scale  for  the  former  is  that  of 
a  mission,  for  example.  The  time  scale  of  interest  for  the  latter  is  on  the 
order  of  frame  times.  A  typical  detector  array  will  be  employed  with  compensa¬ 
tion  for  quiescent  background,  or  dark,  signals  and  with  normalization  of  re- 
sponsivity  differences  between  detector  elements.  The  relevant  time  scale  in 
that  event  is  that  with  which  the  reference  values  for  the  detector  elements 
are  updated  in  memory. 

The  detector  material  parameters  which  undergo  permanent  change  under 
neutron  radiation  are  those  of  mobility,  majority  carrier  lifetime,  and  dark 
resistivity.  The  lifetime  and  mobility  impact  the  photoconductive  gain  directly. 
The  mobility  is  governed  by  ionized  impurity  scattering,  neutral  impurity 
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scattering,  and  phonon  scattering.  Neutron  irradiation  introduces  both  donor 
and  acceptor  levels  deep  within  the  gap,  and  these  serve  as  additional  compen¬ 
sation  in  both  resistivity  types,  thus  increasing  the  concentration  of  charged 
centers  which  thus  contribute  to  the  most  effective  scattering  mechanism  at  low 
temperature. 

The  majority  carrier  lifetime  is  determined  by  all  possible  processes 
which  remove  the  carrier  from  the  band,  and  is  more  appropriately  referred  to 
as  a  trapping  time,  since  recombination  across  the  bandgap  is  generally  not  in¬ 
volved.  For  thermal  imaging  detectors  utilizing  the  shallow  III-V  dopants, 

such  trapping  occurs  principally  at  the  ionized  shallow  levels,  via  the  Lax 

18  - 

giant-trap,  or  phonon-cascade,  mechanism.  Trapping  in  excited  states  of 

neutral  impurities  appears  also  to  be  important  at  low  temperatures  (on  the 

order  of  10°K)  where  the  probability  of  thermal  reexcitation  from  such  levels 

24 

becomes  negligible.  For  the  case  of  near-IR  detection,  where  deeper  dopants 
are  involved,  the  importance  of  trapping  at  excited  states  of  neutral  centers 
has  not,  to  our  knowledge,  been  definitively  established.  It  would  appear  that 
at  the  higher  temperature  of  operation  characteristic  of  those  detectors,  ther¬ 
mal  detrapping  from  such  shallow  excited  states  would  render  such  trapping  pro¬ 
cesses  unimportant. 

The  dark  resistivity  of  detectors  may  also  be  altered  by  neutron  irradia¬ 
tion.  In  the  case  of  the  shallow  dopants,  any  neutron-induced  level  will  be 
deeper  than  the  intentional  dopant  level.  For  the  case  of  arsenic,  as  an  ex¬ 
ample,  any  such  levels  which  are  acceptor- type  will  be  compensating.  An  exper¬ 
imental  determination  of  such  compensation  has  already  been  made,  and  a  carrier 
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removal  race  of  about  16  per  neutron/cm  was  determined.  The  effect  of  this 
is  to  increase  the  device  impedance  at  a  given  temperature,  according  to  Equa¬ 
tion  (9),  and  to  decrease  the  lifetime.  The  impedance  change  in  this  case  is 
innocuous.  The  impact  on  lifetime  depends  on  the  pre-irradiation  compensation 

level.  If,  for  example,  the  detector  were  characterized  initially  by  a  compen- 

13  -3  12  2 

sation  of  1.6  x  10  cm  boron,  then  a  fluence  of  10  n/cm  would  result  in 

a  reduction  in  lifetime  by  a  factor  of  two.  In  the  case  of  deep  impurities, 

matters  may  be  considerably  different.  The  neutron-induced  level  spectrum  of 

acceptors  and  donors  may  extend  close1"  to  the  bandedge  than  the  intentional 

dopant  level.  In  that  event,  the  net  carrier  removal  rate  may  be  considerably 

reduced  below  that  which  is  observed  for  the  shallow  dopants.  In  the  extreme 
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case,  however,  it  is  possible  that  neutron- induced  levels  which  are  shallower 
than  the  Intentional  dopant  will  govern  the  resistivity.  In  that  event,  lower 
operating  temperatures  would  be  required  to  reach  the  BLIP  limit  of  carrier  con¬ 
centration  than  those  implied  by  the  Impurity  activation  energy  (Equations  (8) 
and  (10)).  Since  our  interest  extends  to  the  deeper  impurities  as  well  as  the 
shallow  ones,  this  possibility  has  been  experimentally  examined,  and  results 
are  given  below.  An  additional  reason  for  examining  the  permanent  neutron  ir¬ 
radiation  effects  for  each  dopant  which  might  be  employed  in  extrinsic  detection 
is  that  the  possibility  of  complex  formation  between  the  Intentional  dopants  and 
neutron- induced  defects  may  be  revealed. 

3.3  EXPERIMENTAL  RESULTS 

The  first  experimental  task  undertaken  consisted  of  the  determination  of 
neutron  effects  on  silicon  doped  with  deep-level  dopants  as  a  complement  to 
measurements  previously  performed  to  determine  the  carrier  removal  rate  in  ar- 

oc 

senic-doped  silicon.  The  dopants  studied  were  zinc,  sulfur,  and  gold.  The 
case  of  zinc,  which  has  already  been  mentioned  (see  Figure  3-1),  is  that  of  a 
relatively  deep  acceptor.  The  sulfur  level  of  interest  is  the  donor  level  at 
about  0.19  eV  from  the  conduction  band,  and  the  gold  level  of  Interest  is  the 
acceptor  level  near  midgap.  The  case  of  gold  is  more  of  interest  from  a  theo¬ 
retical  standpoint  rather  than  as  a  practical  detector  material. 

The  carrier  removal  rate  has  been  determined  in  these  materials  by  Hall 
measurement  of  the  carrier  concentration  as  a  function  of  temperature  In  the 
freeze-out  regime.  This  requires  measurement  under  conditions  in  which  the 
sample  Impedance  is  quite  high,  and  for  that  purpose  the  conventional  Hall 
measurement,  method  has  been  suitably  modified.  The  apparatus  is  shown  In 
Figure  3-2.  Pog-eared  samples  are  used,  in  preference  to  Van  der  Pauw  geome¬ 
try,  because  the  field-strength  in  the  active  region  is  more  uniform,  and  un¬ 
ambiguously  known.  This  is  of  potential  significance  for  responsivity  measure¬ 
ments  since  majority  carrier  lifetimes  may  exhibit  field  dependences.  The  Hall 

and  resistivity  voltages  are  sampled  Dy  means  of  electrometer  operational  ampli- 
13 

fiers  o+'  10  -ohm  input  impedance.  These  amplifiers  are  maintained  under  con¬ 
stant-temperature  conditions  within  the  cryostat,  minimizing  input  capacitance, 
microphonics,  and  pick-up.  Measurement  of  the  carrier  concentration  as  a  func¬ 
tion  of  temperature  ever  the  range  from  that  at  which  carriers  are  totally 
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ionized  to  that  at  which  they  are  totally  frozen  out  permits  determination  of 
the  electrically  active  fraction  of  the  intentional  dopant,  the  concentration 
of  compensating  species,  and  the  statistical  weight,  by  fitting  of  Equations 
(7)  or  (9)  to  the  data.  In  general,  such  determinations  are  subject  to  con¬ 
siderable  uncertainties  with  regard  to  the  compensation  level  and  the  statisti¬ 
cal  weight.  Typical  data  taken  on  a  zinc-doped  sample  are  shown  in  Figure  3-3, 
and  are  useful  to  illustrate  the  point.  The  figure  shows  carrier  concentration 
over  the  range  of  about  10  to  more  than  10  cm"  .  The  low- temperature  car¬ 
rier  concentration  is  determined  by  both  the  compensation  and  the  statistical 
weight.  The  steepness  of  the  curve  in  this  region,  however,  makes  manifest  the 
experimental  difficulty  of  obtaining  precise  values.  The  curve-fit  Is  sensitive 

to  the  choice  of  statistical  weight  only  in  the  intermediate  regime,  i.e.,  where 

-3  -1 

fractional  ionization  ranges  between  10  and  10  .  A  fit  to  this  curve  was 

16 

obtained  with  a  choice  of  statistical  weight  of  unity,  and  the  values  N7  =  10 
-3  14  -3 

cm  and  Nn  =  10  cm  were  determined.  With  a  choice  of  statistical  weight  of 

1/  1 C 

two,  a  slightly  better  fit  was  obtained,  with  the  parameters  N,  =  5.6  x  10 

_  7  1C  _-1  If 

cm  ~  and  Nn  =  1.6  x  1013  cm  .  With  a  choice  of  p  -  4,  we  find  N7  =  7.5  x 
IQ  ,  and  Np  =  8.8  x  10  .  An  additional  variable  in  the  fit  is,  of  course, 
the  activation  energy.  A  value  of  0.31  eV  was  employed  in  all  cases  above.  It 
is  apparent  that  the  availability  of  data  over  13  orders  of  magnitude  does  not 
assure  the  absence  of  ambiguities  in  the  fit,  most  particularly  with  regard  to 
the  very  important  compensation  level.  Since  the  experimentally  observed  quan¬ 
tity  is  the  fractional  change  in  compensation,  our  purpose  is  served  only  if 
the  initial  compensation  level  is  well-known  in  our  specimens. 

The  difficulties  of  establishing  an  accurate  value  of  compensation  level 
are  not  entirely  experimental.  The  fit  presupposes  that  o..(y  a  single  level  is 
important  in  the  determination  of  the  carrier  concentration,  in  addition  to  any 
levels  which  are  compensating.  The  existence  of  excited  states  of  the  impurity, 
and  of  other  impurity  levels  in  lower  concentrations,  impact  this  simple  pic¬ 
ture.  For  this  reason,  it  is  preferred  that  bounds  may  be  placed  on  the  per- 
missible  values  of  the  compensation  level  and  the  statistical  weight  Initially. 
This  is  done  by  studying  material  over  a  rrange  of  initial  resistivities.  If 
these  resistivities  are  aufficiently  low,  then  it  may  be  assumed  that  compen¬ 
sation  levels  will  not  undergo  significant  change  in  the  process  of  diffusion 
of  the  Intentional  impurity,  unless  complexes  form  with  that  impurity.  In  this 
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manner,  a  statistical  weight  of  two  had  previously  been  determined  for  zinc. 

This  approach  was  taken  in  the  present  work,  and  the  material  studied  was  there¬ 
fore  heavily  compensated,  and  not  typical  of  detector- grade  material . 

Results  of  measurements  of  the  temperature  dependence  for  a  zinc-doped 
specimen,  both  before  and  after  neutron  irradiation,  are  shown  In  Figure  3-4. 

The  zinc  concentration  had  been  determined  to  be  9  x  10  cm’'"'  for  this  materi¬ 

al  by  previous  measurements  of  carrier  concentration  at  high  temperature.  The 

pre-diffusion  carrier  concentration  was  determined  from  resistivity  measure- 

15  -3 

msnts  to  be  2.3  x  10  cm  of  phosphorus.  The  irradiation  was  performed  at 
dry  ice  temperature,  and  the  sample  kept  at  that  temperature  except  for  one 
hour  at  room  temperature  just  prior  to  measurement. 

The  observed  changes  in  sample  resistivity  under  irradiation  with  up  to 

6  x  1014  n/crn^  are  seen  to  be  quite  small.  If  the  pre-Zn-diffusion  carrier 

concentration  is  taken  to  be  the  compensation  level,  a  carrier  removal  rate  of 
2 

about  2  cm  /n  is  indicated  for  the  first  irradiation,  and  a  value  of  about 
2 

5  cm'/n  for  the  second  irradiation.  Since  the  observed  change  is  not  large 
compared  to  our  experimental  resolution,  there  remains  some  uncertainty  in 
these  values.  Nevertheless,  it  appears  that  the  carrier  removal  rate  is  con¬ 
siderably  smaller  for  this  deeper  level  than  for  a  shallow  level  like  arsenic. 

A  better  comparison  would  of  course  be  provided  by  boron,  but  we  know  of  no 
such  data. 

Data  for  neutron-irradiation-induced  compensation  of  the  sulfur  donor 
level  at  0.19  eV  below  the  conduction  band  is  shown  in  Figure  3-5  for  two  ir¬ 
radiations.  The  sulfur  concentration  in  this  sample  had  been  determined  to  be 
15 

8.5  x  10  by  resistivity  measurements  at  room  temperature,  where  the  sulfur 

level  is  completely  ionized,  and  the  pre-diffusion  carrier  concentration  had 

15  -3 

been  determined  to  be  1.4  x  10  cm  of  boron,  likewise  from  resistivity  meas¬ 
urements.  Conditions  of  irradiation  were  as  before.  Once  again,  resistivity 

changes  are  found  to  be  small,  and  the  calculated  values  of  compensation  rate 

2 

thus  subject  to  some  error.  The  value  determined  for  the  first  run  is  28  cm  /n 

2 

and  for  the  second,  9  cm  /n.  The  second  value  is  clearly  the  more  reliable. 

2 

It  is  roughly  comparable,  though  somewhat  less,  than  the  value  of  16  cm  /n 

25 

found  for  the  shallow  arsenic  level  by  Kalma.  If  the  difference  is  indeed 
real,  then  it  implies  that  the  effective  acceptor  concentration  is  less  for 
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Figure  3-4.  Results  of  measurements  of  sample  resistance 
as  a  function  of  temperature  before  and  after  neutron 
irradiation  for  a  zinc-diffused  sample.  The  zinc  con¬ 
centration  is  9  x  1016  cm"'*  and  the  phosphorus  con¬ 
centration  is  2.3  x  1015  cm"3. 
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Figure  3-5.  Results  of  measurements  of  sample  resistance  as  a 
function  of  temperature  for  a  sulfur-diffused  sample  be¬ 
fore  and  after  neutron  irradiation.  The  sulfur  concentra¬ 
tion  is  8.5  x  10*5  cm“3  and  -  2.3  x  10l5  cm~3. 
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sulfur  than  for  arsenic,  l.e.  a  fraction  of  the  neutron-induced  acceptor  levels 
lie  shallower  than  the  sulfur  level. 

A  more  extreme  case  of  a  deep  level  Is  represented  by  the  gold  acceptor 
state,  which  lies  near  midgap.  The  gold  level  structure  is  illustrated  in 
Figure  3-6.  A  complication  arises  in  this  case  because  of  the  existence  of  a 
tri-gold  shallow  acceptor,  in  a  concentration  proportional  to  (N^)  .  Suffi¬ 
cient  donors  have  to  be  available  to  compensate  this  shallow  level.  A  gold- 
diffused  sample  was  prepared  in  which  these  conditions  were  met,  and  the  level 

governing  the  resistivity  temperature  dependence  was  therefore  the  deep  accep- 

16  -3 

tor  state.  The  gold  concentration  was  determined  to  be  8  x  10  cm  . 

Results  of  measurements  of  the  carrier  concentration  as  a  function  of 
temperature  before  and  after  neutron  irradiation  are  shown  in  Figure  3-7.  The 
effective  activation  energy  governing  the  temperature  dependence  is  seen  to  be 
reduced  with  fluence.  Since  the  apparent  low- temperature  activation  energy  is 
different  after  the  two  runs,  it  is  probably  not  reasonable  to  associate  it 
with  a  particular  level .  If  such  material  were  to  be  used  in  a  detector  ap¬ 
plication  in  a  neutron  environment,  the  operating  temperature  would  have  to  be 
chosen  appropriate  to  the  observed  unpinning  of  the  Fermi  level  from  the  deep 
gold  acceptor  level. 

3.4  DISCUSSION  AND  SUMMARY 

The  essential  finding  of  the  first  experiments  undertaken  in  this  study 
of  extrinsic  silicon  detector  materials  is  that  the  effective  carrier  removal 
rate,  or  equivalently  the  introduction  rate  cf  compensating  defects  or  defect 
complexes,  depends  strongly  on  the  deoth  within  the  gap  of  the  intentional 
dopant  level.  The  neutron-irradiation-produced  level  spectrum  is  distributed 
throughout  the  bandgap,  so  that  essentially  all  of  the  defects  of  the  appropri¬ 
ate  type  (acceptor  or  donor)  are  compensating  with  respect  to  shallow  dopants, 
and  only  a  fraction  is  effective  for  the  deeper  levels. 

Specifically,  the  evidence  thus  far  accumulated  is  to  the  effect  that  a 
donor  level  at  0.19  eV  from  the  conduction  band  (sulfur)  is  compensated  nearly  as 
effectively  as  a  shallow  donor  (arsenic).  Thus,  few  acceptor  states  are  intro¬ 
duced  within  the  range  of  0.045  to  0.15  eV  from  the  conduction  band,  a  result 
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Figure  3-6.  Level  diagram  for  gold  in  silicon,  showing  the  deep  acceptor 
level,  the  deep  donor  level,  and  the  tri-gold  shallow  acceptor  state. 
Values  for  respective  concentrations  for  the  sample  employed  in  this 
study  are  indicated. 
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Figure  3-7.  Results  of  measurements  of  carrier  concentration 
as  a  function  of  temperature  for  a  gold-diffused  sample 
before  and  after  neutron  irradiation.  The  gold  concentra 
tion  is  8  x  10 16  cm- 3  and  the  phosphorus  compensation  is 
1  x  10^6  cm~3. 


which  Is  not  too  surprising.  The  carrier  removal  rate  Is  found  to  be 
considerably  smaller  for  the  case  of  an  acceptor  level  at  0.31  eV  from  the 
valence  band  (zinc).  It  would  be  helpful  in  the  Interpretation  of  these 
data  to  have  available  the  comparable  data  for  a  shallow  acceptor,  which  will 
hopefully  be  acquired  later  in  the  program.  Under  the  circumstances,  one  can 
say  only  that  either  the  donor  introduction  rate  is  considerably  less  than 
the  acceptor  Introduction  rate,  or  that  the  neui**on- induced  donor  distribution 
extends  to  levels  shallower  than  0.31  eV  in  considerable  measure. 

The  results  for  the  gold  level  near  midgap  indicate  that  the  net  effect 
of  both  the  donor  and  acceptor  distributions  (in  energy)  introduced  by  neu¬ 
trons  is  to  unpin  the  Fermi  level  from  midgap.  The  Fermi  level  appears  to  be 
moved  toward  the  conduction  band,  as  may  be  seen  by  the  following  argument. 

The  conductivity  of  the  specimen  employed  was  n-type,  even  though  the  domi¬ 
nant  level  was  the  acceptor  level  near  midgap.  This  is  consistent  with  calcu¬ 
lations  of  resistivity  at.  room  temperature  for  gold-  and  phosphorus-doped  sili¬ 
con  bv  Thurber  et  al.^®  There  it  is  shown  (following  Bullis^)  that  as  the 
gold  concentration  is  increased,  for  a  given  phosphorus  concentration,  the 
resistivity  goes  through  a  maximum  near  which  the  conductivity  converts  to 
p-type  (by  virtue  of  the  tri-gold  acceptor).  In  the  present  case,  the  room 
temperature  resistivity  remains  approximately  the  same  during  the  irradiation, 
and  the  specimen  remains  n-type.  A  shift  in  the  Fermi  level  toward  the  valence 
band  would  have  led  to  an  increase  in  resistivity  before  a  decrease  would  have 
been  observed  at  higher  fluences.  A  decrease  in  resistivity  to  pre-irradiation 
values  would  also  have  been  accompanied  by  type  conversion,  which  was  not  ob¬ 
served. 

There  are  other  methods  of  obtaining  information  about  the  distribution 
of  neutron- induced  levels  within  the  energy  gap.  These  measure  changes  in  oc¬ 
cupancy  of  deep  traps,  as  in  the  depleted  region  of  a  p-n  junction,  to  indicate 
trap  properties  such  as  energy  level  position  and  capture  rates,  as  well  as 
trap  concentrations.  These  techn'gues  have  the  drawback  that  they  are  not 
equally  sensitive  experimentally  to  all  types  of  traps.  Typically,  such 
methods  are  sensitive  to  traps  in  only  one-half  of  the  bandgap,  making  it 
necessary  to  study  the  other  ha1 f  in  a  sample  which  is  fabricated  differently, 
and  thus  not  directly  comparable.  It  is  therefore  usefu",  to  have  available 
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information  such  as  obtained  here.  Another  point  of  difference  is  that  other 
techniques  look  at  the  dominant  traps,  whereas  in  the  present  experiment,  tails 
of  distributions  may  have  profound  effects. 

The  simplest  model  which  accounts  for  the  fact  that  neutron-irradiation 
is  compensating  in  its  effect  on  both  n-type  and  p-type  (shallow-dopant)  ma¬ 
terial  is  that  of  an  acceptor  level  slightly  above  midgap,  and  a  donor  level 
30 

slightly  below.  The  acceptor  states  would  thus  always  be  available  to  com¬ 
pensate  shallow  donors,  and  vice  versa.  By  virtue  of  the  complex  nature  of 
neutron  damage,  such  levels  are  likely  to  be  distributed  somewhat  in  energy. 

If,  then,  the  donor  distribution  extends  significantly  above  the  midgap  gold 
level,  the  effect  would  be  to  unpin  the  Fermi  level.  An  alternative  explana¬ 
tion  postulates  the  existence  of  a  complex  involving  the  gold  center,  ir  which 
case  our  conclusions  about  the  distribution  of  neutron-induced  levels  would 
have  less  than  general  validity. 

Having  dwelt  on  the  interpretation  of  the  data  for  the  gold  level,  we 
should  recall  that  this  is  unlikely  to  be  considered  a  practical  detector  ma¬ 
terial,  and  that  our  interest  here  is  in  forming  a  coherent  picture.  The  data 
taken  on  the  other  samples,  however,  may  have  more  re1  /ance  to  actual  detec¬ 
tors.  In  practical  detectors,  levels  of  compensation  are  typically  much  smaller 
than  those  characterizing  the  specimens  employed  here  (~10^  cm"3).  Detectors 
for  low  background  applications  must  operate  with  the  minimum  achievable  com¬ 
pensation  consistent  wich  uniformity  requirements.  For  high  background  appli¬ 
cations,  matters  are  different  in  that  one  may  wish  to  design  *or  a  specific 
value  of  compensation,  or  dark  resistance,  not  necessarily  the  minimum.  Even 

1  O  O 

in  that  event,  however,  compensation  levels  would  be  much  closer  to  10  cm 
than  to  10  cm  .  If  a  value  of  1013  cm"3  is  assumed  to  prevail,  then  at  a 

p 

carrier  removal  rate  of  10  cm  /n,  the  responsivity  of  the  detector  would  de- 

12  2 

grade  by  a  factor  of  two  under  a  fluence  of  10  n/cm'.  The  fact  that  the 
given  specimens  did  not  appear  to  undergo  significant  change  under  much  higher 
neutron  fluences  than  this  is  strictly  due  to  the  manner  in  which  we  chose  to 
do  the  experiment,  for  the  reason  discussed  earlier. 

It  is  noteworthy  that  the  carrier  removal  rates  measured  here  for  sulfur 

25 

and  previously  for  arsenic  are  much  higher  than  the  carrier  removal  rates 
measured  at  room  temperature  for  shallow  dopant  material.  The  latter  were 
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measured  to  be  between  about  1.5-3  for  n-type  material,  and  about  3  for  p-type 
material,  using  the  same  irradiation  facility  employed  in  the  present  measure- 
ments.  Room  temperature  irradiations  were  employed,  whereas  dry  ice  irradi¬ 
ation  was  used  here,  and  10°K  irradiations  were  used  for  arsenic.  The  data  for 
sulfur  were  analyzed  without  regard  for  mobility  chanties  under  irradiation. 

These  could  have  modest  impact  on  the  values  determined  for  not  carrier 
removal  rate.  In  subsequent  experiments,  it  is  anticipated  that  better  -alues 
will  be  obtained  for  carrier  removal  rate,  since  the  carrier  concentration  and 
the  mobility  will  both  be  measured,  and  since  rasponsivity  measurements  will  be 
made  that  will  yield  corroborative  information  about  the  rate  of  introduction 
of  compensating  defects.  The  measurements  will  also  be  made  upci  irradiation 
at  detector  operating  temperature,  so  that  immediate  comparison  'f  pre-  and 
post-irradiation  values  will  be  possible  without  an  intervening  change  of 
sample  temperature. 

The  above  data  will  be  taken  as  a  matter  of  course  on  all  the  detector 
materials  to  be  studied.  However,  the  priirarv  focus  of  subsequent  studies  will 
be  the  response  to  transient  illumination,  rather  than  permanent  damage  effects. 
Regarding  the  ’atter,  the  present  work  has  achieved  its  basic  purpose  in  that 
there  appears  to  be  no  profound  consequence  of  neutron  irradiation  on  near- 
infrared  extrinsic  detectors  employing  relatively  deep  dopants  that  is  nAt  ex¬ 
plainable  in  terms  of  the  introduction  of  compensating  species.  The  exception 
is  the  gold  acceptor  level  lying  near  the  midgap,  where  an  unpinning  of  the 
Fermi  level  is  observed. 

The  interest  in  pursuing  studies  of  response  to  transient  extrinsic  illum¬ 
ination  of  the  deep  impurities  in  addition  to  that  of  the  shallow  dopants  stems 
from  the  fav.t  that  neutron-induced  defects  may  result  in  profound  changes  in  the 
recombination,  or  trapping,  properties  of  detectors.  Any  levels  which  are  shal¬ 
lower  than  the  intentional  dopant  level  will  be  subject  to  much  larger  probabili¬ 
ties  of  thermal  reexcitation  of  trapped  carriers  than  would  hold  for  the  latter. 
Such  trapped  carriers  may  therefore  be  read  out  during  subsequent  frame  times, 
leading  to  contrast  reduction  and  image  retention.  This  may  be  thought  of  in 
terms  of  "temporal  cross-talk",  by  analogy  to  spatial  cross-talk  in  a  two-dimen¬ 
sional  array.  In  the  case  of  shallow  dopants,  this  problem  does  not  arise,  since 
bandwidth-1 imitinq  of  the  optical  input  signal  can  be  used  to  insure  that  the 
deeper  levels  are  never  excited. 
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SECTION  4.0 


MECHANISMS  OF  RADIATION  EFFECTS  ON  CHARGE-COUPLED  DEVICES 

This  section  describes  results  of  studies  of  radiation  effects  on  charge- 
coupled  devices.  Emphasis  is  placed  on  determining  the  basic  mechanisms  under¬ 
lying  the  various  observed  phenomena.  Work  performed  during  the  first  half  of 
this  program  is  described  in  Section  4.1;  Section  4.2  presents  results  of 
additional  work  which  was  performed  during  the  second  half.  The  primary  focus 
was  on  examination  of  radiation-induced  increases  in  dark  current  density.  Test 
vehicles  employed  include  MOS  capacitors  and  CCDs  fabricated  by  both  Hughes 
Aircraft  and  Fairchild  Semiconductor.  The  MOS  capacitors  were  used  to  obtain 
a  neutron  damage  coefficient  appropriate  for  describing  increases  in  the  gener¬ 
ation-center  density  in  silicon  depletion  regions.  The  rate  of  increase  in 
dark-current  density  for  neutron-irradiated  CCDs  was  measured  and  comparison 
then  made  with  results  of  calculations.  The  effect  of  a  moderate  reduction  in 
operating  temperature  on  the  response  of  CCDs  to  neutron  bombardment  is  des¬ 
cribed.  Examination  of  the  effects  of  ionizing  radiation  on  CCDs  revealed  an 
anomalous  increase  in  dark  current  density  for  Hughes  Aircraft  devices  at 
relatively  low  total  doses.  At  high  doses,  observed  increases  in  dark  current 
are  attributed  to  introduction  by  Co60  gamma  rays  of  generation  centers  in  the 
depletion-region  bulk  and  at  the  SiO^-Si  interface.  Results  of  short-term 
annealing  studies  on  CCDs  following  pulsed  neutron  bombardment  are  alsc  pre¬ 
sented. 
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4.1.1  Introduction 

Charge  transfer  devices  (CTDs)  are  finding  increased  application  for 
various  electronic  functions,  including  analog  signal  processing,  digital  data 
storage,  and  image  sensing.  Included  in  this  class  of  integrated  circuits  are 
charge-coupled  devices  (CCDs)  and  charge- injection  devices  (CIDs).  Charge 
transfer  devices  are  fabricated  using  silicon  MOS  technology,  and  thus  are 
susceptible  to  the  transient  and  permanent  effects  of  ionizing  radiation.  Cer¬ 
tain  electrical  properties  of  CTDs  will  also  degrade  due  to  displacement  damage 
produced  by  neutron  bombardment.  In  this  section,  results  of  a  study  of  the 
basic  mechanisms  of  neutron  effects  on  CTDs  are  presented. 

Several  investigations  of  the  influence  of  neutron  bombardment  on  the 
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operation  of  CCDs  have  been  made  in  the  past  few  years.  The  primary  effects 
are:  (a)  an  increase  in  dark  current  due  to  the  introduction  of  generation 
centers;  (b)  a  decrease  in  charge  transfer  efficiency  due  to  the  introduction 
of  trapping  centers.  Emphasis  in  the  present  investigation  was  placed  on  ex¬ 
amining  the  mechanisms  of  dark  current  increases.  As  an  oid  to  obtaining 
mechanistic  information,  neutron  effects  on  MOS  capacitors,  which  are  the  basic 
cells  of  a  CYD,  were  studied  extensively.  The  MOS  capacitors  were  used  as  test 
vehicles  in  the  determination  of  damage  coefficients  that  are  applicable  to  the 
radiation  response  of  CCDs.  Also  discussed  in  this  section  is  the  degree  of 
radiation  tolerance  achievable  in  a  neutron-irradiated  CCD  by  a  moderate  re¬ 
duction  in  device  operating  temperature. 

A. 1.2  Experimental  Considerations 

4. 1.2.1  Devices 

MOS  capacitors  used  in  this  study  were  fabricated  by  Hughes  Aircraft 
(HAC)  using  radiation  hardening  procedures.  Devices  were  prepared  on  both 
n-type  (~3fi-cm)  and  p-type  (<-0.8  n -cm)  <100 > silicon  substrates  with  both  wet 
and  dry  oxides.  Wet  (pyrogenic)  oxidation  was  performed  at  925°C,  followed  by 
a  N9  anneal  at  925°C.  Cry-oxide  growth  was  performed  at  1000°C.  Circular  A1 
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electrodes  (~25-mi1  diameter)  were  evaporated  from  a  graphite  crucible;  an 
anneal  was  subsequently  performed  at  500°C  in  Ng.  The  oxide  thickness  in  the 
wet-oxide  case  was  ~1000i  for  the  three  lots  studied,  which  are  referred  to  as' 
KA-09  (n-type),  WET  N  (also  n-type),  and  WET  P.  The  thickness  of  the  dry 
oxides  was  ~900A  for  the  two  lots  studied  (DRY  N  and  DRY  P).  The  CCDs  Investi¬ 
gated  were  n-bcried-channel  devices  also  fabricated  by  Hughes  Aircraft.  These 
units  are  128-bit  analog  shift  registers  employing  a  four-phase  clock.  They 
are  hardened  to  the  effects  of  total  ionizing  dose;  processing  details  have  been 
described  previously."^  (Devices  studied  were  from  two  wafers  (#15  and  #19) 
of  lot  BN-09  which  were  fabricated  using  Mask  #2128.) 


4. 1.2.2  Measurement  Techniques  for  MOS  Capacitors 

Dark  current  in  a  charge  transfer  device  is  governed  by  thermal  genera¬ 
tion  of  carriers  at  centers  in  the  depletion  region  and  at  the  oxide-semicon¬ 
ductor  interface.  As  noted  above,  the  fundamental  cell  of  a  CTD  is  an  MOS 
capacitor.  Schroder  and  Nathanson  pointed  out  that  for  an  MGS  capacitor  biased 
into  depletion  there  exists  a  laterally  depleted  region  in  addition  to  that 
directly  under  the  gate  electrode.  Figure  4-1  depicts  this  situation  for  an 
MOS  capacitor  on  n-type  silicon  with  an  applied  bias  Vg  which  is  more  negative 

than  the  flatband  voltage  V*.  .  Effects  associated  with  the  laterally  depleted 

fb  ™  an  di 

region  have  been  discussed  by  Pierret  and  Small,  ’  and  by  Takino.  In  the 

study  of  silicon  properties  using  an  MOS  capacitor,  there  are  three  ''times" 
of  interest:  (1)  depletion-region  storage  (or  charge-up)  time  t$;  (2)  minority  - 
carrier  generation  lifetime  Tg;  (3)  minority-carrier  recombination  lifetime  Tf. 
Storage  time  is  defined  as  the  time  required  for  an  inversion  layer  to  be  es- 
tabished  ir.  an  MOS  capacitor  driven  into  deep  depletion.  This  quantity  is 
readily  obtained  from  measurements  of  capacitance-versus-time  (C-t)  following 
application  of  a  step  voltage  which  drives  a  device  into  a  deeply-depleted 
state.  Establishment  of  an  inversion  layer  in  such  a  transient  measurement 
occurs  via  thermal  generation  of  carriers  in  the  space-charge  region  and  at  the 
interface.  Generation  lifetime  is  a  characteristic  time  associated  with  the 
thermal  generation  of  carriers  in  the  depletion  region  bulk.  Recombination 
lifetime  in  the  neutral  bulk  has  its  usual  definition.  In  an  unirradiated  MOS 
capacitor  at  room  temperature,  and  are  typically  on  the  order  of  micro¬ 
seconds  whereas  t  is  on  the  order  of  seconds. 


Figure  4-1.  Schematic  of  an  MOS  capacitor  cn  n-type  silicon 
with  an  applied  bias  Vg  more  negative  than  the  flatband 
voltage.  The  laterally  depleted  region  is  illustrated. 

Generation  lifetime  is  of  fundamental  importance  because  it  provides  a 

direct  measure  of  properties  of  bulk  generation  centers,  and  thereby  provides 

information  regarding  radiation-induced  dark  current  increases  in  a  CTO.  To 

obtain  t  ,  a  typical  approach  has  been  to  make  a  C-t  measurement  and  interpret 
"  no 

the  results  using  Zerbst  analysis.  ”  In  principle,  this  analysis  yields  the 
generation  lifetime  and  the  surface  generation  velocity  at  the  Si0„-5i  inter¬ 
face.  In  practice,  there  are  potential  problems  with  Zerbst  analysis,  as  dis- 
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cussed  by  various  workers.  ’’’’  A  major  problem  is  the  proper  separation 

of  surface  and  bulk  generation  processes.  Surface  generation  velocity  s  (and 
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surface  recombination  velocity)  is  a  function  of  surface  potential.  The  max¬ 
imum  value  of  s,  denoted  as  s0,  will  occur  when  a  surface  is  depleted.  The 
value  of  s  in  either  accumulation  or  inversion  will  be  significantly  less  than 
sQ.  During  a  C-t  measurement  in  which  a  device  is  pulsed  from  accumulation  in¬ 
to  deep  depletion,  the  surface  generation  velocity  under  the  gate  electrode 
after  pulsing  should  initially  equal  sq  but  will  then  decrease  as  a  function  of 
time  as  the  inversion  layer  forms.  This  layer  does  not  form  in  the  laterally 
depleted  region,  however.  Thus,  the  situation  shown  in  Figure  4-1  exists  in 
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Which  the  lateral  region  is  character!  zed  by  sQ  and  the  pate,  retd  on  is  charac¬ 
terized  by  s. 

In  a  C-t  measurement,  carrier  generation  at  interface  states  under  the 

gate  electrode  can  be  made  neglibible  compared  to  other  generation  components 

by  pulsing  a  device  from  inversion  to  stronger  inversion  instead  of  pulsing  from 

40 

ajcumulation  to  inversion.  In  this  situation,  one  obtains  an  effective  gener¬ 
ation  lifetime  t  from  analysis  of  C-t  data  since  generation  at  the  laterally 
depleted  surface  car.  still  be  important.  The  effect  of  the  lateral  component 
diminishes  as  the  size  of  an  MOS  capacitor  increases.  In  the  present  study, 
devices  employed  were  large  enough  such  that  generation  in  the  lateral  region 
was  unimportant.  Thus,  generation  lifetime  associated  purely  with  the 
depletion-region  bulk  was  obtained. 


Measurements  of  storage  time  were  made  using  a  PAR  Model  410  1-MHz  CV 
system.  Capacitors  on  an  n-type  substrate  were  pulsed  from  -5  to  -10V;  devices 
011  a  p-type  substrate  were  pulsed  from  +5  to  +10V.  Storage  time  was  defined 
here  as  the  time  reauired  to  reach  the  [l-(l/e)]  recover1.'  point  on  a  C-t  curve. 
That  is,  if  the  inversion  capacitance  is  denoted  by  Cdf  and  the  initial  Capaci¬ 
tance  in  deep  depletion  by  Cdo,  then  is  defined  here  as  the  tir.ie  required 
for  the  capacitance  to  recover  to  a  value  equal  to  Cdc  +  0.53(Cdf.  -  Cdo).  (The 
implications  of  this  t$  definition  on  the  relationship  between  t$  and  are 
discussed  below.) 


Generation  lifetime  t  was  determined  by  Zerbst  analysis  of  capacitance- 
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versus-tin«e  data  obtained  in  the  above  manner.  Schroder  and  Nathans  on"'0  modi¬ 
fied  the  analysis  of  Zerbst  to  account  for  carrier  generation  in  the  laterally 
depleted  region  shown  in  Figure  4-1.  The  differential  equation  obtained  is 
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In  this  expression  C  is  time-dependent  capacitance,  C  is  oxide  capacitance, 

vA 

is  donor  concentration,  sQ  and  s  are  defined  in  Figure  4-1,  find  r  is  the 
radius  of  the  MOS  capacitor.  The  factor  [(I/t^)  +  (4sQ/r) 3  is  termed  the  ef¬ 
fective  generation  lifetime  t  .  The  second  term  was  determined  to  be  negli¬ 
gible  compared  to  the  first  for  the  present  experimental  conditions,  so  that 

v (=  t  )  was  obtained, 
g  ge 


Zerbst  analysis  involves  obtaining  C-t  data  and  then  determining  the 
quantities 


A  computer  was  used  here  to  determine  these  quantities  based  on  closely-spaced 
discrete  data  points  from  C-t  plots.  Figure  4-2  shows  an  example  pre-irradiation 
Zerbst  plot  obtained  in  the  present  study.  The  slope  of  the  straight-line  fit 
shown  yields  Tg  and  the  y- intercept  yields  s.  Note  that  although  generation  at 
the  Si -SiO^  interface  is  suppressed  by  driving  the  device  from  inversion  to 
stronger  inversion  in  a  C-t  measurement,  finite  values  for  s  are  still  obtained. 
However,  they  are  quite  small,  as  anticipated  (e.g.,  the  results  of  Figure  4-2). 


Several  comments  are  in  order.  Zerbst  analysis  as  presented  above  as¬ 
sumes  a  uniform  doping  concentration  in  the  semiconductor.  Takino^  has  modi¬ 
fied  the  Zerbst  analysis  to  take  into  account  the  effects  of  a  nonuniform 


Figure  4-2.  Zerbst  plot  for  an  unirradiated  MOS  capacitor.  Values 
for  generation  lifetime  and  surface  generation  velocity  are 
shown. 
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doping  profile.  We  determined  doping  profiles  for  the  present  samples  by  ap¬ 
propriate  computer-aided  analysis  of  C-V  curves  and  found  them  to  be  reasonably 
uniform.  Even  if  they  were  somewhat  nonuniform,  Takino's  modification  would 
only  have  a  second-order  effect  on  the  present  experimental  results. 

44 

Simmons  and  Wei  noted  that  the  carrier  generation  rate  during  C-t  meas¬ 
urements  is  not  necessarily  constant  and  that  the  electron  and  hole  generation 
rates  can  differ.  In  particular,  if  the  bulk  generation  center  density  (N^)  is 
high  and  if  the  voltage  step  (aV )  applied  is  small,  non- steady-state  generation 
will  result  and  Zerbst  analysis  will  not  necessarily  yield  an  accurate  genera¬ 
tion  lifetime.  For  the  present  experiments,  where  aV  =  5V,  using  the  analysis 
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jf  Simmons  and  Wei  we  have  determined  that  as  long  as  Ng  ^  10  cm  ,  the 

effects  they  discuss  are  unimportant.  It  is  quite  reasonable  to  assume  that 

this  condition  holds  before  irradiation  and  after  relatively  low  fluences  of 
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neutron  bombardment  since  10  cm  is  a  rather  high  generation  center  density 
for  silicon.  As  shown  below,  data  obtained  at  all  fluences  are  mutually  con¬ 
sistent  and  yield  a  unique  value  of  generation-lifetime  damage  coefficient.  The 
indication  is  that  we  can  safely  ignore  the  effects  discussed  by  Simmons  and 
Wei  in  the  present  work. 

4. 1.2. 3  Measurement  Techniques  for  CCDs 

Two  methods  were  used  in  the  measurement  of  dark  current:  (1)  integra¬ 
tion  method;  (2)  continuous  clocking  method.  In  both  approaches,  no  charge  is 
input  to  the  CCD;  the  input  diode  is  reverse  biased.  In  the  integration  tech¬ 
nique,  the  shift  register  clocks  are  stopped  for  a  fixed  integration  time 
(typically  10  to  300  msec),  and  then  the  thermally-generated  dark  charge  is 
clocked  to  the  output  transistor  at  a  relatively  fast  rate  (typically  500  kHz). 
In  a  four-phase  clock  system,  two  adjacent  phases  are  held  high  while  the  other 
two  are  held  low.  The  reset  line  is  active  (i.e.,  high)  during  the  integration 
period.  The  total  integration  time  is  the  sum  of  the  time  during  which  the 
clocks  are  static  and  the  time  it  takes  to  clock  the  charge  to  the  output.  Use 
of  a  high  clock  frequency  makes  the  latter  contribution  negligible. 

The  continuous  clocking  method  involves  measuring  the  output  signal  with 
the  clocks  freely  running  at  some  frequency  f  .  To  compare  this  approach  with 
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Figure  4-3.  Dark  output  voltage  versus  integration  time  for  an  un¬ 
irradiated  CCD.  Measurements  obtained  by  two  methods  are  il¬ 
lustrated  (t^  =  128/f c  for  continuous  clocking  method). 

the  integration  method,  one  can  convert  to  an  equivalent  integration  time  by 
the  relation  128/f(.  for  the  HAC  CCD.  That  is,  the  integration  time  in  this 
method  is  the  time  it  takes  to  clock  a  single  bit  through  the  CCD.  Typical  data 
obtained  by  both  dark  current  measurement  techniques  are  shown  in  Fig¬ 
ure  4-3  for  an  unirradiated  CCD.  Plotted  is  the  measured  output  voltage 
(  =  dark  voltage  V^)  versus  integration  time,  where  this  time  is  defined  in  the 
manner  described  above  for  each  method.  The  expected  linear  variation  with  in¬ 
tegration  time  is  exhibited,  and  the  agreement  between  the  two  methods  is  good. 
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The  quantity  of  practical  interest  in  a  CTD  is  the  dark  current  density 
Jd>  This  quantity  is  given  by  Id/A,  where  Id  is  dark  current  and  A  is  the  area 
of  one  bit  (For  a  four-phase  CCD,  this  area  is  that  under  four  adjacent  gate 
electrodes.)  Dark  current  density  can  be  expressed  as 

Jd  -  KVd/Atl  ,  (?) 


where  t^  is  integration  time  and  K  is  a  constant  for  a  given  device.  To  evalu¬ 
ate  K,  one  needs  to  measure  current  through  the  CCD  and  the  corresponding  output 
voltage  for  some  reference  time  trg^.  For  twelve  HAC  CCDs,  Y.  varied  from  7  to 
9  x  10“13  A-sec/V.  The  effect  of  neutron  bombardment  on  this  parameter  was 

slight  for  the  fluences  employed.  The  data  of  Figure  4-3  yield  a  value  for 

7  ° 

dark  current  density  of  15  nA/cm  (A  =  38.7  x  10  cm*'). 


All  devices  were  irradiated  at  room  temperature,  which  ranged  from  ~24 
to  ~33°C.  For  purposes  of  comparison,  all  CCD  data  were  corrected  to  ?5°C  using 
a  measured  dark-current  temperature  dependence  of  exp(-0.63  e//kT). 


4. 1.2.4  Irradiation  Procedures 

All  irradiations  were  performed  using  the  TRIGA  Mark  F  reactor  at 
Northrop  Research  and  Technology  Center.  For  MOS  capacitors,  the  reactor  was 
operated  in  a  steady-state  mode  at  power  levels  of  both  5  and  50  kW.  Devices 
being  irradiated  were  unbiased,  with  their  leads  imbedded  in  conductive  foam. 

A  2.3-inch  lead  shield  plus  a  0.25-inch  boral  shield  were  in  place  in  the  dry 
exposure  room.  An  additional  2  inches  of  lead  shielding  was  placed  in  front  of 
each  device  for  the  purpose  of  enhancing  the  neutron- to-gamma  ratio.  Neutron 
fluences  were  determined  using  sulfur  dosimeters;  a  pi utonium-to- sulfur  ratio 
of  11.9  was  used  to  obtain  >  10  keV  fluences.  Gamma  dosimetry  was  performed 
using  CaF2  thermoluminescent  dosimeters.  For  CCDs,  the  reactor  was  operated  in 
both  pulsed  and  steady-state  modes.  Shielding  was  the  same  as  given  above,  as 
were  the  steady-state  reactor  operating  conditions.  For  pulsed  operation,  a 
pulse  with  $1.90  reactivity  was  used  and  the  pulsewidth  was  ~20  msec. 

In  order  to  minimize  long-term  annealing  effects  between  the  time  of 
irradiation  and  the  time  of  post-irradiation  measurements,  irradiated  MOS 
capacitors  were  placed  in  storage  at  dry-ice  temperature  within  one-half  hour 
following  irradiation,  A  device  to  be  characterized  post-irradiation  was 
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removed  from  dry  ice  and  left  at  room  temperature  for  one  and  one-half  hours 
in  order  to  achieve  stability.  Measurements  of  C-t  were  then  made. 


4.1.3  MOS  Capacitor  Results 

Measurements  of  capacitance  versus  time  were  made  on  MOS  capacitors  be¬ 
fore  and  after  neutron  bombardment,  and  these  data  were  used  to  obtain  storage 
time,  generation  lifetime,  and  surface  generation  velocity.  The  following  ex¬ 
pressions  were  employed  to  define  specific  damage  coefficients: 


1.1.  ft 

T9  Tgo  \ 
s (post)  =  s(pre)  +  K*  Y  • 


(4) 

(») 


In  Eq.  (3),  rso  and  ly  are  values  of  storage  time  measured  pre-  and  post¬ 
irradiation,  respectively,  <f>  ^s  neutron  fluence,  and  K$  is  storage-time  neutron 
damage  coefficient.  Similar  definitions  apply  to  Eq.  (4).  Equation  (5)  relates 
pre-  and  post-irradiation  values  of  surface  generation  velocity;  Y  is  the  total 
ionizing  dose  received  during  neutron  bombardment  (primarily  due  to  gamma  rays) 
and  is  the  damage  coefficient  associated  with  the  increase  in  r.  due  to  the 
effects  cf  ionizing  radiation.  (Typical  pre- irradiation  values  of  storage  time, 
generation  lifetime,  and  surface  generation  velocity  are  shown  in  Figures  4-2 
ano  4-7. ) 


Results  of  storage-time  measurements  for  devices  on  both  n-  and  p-type  sub¬ 
strates  are  shown  in  Figure  4-4.  In  each  case,  the  data  are  fit  quite  well  by  a 
unity-slope  line  as  shown,  and  the  damage  coefficients  obtained  are  as  follows: 
Ks(n-type)  =  1.2  x  10^  n-sec/cm2;  K$(p-type)  =  5.0  x  10^  n  sec/cm2.  The  cor¬ 
responding  generation  lifetime  data  obtained  using  Zerbst  analysis  are  presented 
in  Figure  4-5.  Linear  degradation  is  again  observed,  data  for  devices  on  n-  and 
p-type  substrates  are  in  agreement,  and  the  following  damage  coefficient  is  ob- 
tained:  KQ  =  7.0  x  10  n-sec/cm  .  The  dashed  lines  shown  in  Figure  4-5  correspond 
to  expected  behavior  for  recombination  lifetime,  ly,  in  neutral  bulk  material  (and 
the  associated  damage  coefficient  Kf)  based  on  the  following  low-injection-level 
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Figure  4-4.  Storage  time  degradation  in  neutron- irradiated  MOS 
capacitors.  The  unity-slope  fits  yield  storage-time  damage 
coefficients. 

values:1'  Kr(n-type)  ^  1  x  103  n-sec/cm  ;  Kr(p-type)  ^  2  x  10  n-sec/cm  .  Thus, 
values  for  Kg/Kr  are  on  the  order  of  70  and  35  for  n-  and  p-type  material,  re¬ 
spectively.  Note  in  Figure  4-5  that  there  is  no  dependence  o*  the  generation 
lifetime  data  on  oxide  growth  condition.  No  such  dependence  is  expected  for 
damage  introduction  in  the  depletion  region  bulk.  Also,  the  fact  that  results 
for  n-  and  p-type  material  are  in  agreement  indicates  that  the  dominant  neutron 
induced  generation  center  is  the  same  in  both  cases. 


Figure  4-5.  Generation  lifetime  degradation  in  neutron-irradiated 
MGS  capacitors.  The  unity-slope  fit  yields  generation¬ 
lifetime  damage  coefficient. 

Surface-generation-velocity  results  are  shown  in  Figure  4-6.  Some  scatter 
is  evident  but  the  trends  are  clear,  and  the  following  damage  coefficients  re¬ 
sult:  K^(n-type)  =  6.5  x  10~3  cm/s-rad;  K^(p-type)  =  1.5  x  IQ”3  cm/s-rad.  It  is 
interesting  to  note  that  significant  changes  in  surface  generation  velocity  are 
produced  by  a  total  dose  of  less  than  100  rads.  The  observed  increase  in  s  can 
most  likely  be  attributed  to  the  introduction  by  ionizing  radiation  of  interface 
states.  In  previous  studies  by  various  workers  on  radiation-hardened  oxides, 
such  as  those  contained  in  the  present  devices,  substantially  higher  total 
doses  were  required  before  significant  changes  in  interface-state  density  were 
observed.  This  suggests  that  the  present  experimental  approach  may  be  a  sensi¬ 
tive  tool  in  basic  studies  of  interface  states  introduced  by  ionizing  radiation.* 


*It  is  assumed  here  that  the  observed  change  in  s  (Figure  4-6)  Is  an  effect  of 
ionizing  radiation.  However,  experiments  and  calculations  performed  after  the 
completion  of  this  program  revealed  that  this  is  actually  a  neutron  effect. 
Information  on  this  subject  is  contained  In  a  paper  by  J.R.  Srour  et  al .  that 
will  appear  in  the  December  1979  issue  of  the  IEEE  Transactions  on  Nuclear  Science. 
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Figure  4-6.  Surface  generation  velocity  increase  in  irradiated  MOS 
capacitors.  The  unity-slope  fits  yield  surface-generation- 
velocity  damage  coefficients. 


Table  4-1  lists  the  damage  coefficient  values  obtained  here  and  also 

gives  previously  obtained  values^®  of  K  .  There  is  more  than  six  orders  of 

r 

magnitude  difference  between  values  for  K$  and  and  nearly  two  orders  of  mag¬ 
nitude  difference  between  values  for  and  (for  n-type).  The  versus  Kf 
difference  is  considered  first. 


Table  4-1.  Experimentally  determined  carnage  coefficients.  Values 

for  K  are  from  Reference  46. 
r 


n-type 

p-type 

Ks 

1.2  x  1013  n-sec/cm^ 

5.0  x  10*3  n-sec/cm^ 

_ K1 

7.0  x  106 

7.0  x  106 

K 

r 

LD 

O 

rH 

X 

«“H 

2  x  105 

KS 

Y 

~  6. 5  x  10  J  cm/s-rad 

_3 

—  1.5  x  10  cm/s-rad 

47 

Heiman  performed  a  detailed  analysis  of  a  capacitance-versus-time 
experiment  and  obtained  the  following  expression: 


( 


(CHf/C)  - 
<Cdf/Cdo>  ' 


\ . .  iidf 

Cdo '  Nd  Cox 


(6) 


In  this  equation,  G  is  the  carrier  generation  rate  in  the  depletion  region  and 
the  other  terms  were  defined  above.  (Equation  (6)  differs  from  Eq.  (4)  of 
Heiman  by  the  factor  on  the  right-hand  side,  presumably  due  to  a  typo¬ 
graphical  error  in  his  paper.)  If  the  second  term  on  the  left-hand  side  of 
Eq.  (6)  is  negligible  compared  to  the  first,  a  simple  exponential  solution  for 
capacitance  as  a  function  of  time  results,  with  the  time  constant  being  equal 
to  Cqx/Q  C^.  For  the  present  conditions,  the  first  term  dominates  only 
during  the  last  third  of  the  capacitance  recovery  process.  However,  as  a 
first-order  approximation,  we  neglect  the  second  term. 
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As  noted  above,  t$  was  determined  experimentally  from  the  [1  -  ( 1/e) J 
ocint  in  capacitance  recovery.  However,  Eq.  (6)  is  based  on  recovery  (or  re¬ 
duction)  of  the  depletion  width,  i.e.,  on  the  generation  of  charge,  This  con¬ 
sideration  results  in  a  slight  modification  when  comparing  the  present  experi¬ 
mental  results  with  analysis.  For  our  conditions,  recovery  of  the  depletion 
width  to  the  63"  point  corresponds  to  recovery  of  capacitance  to  the  45%  point. 
This  consideration  results  in  a  reduction  of  the  experimental  values  of  t&  and 
«s  by  a  factor  of  ~0.75  if  they  are  to  be  compared  with  calculations  using 
Eq.  (6). 

The  carrier  generation  rate  in  a  depletion  region  is  given  by  n-j/2Tg. 

'.-/e  therefore  obtain 


_!k 


/i?d\ 

\"i  /  9  * 


and  thus 


Ks =  £ 


For  the  n-type  devices  used  here,  r.q.  (3)  yields  <  value  for  K  /K  of  1.0  x  10° 

6  s  9 

as  comparea  to  the  experimental  value  of  1.7  x  10  .  Upon  reducing  the  latter 
value  by  a  factor  of  0.75,  we  obtain  1.3  x  10®  which  is  in  reasonable  agreement 
with  the  calculated  value.  For  MOS  capacitors  on  a  p-type  substrate,  the  ex¬ 
perimental  value  of  K  /K  is  ~7  x  10®,  which  reduces  to  ~5  x  10®  when  the  0.75 

s  9 

factor  is  appli^  :.  The  calculated  value  for  K  /K  in  the  p-type  case  is  ~7  x 
6  s  9 

10  ,  which  agrees  reasonably  wel 1  with  the  (reduced)  experimental  resulc. 


We  now  discuss  the  reason  for  the  large  differences  between  generation¬ 
lifetime  and  recombination-lifetime  damage  coefficients.  As  noted  above,  values 
for  Kg/«r  are  ~7G  (n-tyoe)  and  ~35  (p-type).  Qualitatively,  this  discrepancy 
is  expected  and  two  contributing  factors  can  be  citrd  to  account  for  it.  The 
first  factor  follows  directly  from  consideration  of  a  one-level  Hall-Shockley- 
Read  (HSR)  model.  It  can  be  shown  that  for  equal  concentrations  of  recombination 
centers  (in  t1,0  neuttal  bulk)  end  generation  centers  (in  the  depletion  region), 

the  inequality  <  ^  must  hold.  If  generation  in  the  depletion  region  occurs 

9  • 
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through  the  same  level  as  recombination  in  the  neutral  bulk,  and  if  this  level 
is  at  midgap,  then  =  t^.  In  general,  one  expects  the  dominant  generation 
level  to  be  at  or  near  midgap  but  recombination  is  likely  to  be  dominated  by 
some  other  level,  so  that  >  t^,  and  thereby  Kg  >  K.„.  A  second  reason  for 
expecting  Kg  to  be  greater  than  Kr  involves  the  nature  of  neutron  damage.  For 
bulk  material.,  neutron-produced  disordered  regions  act  as  potential  wells  for 
minority  carriers  and  enhance  their  recombination  rate.^®»^9  The  interior  of 
such  a  region  is  near  intrinsic,  and  thus  its  associated  potential  well  will  be 
nearly  negligible  in  a  depletion  region.  Defect  levels  introduced  by  neutrons 
will  increase  the  carrier  generation  rate  in  a  depletion  region  and  thus  de¬ 
crease  t  However,  the  enhancing  effect  of  disordered-region  potential  wells 
on  bulk  recombination  appears  to  have  no  direct  counterpart  in  a  space-charge 
region.  That  is,  no  corresponding  enhancement  of  carrier  generation  will  occur 
beyond  the  simple  addition  of  new  generation  centers.  Thus,  once  again  we  ex¬ 
pect  Kg  to  be  greater  than  Kr.  At  low  injection  levels,  a  defect  within  a  dis¬ 
ordered  region  is  on  the  order  of  a  factor  of  5  to  10  times  more  effective  for 
recombination  than  if  it  were  isolated.^®  Clustered  defects  in  a  depletion 
region  should  tend,  in  terms  of  their  effect  on  electrical  properties,  to 
resemble  isolated  defects.  This  recombination  information  thus  suggests  that 
the  difference  in  values  of  Kg  and  Kr  will  be  up  to  an  order  of  magnitude  larger 
than  that  expected  on  the  basis  of  the  above  HSR  argument.  Taking  these  two  ef¬ 
fects  together,  the  large  difference  noted  here  between  Kg  and  K^  appear 
reasonable. 

Measurements  of  storage  time  in  MOS  capacitors  were  made  as  a  function  of 
temperature  for  the  purpose  of  gaining  additional  insight.  Figure  4-7  shows 
storage  time  versus  reciprocal  temperature  for  two  devices  from  lot  KA-09.  In 
the  pre-irradiation  case,  band-to-band  excitation  is  the  dominant  mechanism  at 
high  temperatures  (1000/T<  3.0)  as  evidenced  by  an  activation  energy  approxi¬ 
mately  equal  to  the  bandgap  energy.  At  lower  temperatures,  the  temperature  de¬ 
pendence  is  dominated  by  that  of  .  The  activation  energy  at  such  temperatures 
is  the  same  before  and  after  bombardment  which  suggests  that  the  energy  level 
position  of  the  dominant  generation  center  is  nearly  the  same  in  both  cases. 

To  explore  this  point  further,  Zerbst  analysis  was  performed  on  the  C-t 
plots  used  to  obtain  the  Tg  data  in  Figure  4-7.  Results  of  this  analysis  are 
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Figure  4-7.  Storage  time  versus  reciprocal  temperature  for  irradiated 
and  unirradiated  MOS  capacitors  on  an  n-type  substrate. 


shown  in  Figure  4-8  where  generation  lifetime  is  observed  to  be  only  very 

weakly  dependent  on  temperature  both  before  and  after  irradiation.  Such  be- 

.  .  51 

havior  is  reasonable  based  on  a  one-level  expression  for  t  : 


T 

g 


cDPi  *  c,ini 

2cP  cn  Ng  ni 


(9) 


(All  terms  have  their  usual  definition.)  If  the  dominant  generation  level  is  at 
midgap,  then  n^  2:  p^  and  the  temperature  dependence  of  n^  (or  p^)  will  be  equal  to 
that  of  n...  Therefore,  this  strong  temperature  dependence  cancels,  leaving  only 
that  of  the  capture  probabilities.  Capture  probability  temperature  dependences 
typically  are  rather  weak.  The  weak  temperature  dependences  shown  in  Figure  4-8 
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Figure  4-8.  Generation  lifetime  and  surface  generation  velocity  | 

versus  reciprocal  temperature  for  irradiated  and  unirradiated  ) 

MOS  capacitors  on  an  n-type  substrate.  \ 


are  thus  attributable  either  to  capture  probabilities  or  to  a  slight  deviation  J 
of  the  dominant  generation  level  from  midgap.  Since  there  is  no  strong  differ-  * 
ence  between  the  temperature  dependence  pre-  and  post-irradiation,  the  indica-  ’ 
tion  is  that  the  dominant  radiation-induced  generation  level  is  coincident  with  j 
(or  very  near)  the  corresponding  pre-irradiation  level.  The  results  of  Fig-  j 
ure  4-8  also  reveal  that  is  only  very  weakly  temperature  dependent  over  the  ; 
range  studied  and  that  surface  generation  velocity  is  temperature  dependent  | 
over  portions  of  the  range  examined.  The  1-tter  behavior  can  most  likely  be  j 
accounted  for  in  terms  of  the  temperature  dependence  of  parameters  associated  \ 


with  a  generation  level  (or  levels)  at  the  SiC^-Si  interface. 

4.1.4  CCD  Results 

Figure  4-9  shows  change  in  dark  output  voltage  for  several  CCDs  follow¬ 
ing  both  pulsed  and  steady-state  neutron  bombardment.  These  devices  were  un¬ 
powered  during  irradiation  and  dark  voltage  measurements  were  made  using  the 
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Figure  4.9.  Change  in  dark  output  voltage  due  to  both  pulsed  and 
steady-state  neutron  bombardment  for  CCDs  irradiated  in  an 
unpowered  mode. 


continuous  clocking  method  (f  *  10  kHz)  within  one  hour  following  bombardment. 
An  approximately  linear  dependence  of  change  in  dark  voltage  on  fluence  is  ob¬ 
served  and  a  value  for  dark-current  introduction  rate  -  JCo)/4^  °f  » 
10"11  nanoamps  per  neutron  results.  Additional  data,  obtained  by  the  integra¬ 
tion  method,  yielded  an  introduction  rate  in  agreement  with  this  value. 


An  anomalous  sensitivity  for  the  present  CCD;  to  radiation-induced  in¬ 
creases  in  dark  current  was  observed  when  devices  were  irradiated  in  a  powered 
mode  (continuous  clocking,  no  input  charge).  Dark  current  increases  were  ob¬ 
served  to  be  proportional  (apparently)  to  neutron  fluence  to  the  3.5  power. 
Further  study  revealed  that  this  sensitivity  was  due  to  ionizing  radiation  and 
was  not  a  neutron  effect.  The  ionizing  dose  accompanying  neutron  bombardment  in 
this  study  ranged  from  10  to  500  rads  (Si).  Irradiation  of  an  unpowered  CCD  in 
the  gamma-ray  environment  of  a  Co®®  source  to  a  total  dose  of  100  rads  produced 
no  observable  change  in  dark  current.  However,  irradiation  of  a  powered  device 
in  the  same  source  to  about  the  same  total  dose  produced  a  dark  current  increase 
comparable  to  that  produced  during  neutron  bombardment.  (These  data  are  pre¬ 
sented  in  Section  4.2.3.)  In  general,  the  data  indicate  a  strong  bias  dependence 
of  ionization-induced  increases  in  measured  dark  current.  It  is  not  presently 
clear  which  region  in  the  CCDs  investigated  exhibits  an  extreme  total-ionizing- 
dose  sensitivity.  (One  possibility  is  charge  buildup  in  the  radiation-sensitive 
field  oxide  which  would  then  deplete  the  silicon  and  give  rise  to  an  additional 
source  of  dark  current.  Further  discussion  of  this  possibility  is  given  below.) 
However,  this  effect  is  probably  related  to  the  present  device  design  and/or 
processing  procedures  and  therefore  is  most  likely  not  of  general  interest. 

The  results  of  Figure  4-9  are  of  general  applicability  in  that  the  basic  effect 
of  neutrons  on  dark  current  in  a  CTD  is  illustrated.  (See  additional  data  in 
Section  4.2.4. ) 


Short-term  annealing  measurements  of  dark  current  in  CCDs  bombarded  by 
a  pulse  of  neutrons  were  attempted  in  this  investigation.  To  perform  such  meas¬ 
urements,  it  is  necessary  to  operate  a  given  device  with  power  applied  during 
irradiation  since  active  measurements  are  made  as  a  function  of  time  following 
the  pulse.  Because  devices  were  powered,  the  observed  response  was  dominated 
by  the  ionizing  radiation  effect  mentioned  above.  Thus,  information  regarding 
the  transient  annealing  behavior  of  generation  centers  in  the  depletion  region 
of  HAC  CCDs  was  not  obtained.  (Such  data  for  Fairchild  devices  are  presented 
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In  Section  4.2.4.)  Since  transient  annealing  of  neutron  damage  In  silicon  Is 

52 

strongly  dependent  on  the  electron  density  present,  one  might  expect  the 
rate  of  dark-current  recovery  to  be  relatively  slow  since  the  carrier  concen¬ 
tration  in  a  depleted  region  Is  quite  low. 

Before  having  determined  that  the  observed  anomalous  response  of  the  pre¬ 
sent  HAC  CCDs,  when  irradiated  under  power,  was  an  ionizing  radiation  effect,  we 
performed  an  experiment  to  rule  out  the  possibility  that  introduction  of  neutron 
damage  in  a  depletion  region  is  bias  dependent.  The  following  measurements  were 
made  on  MOS  capacitors  on  an  n-type  substrate.  Devices  were  irradiated  with  a 
neutron  burst  (5.4  x  10  n/cm  )  under  three  bias  conditions:  (1)  DC  (-8.5V); 

(2)  clocked  to  to  -8.5 V,  5  kHz  square  wave);  (3)  grounded  gate.  These  conditions 
were  repeated  in  a  second  experiment  with  unirradiated  devices  using  a  steady- 
state  neutron  flux  to  achieve  the  same  fluence.  Pre-  and  post-irradiation  deter¬ 
minations  of  t$  and  Tg  were  performed  and  damage  coefficients  then  determined. 
Values  for  K$  were  in  good  agreement  with  each  other  (and  with  the  value  given 
in  Figure  4-4)  for  all  conditions.  The  same  statement  also  applies  reasonably 
well  for  Kg  values.  We  conclude  that  introduction  of  generation  centers  in  a 
depletion  region  is  independent  of  biasing  condition,  and  in  fact  no  such  de¬ 
pendence  was  anticipated. 


4.1.5  Discussion 

It  is  of  interest  to  calculate  the  expected  increase  in  dark  current 

for  a  neutron- irradiated  charge  transfer  device.  Dark  current  density  0d  in  a 

CCD  due  to  generation  of  carriers  in  a  depletion  region  of  width  x  .  is  given 
.  53  0 

by 

Jd  =  qni  xd/2Tg  .  (10) 


Using  Eq.  (4),  we  obtain 


J 


d 


q  ni  xdtf> 
2K9 


and  thus 


d<Jd  -  -W  l"iid 

d  <t>  2Kg 


(ID 


(12) 
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At  room  temperature,  n4  =  1.4  x  10^  cm  ;  from  Figure  4-5,  K  ■  7.0  x  10® 

?  1  9 

n-sec/cm  .  The  depletion  region  width  Is  estimated  to  be  on  the  order  of  2  to 

3  pm,  which  Includes  all  of  the  0.3  m  n-type  epl  layer  plus  a  portion  of  the 

30  o-cm  p-type  substrate*  This  yields  a  calculated  range  for  0(0^  - 

of  3.2  to  4.8  x  10  11  nA/n,  In  good  agreement  with  the  experimental  value  of 

~4  x  10" 11  nA/n  for  devices  irradiated  In  an  unpowered  condition.  Thus,  we 

have  employed  M0S  capacitors  to  obtain  the  basic  quantity  (Kg)  needed  to  predict 

neutron- induced  increases  in  0d  in  a  CCD,  and  the  resulting  agreement  between 

calculation  and  experiment  is  satisfying.  We  further  note  that  should  be  of 

general  usefulness  in  predicting  neutron  damage  effects  in  depletion  regions  for 

any  silicon  device  bombarded  by  fission  neutrons.  (See  Section  5.0.) 


Table  4-2  compares  values  of  d(J(j  -  Jcj0)/d^  derived  from  data  obtained 
by  various  workers  with  the  present  result.  It  should  be  borne  in  mind  that 
values  for  this  quantity  are  expected  to  be  in  agreement  only  if  x^  is  the  same, 
based  on  Eq.  (12),  and  also  if  all  experimental  conditions  are  the  same.  For 
purposes  of  comparison,  we  assume  that  %.  is  approximately  the  same  in  all 

Q  oo  35 

cases.  The  present  data  and  that  of  Hartsell  and  Chang  were  obtained  using  a 
fission  reactor,  and  their  values  for  d( agree  within  a  factor  of 
~2  with  the  current  value.  The  value  obtained  by  S?  cs^  is  a  factor  of  5  higher 
than  ours,  which  is  consistent  with  the  more  effective  introduction  of  genera¬ 
tion  centers  by  the  higher  energy  neutrons  he  employed  than  the  fission  neutrons 
used  here.  In  summary,  the  present  value  of  d(Jd  -  is  in  reasonable 

agreement  with  values  derived  from  data  obtained  by  Hartsell  and  Chang,  _  the 

two  cases  for  which  comparison  appears  appropriate.  One  further  comparison  can 

32  13  2 

be  made.  Williams  and  Nelson  obtained  a  value  for  K$  of  3.3  x  10  n-sec/cm 

using  a  p-surf ace-channel  device  (n-type  substrate).  This  value  falls  midway 
between  those  shown  in  Figure  4-4  and  thus  is  in  reasonable  agreement  with  the 
present  work.  However,  this  agreement  may  be  fortuitous  since  (and  there¬ 
fore  K$)  depends  on  the  magnitude  of  the  voltage  step  applied  in  a  C-t  measure¬ 
ment  and  on  the  doping  concentration.  Storage  time  and  its  associated  damage 
coefficient  do  not  have  the  fundamental  physical  significance  that  and  Kg 
have. 


Lowering  the  operating  temperature  of  a  charge  transfer  device  causes  a 
dramatic  reduction  in  dark  current  due  to  the  strong  temperature  dependence  of 
carrier  generation  rate  in  the  depletion  region.  Neutron  bombardment  introduces 

*In  subsequent  work,  described  in  Section  5.0  of  this  report,  we  concluded  that 
the  depletion  region  width  is  most  likely  larger  than  that  used  here. 
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Table  4-2.  Comparison  of  values  for  d(Jrf  -  0do  )/<*«£  derived  from 
data  obtained  by  various  workers.  ao 


Reference 

1 

d(Jd  -  Jd0>/d* 

(nA/neutron) 

Radiation 

Source 

Comments 

Hartsell 33 

~9.6  x  10-11 
~5.6  x  10"11 

Fast-burst  reactor 

Fast-burst  reactor 

Buried  channel 

Surface  channel 

c  u  34 

Saks 

1.9  x  id"10 

Cyclotron  (0  to 

30  MeV  neutrons, 
15-MeV  peak) 

Buried  channel 

Chang35 

~2  x  10"11 

Fast-burst  reactor 

Buried  channel 
and 

Surface  channel 

Present  work 

~4  x  10'11 

TRI6A  reactor 

Buried  channel 

_ 

generation  centers  in  that  region  which  give  rise  to  an  increase  in  dark  cur¬ 
rent.  The  impact  of  this  radiation- induced  increase  can  be  diminished  by  lower¬ 
ing  the  device  operating  temperature.  We  consider  such  a  hardening  approach 
here  and  present  a  graph  that  allows  the  selection  of  an  operating  temperature 
which  should  prevent  post-irradiation  dark  current  from  being  larger  than  some 
specific  value. 


The  temperature  dependence  of  Jd  [Eq.  (10)]  is  dominated  by  that  of  n.. . 

where 


n.  =  CxT3/2  exp(-Eg/2kT)  =  Cj  f(T) 


(13) 


(C^  is  a  constant  and  Eg  is  bandgap  energy.)  Upon  defining  a  second  constant, 
C^,  as  qxy2  ,  we  obtain 


C,  C„f(T)<f> 
Jd  -  JdO  +  SL - 


(14) 


Plotting  log  f(T)  versus  reciprocal  temperature  yields  a  slope  (i.e.,  an  activa¬ 
tion  energy)  of  -0.58  eV.  (For  this  determination,  we  employed  Eg  =  1.11  eV 
and  the  bandgap  temperature  dependence  was  neglected.  This  is  reasonable  over 
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the  temperature  range  considered.)  As  noted  above,  Tg  is  only  weakly  temperature 
dependent.  Using  the  pre-irradiation  value  for  the  activation  energy  of  tg  in 
Figure  4-8  (0.05 eV)  results  in  an  overall  temperature  dependence  for  JdQ  of 
exp (-0.63  eV/kTl.  This  dependence  is  in  good  agreement  with  the  ts  data  of  Fig¬ 
ure  4-7  and  in  exact  agreement  with  pre-  and  post-irradiation  data  obtained  for 
HAC  CCDs  (not  snown  here).  For  simplicity,  we  assume  that  the  second  term  in 
Eq.  (14)  also  varies  as  exp(-0.63  eV/kT)  and  thus  that  Jd  has  this  same  dependence. 
This  is  a  reasonable  assumption  since  Kg  is  only  weakly  temperature  dependent,  as 
discussed  above. 

> 

Upon  normalizing  Jd  to  JdQ,  we  obtain 


J .  t  d> 

-1=  i 


do 


(15) 


This  expression  was  used  at  303°K  (i.e.,  J.  =  1  before  irradiation  at  1000/T  = 

o  00 

3.3  K  ")  to  obtain  values  of  normalized  dark  current  density  as  a  function  of 

6  ? 

fiuence  using  K  -  7.0  x  10  n-sec/cm  .  A  value  of  r  was  selected  (14  psec) 


so  that  AJd/A«i>  would  equal  4  x  10-*1  nA/n  for  Jdo  =  20  m/cmu ,  a  tvpical  dark 
current  density  for  the  CCDs  studied.  Curves  with  a  -0.63  eV  slope  were  then 
drawn  through  the  calculated  room- temperature  values  and  results  are  shown  in 
Figure  4-10.  This  graph  illustrates  the  effectiveness  of  lowering  the  operating 
temperature  of  a  CTO  to  avoid  neutron- induced  dark  current  problems.  As  an  ex¬ 
ample,  assume  that  the  maximum  dark  current  density  that  be  tolerated  is  unity, 
the  value  of  JdQ  at  303°K  in  Figure  4-10.  If  the  expected  fiuence  is  10*3  n/cm^, 
then  this  figure  indicates  that  operating  the  device  at  T  £  269°K  will  keep  Jd  £ 
unity  following  irradiation.  Thus,  a  relatively  moderate  lowering  of  the  oper¬ 
ating  temperature  in  a  CTC  should  resuH  in  a  substantial  radiation  tolerance  in 
terms  of  neutron-induced  increases  in  dark  current. 


go 


4.1.6  Summary 

A  study  of  the  mechanisms  of  neutron- induced  dark  current  increases  in 
charge  transfer  devices  has  been  performed  using  MOS  capacitors  and  CCDs  as  test 
vehicles.  MOS  capacitors  were  used  to  obtain  neutron  damage  coefficients  as¬ 
sociated  with  storage-time  and  generation-lifetime  degradation  («s  and  Kg, 
respectively)  and  gamma-r'V  damage  coefficients  associated  with  increases  in 
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TEMPERATURE  (°K) 


Figure  4-10.  Normalized  dark  current  density  versus  reciprocal 
temperature  with  neutron  fluence  as  a  parameter.  The 
effect  of  a  reduced  operating  temperature  on  radiation 
tolerance  is  illustrated. 

surface  generation  velocity  at  the  SiOg-Si  interface.  In  terms  of  the  radiation 
response  of  a  buried-channel  CCD,  the  quantity  determined  having  the  most  funda- 
menta1  significance  is  K  .  the  generation  lifetime  damage  coefficient.  A  unique 
value  for  this  coefficient  of  7.0  x  10°  n-sec/cm  was  determined  for  devices  on 
both  n-  and  p-type  substrates.  This  value  was  then  used  to  calculate  the  ex¬ 
pected  change  in  dark  current  density  with  neutron  fluence  in  a  CTD.  The  cal¬ 
culated  value  is  in  good  agreement  with  the  present  experimental  value  of  ~4  x 
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10-11  nA/n  for  Irradiated  CCDs.  (See  additional  discussion  of  this  point  In 

Section  5.0.)  Experimental  values  for  the  ratio  of  K$/Kg  were  obtained  and 

found  to  be  grater  than  six  orders  of  magnitude.  These  substantial  differences 

were  accounted  for  reasonably  well  for  both  n-type  and  p-type  material  using 

47 

the  model  of  Helman.  A  qualitative  explanation  was  presented  to  account  for 
the  substantial  difference  between  recombination-lifetime  and  generation-life¬ 
time  damage  coefficients  which  is  based  on  consideration  of  the  Hall -Shockley- 
Read  model  and  also  on  the  nature  of  neutron  damage  In  a  depletion  region.  A 
discussion  was  presented  which  indicates  that  a  moderate  reduction  in  operating 
temperature  of  a  CTD  should  substantially  reduce  its  sensitivity  to  the  pro¬ 
duction  of  dark  current  by  neutron  bombardment. 
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4.2  ADDITIONAL  STUDIES  OF  RADIATION  EFFECTS  ON  CHARGE-COUPLED  DEVICES 

4.2.1  Introduction 

Section  4.1  of  this  report  describes  results  of  studies  of  radiation 
effects  on  CCDs  performed  during  the  first  half  of  the  performance  period  for 
the  present  contract.  In  this  section,  we  present  findings  obtained  during 
the  second  half  of  the  program.  In  addition,  data  are  presented  which  serve 
to  reinforce  and  provide  further  insight  regarding  several  of  the  points  briefly 
made  above.  The  primary  subjects  included  here  are:  e)  effects  of  ionizing 
radiation  on  Hughes  Aircraft  and  Fairchild  CCDs;  b)  transient  and  permanent 
effects  of  neutron  bombardment  on  the  Fairchild  CCD  311.  As  in  Section  4.1, 
emphasis  is  placed  here  on  investigation  of  radiation-induced  increases  in 
dark-current  density. 

4.2.2  Experimental  Procedures 

Techniques  used  to  measure  dark  current  were  the  same  as  those  described 
in  Section  4. 1.2. 3.  The  radiation  facilities  employed  were  the  TRIGA  Reactor 
and  the  Co60  source  at  Northrop  Research  and  Technology  Center.  All  measure¬ 
ments  were  made  at  room  temperature,  which  was  monitored.  Data  were  corrected 
to  298°K  based  on  the  measured  temperature  dependence  of  dark-current  in  CCDs, 
as  mentioned  above.  Two  types  of  n-buried-channel  CCDs  were  studied:  devices 
fabricated  by  Hughes  Aircraft  (described  in  4. 1.2.1  and  in  References  35  and 
54)  and  the  commercially  available  Fairchild  XCD  311.  Figure  4.11  shows  a 
block  diagram  of  and  pin  identification  for  the  CCD  311.  Note  that  this  device 
contains  two  130-bit  analog  shift  registers  which  normally  are  multiplexed. 

We  operated  the  clocks  for  these  two  shift  registers  180°  out  of  phase  so  that 
darx-current  information  for  both  was  obtained  simultaneously.  Also  note  that 
the  CCD  311  uses  a  two-phase  clock  whereas  the  Hughes  Aircraft  device  operates 
with  four-phase  clocking. 

Pertinent  schematic  diagrams  of  circuits  employed  in  conjunction  with 
the  CCD  311  are  now  presented.  Figure  4-12  shows  a  diagram  cf  the  circuit  used 
in  radiation  testing  (both  reactor  and  Co*^).  This  circuit  permits  active 
measurements  to  be  made  during  pulsed  or  steady-state  irradiation.  The  drive 
circuit  illustrated  in  block  form  in  Figure  4-12  is  shown  in  detail  in  Figure 
4-13.  Tigure  4-13a  depicts  the  clock  waveform  and  signal  input  generator.  In 
F-'gurr:  4-13b,  the  circuit  used  to  provide  power  for  the  CCD  311  and  for  the 
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PIN  NAMES 

*1A' 

02A’  ®1B'  02B 

Analog  Shift  Register  Transport  Clocks 

OG 

Output  Gate 

OS 

Output  Transistor  Source 

OD 

Output  Transistor  Drain 

CS 

Compensation  Transistor  Source 

*R 

Reset  Transistor  Gate  Clock 

RD 

Reset  Transistor  Drain 

vss 

Substrate  (Ground) 

VIA* 

V1B 

Analog  Input  Gates 

0SA’ 

0SB 

Analog  Sample  Clocks 

TP1, 

TP2,  TP3 

Test  Points 

Figure  4-11.  Block  diagram  of  and  pin  identification  for  the 
Fairchild  CCD  311. 


93 


circuit  of  Figure  4-13a  is  presented.  A  timing  diagram  showing  w&veforms 
applicable  for  the  circuit  of  Figure  4-13a  is  given  in  Figure  4-M  Details 
of  the  pin  connections  for  the  CCD  311  during  irradiation  are  shown  in  Figw* 
4-15.  Figure  4-16  shows  the  differential  amplifier  circuit  used  to  eliminate 
the  roset  pulse  at,  the  output  (refer  to  Figure  4-12).  In  Figure  4-17,  the  driVe 
amplifiers  used  in  testing  the  CCD  311  are  illustrated. 


4.2.3  Ionizing  Radiation  Effects 


As  mentioned  briefly  in  Section  4.1.4,  dark  current  density  In  Hughes 
Aircraft  CCDs  was  found  to  be  quite  sensitive  to  ionizing  radiation  when 
devices  were  irradiated  in  a  powered  mode.  We  present  here  the  data  which  led 
to  that  conclusion  along  with  additional  results  of  ionizing  irradiation 
studies. 


A  schematic  diagram  of  the  Hughes  Aircraft  buried-channel  CCD  studied  is 

shown  in  Figure  4-18.  (A  similar  design  is  used  in  the  Fairchild  part.)  In 
55 

earlier  work,  we  noted  that  the  apparent,  sensitivity  of  the  HAC  CCD  to 
neutron  bombardment  was  greatly  enhanced  if  the  device  was  irradiated  with 
power  applied.  These  previous  data  are  shown  in  Figure  4-1:'  where  the  change 
in  output  dark  voltage  (proportional  to  dark  current)  is  plotted  versus 
neutron  fluence.  The  open  symbols  represent  data  for  powered  devices  and  the 
solid  symbols  are  for  unpcwered  parts.  {’Towered’1  is  defined  here  as  the 
application  of  typical  ac  and  dc  voltages  to  all  appropriate  terminals  of  the 
CCD  with  the  exception  of  the  input.  No  charge  was  input  to  the  CCD  for  the 
purpose  of  permitting  dark  current  measurements  by  either  cf  the  methods  des¬ 
cribed  in  Section  4. 1.2. 3.  The  "unpowered"  condition  refers  to  irradiation  of 
a  CCD  with  all  of  its  leads  imbedded  in  conductive  foam.) 

An  extreme  sensitivity  of  the  powered  parts  as  compared  to  the  unp’owered 
units  is  evidenced  in  Figure  4-19.  Ionizing  radiation  is  also  incident  on  the 
CCDs  during  neutron  bombardment,  so  we  performed  an  experiment  to  determine 
whether  the  sensitivity  of  dark  current  in  parts  irradiated  under  power  was  a 
neutron  effect  or  an  effect  of  ionizing  radiation.  During  neutron  bombardment, 
CaF2  TLDs  were  used  to  measure  the  concomitant  ’cnizing  doss.  Figure  4-20 
shows  the  data  of  Figure  4-19  replotted  in  terms  of  change  in  dark  voltage 
versus  total  dose.  Also  shown  are  data  points  obtained  for  HAC  CCDs  irradiated 
in  a  Co®^  source.  Device  19-41  had  received  no  previous  irradiation.  It  was 
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Figure  4-14.  Timing  diagram  for  the  CCD  311.  Waveforms 

applicable  for  the  circuit  of  Figure  4-13a  are  shown. 


Figure  4-15.  Pin  connections  for  the  CCD  311  during  irradiation. 
(See  figure  4* 12. ) 


9? 


Figure  4-16.  Differential  amplifier  used  with  the  CCD  311  to 
eliminate  the  reset  pulse  at  the  output.  (See  Figure 


Figure  4-17.  Drive  amplifiers  used  with  the  CCD  311.  (See 
Figure  4-12.) 
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Figure  4-18.  Schematic  diagram  of  the  Hughes  Aircraft  buried 
channel  CCD. 


FIUCNCE  In/cm’,  >  lc  kcV) 


Figure  4-19.  Radiation-induced  change  in  dark  output  voltage  for 
the  Hughes  Aircraft  CCD.  Although  dark  voltage  is  plotted 
versus  neutron  fluence,  only  the  solid  data  points  are 
attributable  to  neutron  damage.  (See  text.') 
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Figure  4-20.  Change  in  dark  output  voltage  versus  ionizing  dose  for 
the  Hughes  Aircraft  CCD.  Except  for  the  two  designated  ”Co°°", 
all  data  points  shown  are  those  given  in  Figure  4-19.  In  the 
present  case,  the  ionizing  dose  received  during  neutron  bom¬ 
bardment  is  plotted  as  the  abscissa . 


first  irradiated  to  130  rads  with  power  applied  and  the  resulting  change  in 
dark  voltage  was  substantial,  as  shown  in  Figure  4-20.  The  130-rad  data  point 
shown  for  this  device  is  in  general  agreement  with  those  shown  for  other 
devices  from  wafer  #19.  This  result  reveals  that  the  extreme  sensitivity  of 
dark  current  in  powered  HAC  CCDs  is  an  effect  of  ionizing  irradiation,  and 
not  of  neutron  bombardment.  After  a  16h  anneal  at  room  temperature,  for 
Device  19-41  decreased  to  65%  of  the  13J-rad  value  shown  in  Figure  4-20.  This 
anneal  was  followed  by  an  additional  irradiation  under  power  to  a  dose  of  75 
rads.  The  change  in  dark  voltage  at  10  kHz  as  a  result  of  this  irradiation 
was  140  mV.  This  result  is  plotted  in  Figure  4-20  in  terms  of  the  incremental 
dose  and  the  incremental  change  i».  (i.e.,  75  rads  and  140  mV).  The  result 

of  this  second  irradiation  of  Device  19-41  is  consistent  with  that  obtained 
in  the  first  bombardment. 
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In  a  second  experiment,  Device  19-33  was  irradiated  in  the  Co®®  source 
with  no  power  applied.  This  device  had  been  previously  neutron  irradiated 
twice  and  had  received  ionizing  doses  of  180  rads  (unpowered)  and  170  rads 
(powered).  Three  months  later  the  Co6®  experiments  were  performed,  and  at  that 
time  the  dark  current  density  for  this  device  was  50%  larger  than  its  pre¬ 
irradiation  (virgin)  value.  The  first  unpowered  gamma  irradiation  was  to  a 
total  dose  of  100  rads,  and  no  change  in  dark  current  was  observed.  This 
result  indicates  that  the  ionizing-radiation  sensitivity  of  HAC  CCDs  to  ionizing 
radiation  arises  only  when  power  is  applied.  However,  a  second  unpowered 
irradiation  of  Device  19-33  to  an  additional  660  rads  caused  to  increase 
significantly,  as  shown  in  Figure  4-20.  (Also  shown  are  results  of  long-term 
room-temperature  anneals  of  90  and  1600  hours.  Significant  recovery  was  noted.) 
This  result  could  cause  one  to  suspect  that  the  solid  data  points  in  Figures 
4-19  and  4-20  (also  shown  in  Figure  4-9),  which  are  results  of  neutron  bom¬ 
bardment  in  an  unpowtred  state,  are  attributable  to  the  effects  of  both 
neutrons  and  gamma  rays.  However,  results  of  neutron  bombardment  of  the 
Fairchild  CCD  311  (discussed  below),  which  is  insensitive  to  ionizing  radiation 
(also  discussed  below),  are  consistent  with  findings  shown  in  Figure  4-9  for 
unpowered  HAC  CCDs.  That  is,  we  judge  the  solid  data  points  in  that  figure  to 
be  attributable  solely  to  the  effects  of  neutrons.  It  is  possible  that  results 
of  Co®®  bombardment  of  Device  19-33  were  influenced  by  the  fact  that  it  had 
previously  been  irradiated  with  power  applied.  Further  study  would  be  required 
on  virgin  devices  to  clear  up  this  issue.  (Exploration  of  this  point  was  not 
of  primary  importance  in  the  present  program.) 

For  comparison  purposes,  Fairchild  CCD  311s  were  irradiated  in  the  Co®® 
source  with  and  without  power  applied.  Results  are  shown  in  Figure  4-21  where 
dark  current  density  is  plotted  versus  total  dose.  Devices  1  and  2  were 
irradiated  unpowered.  Device  1  was  still  functioning  at  10®  rads  (the  highest 
dose  examined  for  this  unit),  and  the  increase  in  dark  current  density  was 

5 

only  ~20%.  Device  2  was  still  operational  at  3  x  10  rads  but  failed  some¬ 
where  between  that,  dose  and  10®  rads.  Devices  2?  and  25  were  irradiated  under 

power.  Device  23  exhibited  a  strong  increase  in  dark  current  (saturated)  at 

4  4 

5  x  10  rads  and  failed  to  operate  after  8  x  10  rads.  Device  25  failed  at 

9  x  10^  rads.  It  is  interesting  to  note  that  the  Fairchild  CCD  was  not  in¬ 
tentionally  fabricated  as  a  radiation-hardened  part  but  it  exhibits  moderate 
tolerance  to  ionizing  radiation. 
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Figure  4-21.  Dark  current  density  versus  ionizing  dose  for 
the  Fairchild  CCD  311. 

To  explore  further  the  sensitivity  of  HAC  CCDs  to  ionizing  radiation, 
additional  measurements  were  made.  Figure  4-22  shows  dark  current  density 
versus  total  dose  for  three  HAC  devices  irradiated  with  power  applied.  The 
results  are  rather  striking  and  deserve  considerable  comment.  Before  doing  so, 
some  experimental  differences  should  be  noted.  Devices  19-62  and  17-65  were 
previously  unirradiated  whereas  Device  19-21  had  received  an  earlier  neutron- 
plus-ganma  irradiation  (see  Table  2  in  Reference  55).  In  that  earlier  irrad¬ 
iation,  dark  current  was  observed  to  recover  to  within  a  factor  of  two  of  its 
pre- irradiation  value  after  a  thirty- day  room- temperature  anneal.  It  is  not 
presently  clear  whether  the  irradiation  history  of  Device  19-21  had  an  effect 
on  the  later  experiment  (Figure  4-2 2),  i.e.,  gave  rise  to  the  difference 
between  results  for  this  unit  and  the  other  two  devices.  However,  this  aspect 
of  the  data  is  unimportant  for  the  present  purposes.  We  concentrate  here  on 
the  trends  evident,  and  results  shown  for  all  three  devices  in  Figure  4-22  are 
qualitatively  similar.  One  additional  difference  was  that  Device  19-21  was 
operated  at  a  constant  c1ock  frequency  (500  kHz)  during  irradiation  whereas 
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Figure  4-22.  Dark  current  density  versus  ionizing  dose  for 
Hughes  Aircraft  CCOs. 

the  clock  frequency  was  varied  for  Devices  19-62  and  17-65.  This  difference  in 
procedure  is  probably  not  of  importance.  A  final  note  regarding  experimental 
procedures:  the  dark  current  measurements  of  Figure  4-22  were  made  with  the 
Co60  source  lowered  into  the  floor  (i.e.,  dose  rate  equal  tc  zero).  With  the 
source  raised,  the  presence  of  ionizing  radiation  enhances  the  dark  current  and 
also  produces  noise  in  the  measurement  system. 

Three  dark-current  regimes  are  evident  in  Figure  4-22.  A  strong  increase 

2  3 

in  •).  is  noted  in  the  range  10  -  10  rads.  This  is  followed  by  a  strong 
decrease  in  Jd  in  the  range  10^  -  10'  rads.  In  the  third  regime  (>Mr  rads), 
once  again  exhibits  an  increase.  This  third  regime  is  reasonably  well  under¬ 
stood  and  is  attributable  to  two  effects:  1)  an  increase  in  the  generation 
center  density  in  the  aepletion-region  bulk  due  to  bombardment  by  the  ~l-MeV 
Co60  ganina  rays;  2)  an  increase  in  the  generation  center  density  at  the  SiO^-Si 
interfcce  due  to  production  of  interface  states  by  the  ionizing  radiation 
(Note  that  the  first  effect  is  production  of  displacement  damage  by  ionizing 
radiation).  Both  of  these  effects  should  increase  linearly  with  dose.  Figure  4 


103 


10* 


102 

i 


10' 

10 3  JO6  107 

DOSE  [rads (Si)] 


Figure  4-23.  Change  in  dark  current  density  versus  ionizing 
dose  for  a  Hughes  Aircraft  CCD.  Data  shown  were  derived 
from  Figure  4-22. 

shows  the  change  in  dark  current  density  (i.e.,  -  JdQ)  versus  dose  for  Device 

19-21  (data  from  Figure  4-22),  and  a  linear  dependence  is  indeed  observed.  The 

^  a  p 

rate  of  change  obtained  is  1.2  x  10  nA/cm  -rad.  To  examine  how  well  that 
rate  describes  data  obtained  by  other  workers,  we  have  made  the  comparison 

35  54 

illustrated  in  Figure  4-24.  Data  derived  from  the  work  of  Chang  and  Aubuchon  ’ 
are  shown,  along  with  a  dashed  line  corresponding  to  the  rate  of  change  obtained 
in  Figure  4-23.  The-e  is  considerable  scatter  in  the  experimental  data,  which 
is  perhaps  due  in  part  to  the  fact  that  some  of  the  data  were  obtained  with 
bias  applied  and  the  rest  without  bias.  Radiation-induced  production  of  inter¬ 
face  states  should  be  enhanced  by  the  application  of  a  field  across  the  gate 
oxide.  There  is  a  slight  indication  in  Figure  4-24  that  powered  devices 
exhibited  a  larger  tncrease  in  dark  current  density,  but  further  study  is 
needed  to  confirm  this  point.  (In  the  present  program,  we  irradiated  a  HAC 
CCD  in  the  Co^°  source  with  no  power  applied  up  to  a  dose  of  3  x  10^  rads. 

Dark  current  exhibited  a  linear  increase  with  dose  at  high  doses.  However, 
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figure  4-24.  Change  in  dark  current  density  versus  ionizing  dose.  35  54 
Data  points  were  obtained  from  the  work  of  Chang  and  Aubuchon-  ’ 

The  dashed  line  is  based  on  the  fit  to  data  shown  in  Figure  4-24. 

the  rate  of  change  observed  was  a  factor  of  ~3.5  larger  than  that  shown  in 

Figure  4-23.  Clearly,  additional  investigation  of  the  bias  dependence  of 

ionizinq  radiation  effects  on  CCDs  at  high  doses  is  needed.  In  this  zero-bias 

5 

experiment,  unfortunately  no  data  points  were  recorded  for  doses  less  than  10 
rads.  In  retrospect,  such  data  would  have  been  useful  in  view  of  the  dis- 
cusion  above  in  conjunction  with  Figure  4-20.)  Note  that  at  10^  rads  in 
Figure  4-24  the  dashed  line  passes  roughly  through  the  middle  of  the  scattered 
data  points,  indicating  some  success  in  predicting  the  experimental  observations. 
At  10°  rads,  the  data  points  are  all  higher  than  the  dashed  line  and  this  may 
be  due  to  remanence  of  a  portion  of  the  anomalous  enhanced  current  evident  at. 
lower  doses  (Figure  4-22-).  We  now  discuss  this  anomalous  current. 

The  strong  increase  in  dark  current  density  at  low  doses,  shown  in  Fiqure 
4-22,  is  qualitatively  consistent  with  the  data  shown  in  Figure  4-20.  In  the 
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former  case,  the  ionizing  dose  was  obtained  from  a  Co^  source  whereas  in  the 
latter  case  the  Ionizing  radiation  was  concomit  nt  with  neutron  bombardment. 

A  starting  point  for  a  model  which  has  potential  for  describing  these  findings 
is  now  given.  The  HAC  CCDs  contain  a  non-radiati on-hardened  field  oxide,  as 
shown  in  Figure  4-18.  Not  shown  there  is  a  channel-stop  diffusion  that  is  also 
present.  Figure  4-25  is  a  schematic  cf  the  field  oxido  and  the  channel  stop. 

A  relatively  small  total  ionizing  dose  can  give  rise  to  significant  positive 
charge  buildup  if.  a  field  oxide.  This  effect  is  enhanced  if  a  bias  is 
applied  across  that  oxide  layer.  Positive  charge  buildup  will  result  in  for¬ 
mation  of  a  depletion  region  in  the  underlying  p-type  substrate.  Thermal 
generation  of  carriers  will  occur  in  this  region  and  these  carriers  could 
substantially  increase  the  CCD  dark  current  if  they  were  able  to  couple  into 
tne  channel.  However,  tne  channel  stop  should  prevent  this  from  occurring.  At 
higher  doses,  negatively  charged  interface  states  will  form  at  the  field  oxide- 
silicon  interface,  as  shown  in  Figure  4-25.  This  negative  charge  will  serve  to 
compensate  for  built-up  positive  charge  in  the  oxide,  thereby  reducing  the 
extent  of  the  depletion  region  in  the  silicon.  Thus,  fewer  carriers  would  be 
thermally  generated  in  that  region  at  higher  doses.  Once  again,  the  problem 
with  this  model  is  that  it  is  not  clear  how  carriers  thermally  generated 
external  to  the  CCD  channel  can  couple  into  the  channel  and  enhance  the  dark 
current.  The  anomalous  response  of  HAC  CCDs  to  ionizing  radiation  is  most 
likely  associated  with  a  device  processing  problem  that  has  yet  to  be  deter¬ 
mined.  Further  study  of  the  effects  of  processing  variations  on  this  response 
will  most  likely  be  required  before  the  observed  phenomena  can  be  understood 
and  counteracted. 
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Figure  4-25.  Illustration  of  charge  buildup  in  the  field  oxide 
for  a  CCD. 
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One  final  point  should  be  noted  regarding  data  shown  in  Figure  4-22.  For 

Device  19-62,  dark  current  density  was  measured  after  a  12-day  post-irradiation 

anneal  at  room  temperature.  Dark  current  increased  by  50%  over  this  period  of 

time,  as  shown  in  the  figure.  This  behavior  indicates  that  annealing  of  the  radi- 

4  5 

ation  "damage"  responsible  for  the  decrease  in  in  the  10  -  10  rads  range  pro¬ 
ceeds  at  a  faster  rate  than  annealing  of  the  "damage"  responsible  for  the  dark- 

2  3 

current  increase  in  the  10  -10  rads  range.  Qualitatively,  such  beiiavior  appears 

inconsistent  with  the  model  depicted  in  Figure  4-25  since  the  amount  of  radiation- 

induced  positive  charge  in  Si 0g  tends  to  decrease  with  time  following  irradiation 

whereas  the  interface-state  density  tends  to  remain  constant  or  increase  over  the 

56 

same  time  regime  for  a  tnermal-oxide/silicon  interface.  (To  determine  whether 
this  statement  applies  for  a  field-oxide/s  11  icon  interface  would  require  further 
study. ) 

4.2.4  Neutron  Effects 

The  relative  absence  of  ionizing  radiation  effects  on  dark  current  density 
in  the  Fairchild  CCD  311  as  compared  to  the  HAC  CCD  makes  if  a  good  test  vehicle 
for  examining  neutron  effects  cn  CCDs  tvowered  during  bombardment.,  Results  of  neu¬ 
tron  irradiations  of  the  CCD  311  are  reported  here.  In  the  first  experiment,  five 
separate  devices  were  irradiated  unpowered  to  five  different  fluences,  and  dark 
current  measurements  were  made  before  and  after  these  steady-state  irradiations. 
(Post-irradiation  measurements  were  performed  within  approximately  one  hour  after 
bombardment.)  Figure  4-26  is  a  plot  of  the  change  in  dark  current  density  versus 
neutron  fluence.  Data  over  a  fairly  wide  fluence  range  are  fit  quite  well  by  a 
unity-slope  line  as  shown  in  the  figure.  The  slope  obtained,  3.5  x  10-11  nA/neu- 
tron,  compares  quite  well  with  the  result  for  the  HAC  CCD  shown  above  in  Figure 
4-9.  (Additional  discussion  of  these  results  is  given  in  Section  5.3  of  this  re¬ 
port.  ) 

Measurements  of  dark  current  density  were  performed  as  a  function  of  time 
following  pulsed  neutron  bombardment  for  the  CCD  311.  In  the  first  experiment 
performed.  Device  J6  was  oDerated  at  a  clock  frequency  of  2  kHz  during  bombard¬ 
ment,  and  the  continuous  clocking  method  was  used  to  monitor  dark  current  as 
a  function  of  time.  This  device  was  bombarded  by  two  neutron  pulses  a  few  hours 
apart.  Both  bursts  were  ~20  msec  vride  ($1.90  pulse)  and  shielding  was  as  fol¬ 
lows:  1/4"  bora1  plus  4.3"  Pb.  In  addition,  for  the  first  pulse  7"  of  water 
shielding  was  employed  to  reduce  the  fluence.  No  water  shielding  was  used  on  the 
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FLUENCE  (n/cm2  >10  kaV) 


Figure  4-26.  Change  in  dark  current  density  versus  neutron 
fluence  for  the  Fairchild  CCD  311. 

second  burst.  Sulphur  pellets  were  used  for  neutron  dosimetry  and  CaF9  TLDs  were 

employed  to  obtain  total  ionizing  dose.  With  no  water  shielding,  the  plutonium-to- 

sulphur  ratio  is  known  to  be  11.9  for  our  shieldirg  conditions.  We  obtained  a  flu- 

ence  of  4.4  x  10  n/cm  (>10  keV)  for  the  second  pulse  (gamma  dose  ^  130  rads ( Si ) ) . 

For  7"  H90  shielding,  the  pi utoniuin-to-sulphur  ratio  is  estimated  to  bp  ~9  and  we 
c  10  2 

thereby  obtained  ~2.6x  10  n/cm  (>10keV)  as  the  fluence  for  the  first  burst 
(gamma  dose  s:  60  rads(Si)). 

Figure  4-27  shows  dark  current  density  versus  time  for  Device  J6  follow¬ 
ing  the  two  neutron  pulses  described  above.  Data  are  displayed  for  the  time 
regime  7-2400  sec  after  each  burst.  Following  the  first  pulse,  no  useable 
information  was  obtained  at  times  less  than  7  sec  on  any  of  the  four  oscillo¬ 
scope  channels  employed.  This  occurred  because  of  inappropriate  choice  of  gain 
settings  (It  was  clear,  however,  that  the  dark  current  signal  exhibited 
erratic  behavior  for  t  <  7  sec.)  The  pre-irradiation  value  of  dark  current 
density  for  Device  J6  at  ?98°K  was  4.6  nA/cm2.  At  2400  sec  after  the  first 
neutron  burst  (~2.6  x  10  n/cm  ),  the  value  of  at  that  same  temperature 
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was  5.4  nA/cifl  ,  so  only  slight  permanent  damage  was  incurred.  At  2400  sec 
following  the  second  burst  (4.4  x  10  n/cm  ),  Jd  =  27.7  nA/cm  as  shown  in 
Figure  4-27.  At  7  sec,  Oj  is  a  factor  cf  2.7  larger  than  this  value. 
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Figure  4-27.  Dark  current  density  versus  time  following  pulsed 
neutron  bombardment  for  the  Fairchild  CCD  311  (Device  J6). 

For  times  less  than  7  sec  after  pulse  #2,  several  phenomena  were  observed. 

As  an  aid  to  describing  the  observations,  let  the  value  of  Jd  at  7  sec  be  equal 

to  x.  For  t<  0.2  sec,  no  definitive  information  was  obtained,  but  it  is 
likely  that  the  CCD  wells  were  saturated  at  such  times  due  to  ionising  radiation 

associated  with  the  burst.  At  t  =  0.2  sec,  the  two  shift  registers  in  Device 

J6  had  a  significantly  different  dark  current.  One  register  had  a  relative 
Jd  of  2.5  x  whereas  the  second  had  a  value  of  12  x.  Between  0.2  and  1.9  sec, 

Jd  decreased  monotomically  (and  significantly).  At  t  -  1.8  sec,  for  the  first 
register  Jd  =  13  x  and  for  the  second  register  Jd  =  6  x.  At  t  =  2  sec,  con¬ 
trol  rod  reentry  occurs,  thereby  reducing  the  amount  of  ionizing  radiation 
incident  on  the  device.  Thus,  dark  current  decreases  at  this  point.  (See 
Reference  55  for  additional  discussion  of  the  contribution  of  ionization  to 
Jd  in  CCDs  before  reentry  of  the  control  rod  occurs.)  At  t  =  2.2  sec,  Jd  -  1.2  x 
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I  and  1.6  x  for  the  first  and  second  registers,  respectively.  Between  t  =  2.2 

and  5.8  sec,  Jd  remained  constant  and  the  difference  between  the  two  registers 
was  maintained.  Between  t  =  5.8  and  7  sec,  a  transition  occurred  in  which  the 
|  dark  current  for  both  registers  became  equal  and  also  Jd  began  to  decrease 

|  with  time.  We  interpret  the  data  as  follows.  For  t  a  7  sec,  Jd  is  dominated 

I  by  thermal  generation  of  carriers  at  neutron-produced  centers  in  the  depletion 

j,  region  hulk.  As  t‘:  r  increases,  Jefect  reordering  occurs,  thereby  reducing 

1  the  effective  generation  center  density  (i.e.,  Jd  decreases).  For  t<  2  sec, 

1  is  dominated  by  ionizing  events  for  this  device.  Between  2  and  7  sec, 

apparently  electrical  stability  of  the  CCD  is  being  restored  subsequent  to  the 
large  signals  experienced  at  earlier  times.  (In  this  regard,  we  performed  an 
experiment  on  an  unirradiated  CCD  311  and  observed  a  recovery  time  constant 
of  1  to  2  sec  when  the  deck  frequency  was  abruptly  doubled.  The  initial 
I  frequency  was  low  enough  so  that  the  dark  output  voltage  was  saturated. 

Doubling  the  frequency  reduced  V^,  but  it  took  1  to  2  sec  before  steady-state 
(i.e.,  constant  Vd)  conditions  were  reached.  This  result  indicates  that  one 
must  exercise  caution  when  interpreting  early-time  data  in  short-term  annealing 
>  experiments  on  CCDs.) 
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Since  we  interpret  the  data  of  Figure  4-27  in  terms  of  transient  anneal¬ 
ing  of  neutron  damage,  further  analysis  was  performed  in  a  manner  analogous  to 

previous  studies  of  short-term  recovery  in  neutron-irradiated  silicon  devi- 
58 

ces.  We  define  a  dark-current  annealing  factor  (AF)  as  follows: 


AF  = 


Jd(t)  -  J 


J7T“T  -  J 


do 


do 


(15) 


In  this  equation  Jd(t)  is  dark  current  density  at  any  time  t  following  pulsed 
neutron  bombardment,  JdQ  is  the  pre-irradiation  value,  and  Jd(°°)  is  the  value 
at  some  long  time  following  bombardment  (2400  sec  for  the  data  of  Figure  4-27). 

Figure  4-28  shows  AF  versus  time  for  Device  J6  based  on  the  data  of 
Figure  4-27.  A  substantial  difference  is  noted  between  results  for  pulse  #1 
and  pulse  #2.  One  may  question  the  validity  of  this  difference  since  the 
amount  of  permanent  damage  induced  by  the  first  pulse  was  relatively  small. 
Thus,  in  Equation  (15)  we  have  a  relatively  small  difference  between  relatively 
large  nurrKers  in  the  denominator,  so  a  moderate  degree  of  experimental  error 
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Figure  4-28.  Dark  current  annealing  facccr  versus  time  for  the  ! 

Fairchild  CCD  311.  Data  shown  were  derived  from  Figure  1 

4-27.  j 

i 

i 

can  have  a  large  impact  on  the  annealing  factor.  However,  this  impact  is  j 

stronger  in  the  direction  of  a  larger  AF  as  opposed  to  the  direction  of  closer  j 

agreement  between  the  two  data  sets  in  Figure  4-28.  To  examine  how  much  ; 

experimental  error  would  be  required  for  the  data  obtained  following  pulse  #1 
to  bring  that  data  set  into  better  agreement  with  the  results  obtained  after  \ 

pulse  #2,  we  varied  Jd(®)  and  0^  .  For  example,  assuming  10%  error,  AF  was  j 

recalculated  by  increasing  0d(“)  by  10%  and  decreasing  JdQ  by  10%.  This  con-  ] 

stitutes  u  worst-case  error  situation  in  which  the  AF  calculated  from  measured 
values  overestimates  the  actual  AF.  An  error  bar  is  shown  in  Figure  4-28  at  j 

7  sec  corresponding  to  10%  error.  Also  shown  arc  20%  and  50%  error  bars.  We  j 

estimate  that  the  error  in  0^  measurements  for  our  experiments  was  no  larger  { 

than  ±  20%,  and  most  likely  was  considerably  smaller.  Even  if  the  error  was  \ 

±  50%,  AF  after  pulse  #1  would  still  be  significantly  larger  than  that  obtained  1 

after  pulse  #2.  Hence,  the  discrepancy  evident  in  Figure  4-2G  appears  to  be  j 

real,  but  additional  low-fluence  data  would  be  required  to  confirm  it.  j 

There  is  also  a  potential  problem  with  the  pulse  #2  AF  data  in  Figure  4-28.  j 
Measured  dark  voltage  at  early  times  (i.e.,  7  sec)  was  very  close  to  the  maximum  ] 
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value  that  can  be  obtained  in  normal  operation  from  the  device  studied.  This 
situation  resulted  because  of  the  relatively  low  clock  frequency  employed. 

The  question  arose  whether  the  pulse  #2  AF  data  were  influenced  by  the 
approach  of  dark  voltage  to  near  its  saturation  va?je.  That  is,  we  suspected 
chat,  at  early  times  the  device  was  not  in  a  linear  operating  range  so  that  the 
variation  of  Vd  with  frequency  was  nonlinear.  To  avoid  the  possible  effects  of 
saturation,  additional  measurements  were  made  on  a  second  CCD  311  (Device  J3)  at 
higher  deck  frequencies. 
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Figure  4-29.  Dark  current  density  versus  time  following  pulsed 
neutron  bombardment  for  the  Fairchild  CCD  311  (Device  J3). 


Figure  4-29  shows  J .  versus  time  following  two  identical  neutron  bursts 
C2.5  x  10  n/cm  (>  10  keV),  (2.90  pulse,  ~10  msec  width,  gamma  dcse  a:  330 
rads  (Si Q  for  Device  J3.  Note  that  the  pre-irradiation  dark  current  density 
(73  nA/cm  )  is  considerably  larger  than  that  for  Dev’ce  J6,  so  that  a  consid¬ 
erably  higher  fluence  was  required  for  J3  to  cause  a  significant  increase  in 
dark  current.  For  pulse  #1,  the  clock  frequency  during  bombardment  was  200  kHz. 
This  frequency  was  sufficiently  high  to  keep  Vd  substantially  less  than  its 
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saturation  value  at  early  times  following  bomba rchnent.  At  t  -  15  sec,  the 
clock  frequency  was  decreased  to  20  kHz  for  the  purpose  of  increasing  the 
signal  amplitude  during  later-time  measurements .  For  pulse  #2,  the  clock 
frequency  was  503  kHz  our* ng  bc.HDardment  and  was  decreased  to  50  kHz  at 
t  =  2  sec. 

Figure  4-30  shows  AF  versus  time  for  the  two  data  sets  of  Figure  4-29 
(Device  J3).  Agreement  between  the  results  for  pulse  #1  and  pulse  #2  is  good. 
Also  shown  in  Figure  4-30  are  AF  results  for  Device  J6  (pulse  #2)  taken  front 
Figure  4-28.  For  t  2  20  sec,  agreement  among  the  three  data  sets  is  good, 
which  suggests  that  the  true  dark-current  AF  at,  for  example,  20  sec  is~1.7 
as  compared  to  the  value  of  30  obtained  for  Device  J6  after  pulse  #1  (Figure 
4-28).  There  is  still  the  possibility  of  a  fluence  dependence  of  the  transient 
annealing  behavior  at  low  fluences,  but  such  a  dependence  appears  unlikely  in 
view  of  the  agreement  shown  in  Figure  4-30  for  Devices  J3  and  J6.  (We  also 
note  that  no  fluence  dependence  of  short-term  annealing  in  neutral  bulk  silicon 
was  observed  previously.59) 


Figure  4-30.  Dark  current  annealing  factor  versus  time  for  the 
Fairchild  CCD  311.  Data  shown  were  derived  from  Figures 
4-28  and  4-29. 
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it  is  of  interest  to  compare  the  observed  change  In  at  long  times 

following  pulsed  bombardment  with  the  predicted  change  based  on  data  in  Figure 
4-26,  and  such  a  comparison  is  made  in  Table  4-3.  The  experimental  data 


- — 1 

Device 

# 

Pulse 

# 

—  -  -  -  * 

FI uence 

(n/cm2) 

Measured 

A  J  j 
d  2 
(nA/'cm  ) 

Calculated 

A 

(nA/cm2) 

06 

1 

~2.6  x  1010 

0.80 

0.91 

J6 

2 

4.4  x  1011 

22.3 

15.5 

J3 

1 

2.5  x  1012 

82.4 

87.5 

J3 

2 

_ 

12 

2.5  x  10 

83.3 

87.5 

Table  4-3.  Comparison  of  measured  and  calculated  "permanent"  changes 
in  dark  current  density  for  irradiated  CCD  311s.  Measured 
values  are  based  on  the  value  for  J.  at  times  2000  sec  follow¬ 
ing  pulsed  neutron  bombardment  (see  Figures  4-27  and  4-29). 
Calculated  values  are  based  on  the  slope  shown  in  Figure  4-26 
(3.5  x  10'11  nA/n). 

are  taken  from  Figures  4-27  and  4-29,  and  the  calculated  values  were  obtained 

using  AJ^/A^  *  3.5  x  10_1*  nA/n.  Agreement  between  measured  and  calculated 

values  for  AJd  is  good,  which  demonstrates  that  permanent  increases  in  Jd  by 

neutron  bombardment  can  be  predicted  quite  well  for  the  CCD  311.  Regarding 

prediction  of  the  transient  behavior  of  J^,  the  AF  at  a  given  time  can  be  used 

for  this  purpose.  For  example,  AF  ^  1.7  at  20  sec  following  bombardment,  to 

the  appropriate  value  of  AJ, JL<t>  to  use  at  that  time  is  (1.7)(3.5  x  10  li), 

-11  ° 

or  6  x  10  nA/n.  The  basis  for  this  calculation  is  obvious  upon  inspection 
of  Equation  (15). 

The  dark-current  annealing  factors  shown  in  Figure  4-30  are  significantly 
larger  at  times  >  10  sec  than  those  observed  for  neutral  bulk  regions  in  bulk 
silicon  and  silicon  devices.58’59  For  example,  in  the  neutral -region  case  at 
20  sec  AF  ranges  from  1.0  to~1.2.  For  the  depletion-region  situation  (i.e., 
dark  current  annealing)  AF  ^  1.7  at  that  time.  Short-term  annealing  is 
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strongly  dependent  on  Injection  level  In  neutral  material.®^  On  this  basis, 
one  might  speculate  that  the  higher  AF  in  the  depletion-region  case  Is  due  to 
the  extremely  low  free-carrier  density  in  that  region.  However,  since  the  CCO 
gates  were  clocked  during  the  transient  recovery  period,  most  of  the  depletion 
region  is  collapsed  (l.e.,  becomes  neutral  material)  for  one-half  the  time.  It 
is  possible  that  the  portion  of  the  depleted  region  which  does  not  collapse  Is 
responsible  for  the  relatively  slow  transient  recovery  observed.  Further  study 
will  be  required  to  clarify  this  point. 

4.2.5  Summary. 

In  this  section,  results  of  an  investigation  of  neutron  and  ionizing 

radiation  effects  performed  during  the  second  half  of  this  program  have  been 

presented.  Principal  findings  are  as  follows.  The  HAC  CCD  is  quite  sensitive 

to  the  effects  of  ionizing  radiation  when  irradiated  under  power.  Dark  current 

2  3 

density  for  this  device  exhibits  a  strong  increase  in  the  range  10  -  10  rads, 

4  S 

followed  by  a  significant  decrease  in  the  range  10  -  10  rads.  At  higher 
doses,  Jd  increase  linearly  with  ^ose  and  this  is  attributable  to  production 
of  generation  centers  in  the  depletion-region  bulk  and  to  production  of  elec¬ 
tronic  staces  at  the  Si02-Si  interface.  Dark  current  in  the  Fairchilu  CCD  311 
was  relatively  insensitive  to  ionizing  irradiation,  both  in  unpowered  and 
powered  conditions,  up  to  ~105  rads,  where  device  failure  was  observed.  Stud¬ 
ies  of  the  transient  and  perm  nent  ef'ects  of  neutron  bombardment  on  the  CCD 
311  were  performed.  Thn  :  ate  of  change  of  dark  current  density  w^th  neutron 
flaence  was  found  to  be  3.5  x  10”*1  nA/neutron  for  measurements  made  approxi¬ 
mately  one  hour  aftei  sceaiy-state  neutron  irradiation.  Short-term  annealing 
measurements  were  performed  in  which  dark  current  density  was  monitored  as  a 
function  of  time  foil  iwing  pulsed  bombardment.  As  an  example,  at  20  sec  after 
irradiation  the  dark  current  anncaii^g  factor  was~1.7.  This  me'os  that  the 

appropriate  value  for  AJ jL<t>  *t  that  time  for  the  CCD  311  is  6  x  lO-11 

0  1 

nA/n  as  opposed  to  tie  long-term  value  of  3.5  x  10  nA/n. 
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SECTION  5.0 


NEUTRON  DAMAGE  COEFFICIENT  FOR  SILICON  DEPLETION  REGIONS 


5.1  INTRODUCTION 

A  damage  coefficient  appropriate  for  the  Introduction  by  neutrons  of 
carrier  generation  centers  in  a  silicon  depletion  region  has  been  determined, 
as  described  in  Section  4.1  of  this  report.  There  it  was  suggested  that  this 
damage  coefficient  might  be  generally  useful  for  predicting  the  effects  of 
neutron  bombardment  on  silicon  devices.  In  the  present  section,  we  explore 
that  possibility.  A  presentation  of  the  analytical  background  Is  given  first, 
followed  by  a  comparison  of  calculations  with  experimental  results  obtained  on 
silicon  devices  at  several  laboratories. 

5.2  ANALYTICAL  CONSIDERATIONS 

The  expression  for  the  steady-state  net  recombination  rate  R  at  a 
single  level  in  neutral  bulk  material  is  given  by  the  Nall -Shockley-Read  mode1 
as 


'  Cn(n  +  V  +  cp(P  +  pl}  '  1 

In  this  expression,  p  and  n  are  hole  and  electron  concentrations.  Cp  and  cn  are 
hole  and  electron  capture  probabilities,  n-  is  the  intrinsic  carrier  concentra¬ 
tion.  Nr  is  the  recombination  center  concentration,  and  n^  is  the  electron  con¬ 
centration  that  would  exist  if  the  Fermi  level  were  located  at  the  recombination 
level.  (A  similar  definition  holds  for  p^.)  In  the  space-charge  region  of  a 
reverse-biased  p-n  junction,  p,n  «  n^.  Also,  for  a  level  near  midgap,  nj»Pj  ** 
n^ ,  so  that  r,p  «  n^p^.  Equation  (1)  thus  reduces  to  the  following  expression 

for  the  thermal  generation  rate  G  of  carriers  at  generation  centers  in  a  deple- 
51 

tion  region: 


G 


c  cN_n_. 
n  p  g  i 

:  n,  +  c  p, 
n  1  pKl 


(2) 


(N  is  the  generation  center  concentration.)  Generation  lifetime  t  is  a  charac- 
y  y^i 

teristic  time  associated  with  the  thermal  generation  of  carriers  and  is  given  by" 
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(3) 


V 


-  Cnnl  *  Vl 
2cncpVi 


Thus,  the  g  neration  rate  is  given  by^ 


G  -  -  n./2Td  .  (4) 

(See  additional  discussion  of  this  relation  in  Section  5.4.)  The  current  density 
in  a  reverse-biased  p-n  junction  due  to  thermal  generation  of  carriers  in  a  deple¬ 
tion  region  of  width  is  given  by 


/d  qn,  x. 

q)G|dx  —  qiGjx^  -  *  (5) 

(The  approximation  is  due  to  the  assumption  of  a  uniform  generation  rate  over  the 
entire  deplete i  region.  See  further  discussion  ir  Section  5.4  on  the  validity  of 
Eq.  (5).) 


The  generation  current  density  typically  dominates  reverse  current. 

9  0  51 

flow  in  a  silicon  p-n  junction  below  ~100  C.  In  certain  device  structures,  it 
is  possible  for  carrier  generation  at  interface  states  to  be  an  important  compon¬ 
ent  of  reverse  current.  We  assume  for  simplicity  here  that  dark  (or  leakage) 
current  in  neutron- irradiated  silicon  devices  is  dominated  by  Jg.  The  primary 
effect  of  r.eutron  bombardment  is  to  decrease  Tg  due  to  introduction  of  generation 
centers  in  the  depletion  region.  As  noted  in  Section  4.1,  this  degradation  proc 
ess  can  be  expressed  as 


1 


T 


9 


(6) 


In  this  expression,  and  are  pre-  and  post-irradiation  values  of  generation 
lifetime,  respectively,  0  is  neutron  fluence,  and  is  generation-1 if etime  damage 
coefficient.  Substituting  (6)  into  (5),  we  obtain 


V 


q  n.  x 


Jid  L L  +±\, 

2  K  v 


go 


qni  v 

2  K_ 


(7) 


go  g  /  g 

and  thus  the  neutron-induced  change  in  generation  current  density  is  given  by 


J9  '  J9°  ’  1  "1  “d  *  n  Kg 
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As  described  In  Section  4.1  of  this  report,  capacitance-versus-time 
measurements  Mere  made  on  NOS  capacitors  before  and  after  neutron  bombard¬ 
ment.  Zerbst  analysis  was  used  to  obtain  generation  lifetime,  and  the  re¬ 
sulting  data  are  repeated  here  for  convenience  (Figure  5-1}.  Linear  degra¬ 
dation  of  Tg  Is  evident  over  a  wide  range  In  fission  neutron  fluence.  Also, 
results  for  devices  fabricated  on  n-  and  p-type  material  are  In  agreement. 

The  value  of  K_  which  describes  the  data  In  Figure  5-1  quitr  ‘•♦ell  Is  7.0  x 
6  So 

10  neutron- sec/cm  .  We  also  found  (see  Section  4.1'  that  TgQ  and  t^,  and 
hence  Kg,  are  only  weakly  temperature  dependent.  Thus,  the  value  obtained 
here  for  Kg  appears  to  apply  for  depletion  regions  In  both  n-  and  p-type 
silicon  at  room  temperature  and  a'lso  Is  a  good  value  for  use  over  a  moderate 
temperature  range  (~240  to  ~370°K,  the  range  studied).  The  Implication  is 
that  the  present  Kg  value  is  of  general  usefulness  In  predicting  neutron 
damage  effects  .  depletion  regions  for  any  silicon  device  bombarded  by  fis¬ 
sion  neutrons.  3elow  we  compare  calculations  with  experimental  data  to  ex¬ 
amine  the  validity  of  this  concept. 


It  should  also  be  noted  that  neutron- induced  carrier  removal  can  have 

an  important  effect  in  silicon  p-n  junctions.  As  stated  above,  the  primary 

effect  of  neutrons  on  a  reverse-biased  p-n  junction  Is  to  degrade  t  and  thereby 

9  1 1 

increase  J  .  This  effect  is  important  at  moderate  fluence  levels  (  >10  n/ cm  ), 

9  13  2 

as  shown  in  Figure  5-1.  At  relatively  high  fluences  (e.g.,  >  10  n/cm  ),  car¬ 
rier  removal  can  cause  the  depletion  region  to  widen,  thereby  causing  to 
increase  due  to  the  larger  volume  for  carrier  generation.  Emphasis  here  is 
placed  on  changes  in  due  to  rq  degradation,  but  mention  of  carrier  removal 
effects  is  made  where  appropriate 


The  rate  of  change  of  generation  current  density  with  neutron  fluence 
can  be  expressed,  using  Eq.  (8),  as 


<Tp” 


V 


*  qn 


V2  *9 


■9) 


6  7 

Assuming  a  temperature- independent  value  of  Kg  (7.0  x  10  n-sec/cm  ),  we  can 
calculate  values  of  this  derivative.  Figure  5-2  shows  results  of  such  calcu¬ 
lations  at  270,  300,  and  330°K  as  a  function  of  depletion  region  width.  As  an 
example,  for  a  1-pm  width,  the  value  of  the  right  side  of  Eq.  (9)  (from 
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Figure  5-1.  Neutron- induced  change  in  reciprocal  generation  lifetime  for  MOS 
capacitors  on  bulk  n-  and  p-type  silicon.  These  data  yield  a  genera¬ 
tion  lifetime  damage  coefficient  K  . 
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Figure  5-2)  is  1.6  x  10“**  nanoamps  per  neutron.  Thus,  a 
13  2 

fluence  of  10  n/cm  (>  10  keV)  should  cause  an  increase 

2 

rent  density  of  160  nA/cm  . 


fission  neutron 
in  generation  cur- 
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5.3  COMPARISON  OF  CALCULATIONS  WITH  EXPERIMENTAL  DATA 

Available  experimental  data  are  now  compared  with  results  of  calcula¬ 
tions  using  the  analytical  approach  described  above.  The  data  considered  here 
are:  (1)  CCD  data  presented  in  Section  4.0;  (2)  PIN  photodiode  results  obtained 
by  Kalma  and  Hardwick;6*  (3)  JFET  data  of  Kern  and  McKenzie.6^ 

Figures  4-9  and  4-26  show  the  increase  in  dark  current  density  in 
neutron-irradiated  buried-channel  CCDs.  For  the  Hughes  Aircraft  device  (Fig¬ 
ure  4-9),  the  rate  of  change  in  dark  current  density  with  fluence  is  ^4x 
10”**  nA/neutron.  For  the  Fairchild  CCD311  (Figure  4-26),  the  value  for  this 
rate  is  3.5  x  10" 11  nA/neutron.  Upon  comparing  these  experimental  results  with 
the  calculations  shown  in  Figure  5-2,  depletion  region  widths  o*  —2.2  and 
-2.5  (jfli  are  indicated  for  the  Fairchild  and  Hughes  devices,  respectively.  How¬ 
ever,  calculated  deep  depletion  widths  for  these  device  types  are  considerably 
larger.*  To  make  such  calculations,  we  use  the  expression  for  the  depletion 
width  of  a  step  junction  (n+p)  with  a  reverse  bias  Vr  applied: 


For  the  Hughes  Aircraft  CCD,  the  substrate  dopant  concentration  N  is  —5  x 

14  o  #  d 

10  cm”  ,  corresponding  to  —30  fi-cm  material.  Also,  gate  electrodes  were 
clocked  between  0  and  10V.  We  assume  that  all  of  this  voltage  appears  across 
the  n+p  junction  and  the t  none  of  it  is  dropped  across  the  gate  oxide.  Equa¬ 
tion  (10)  then  yields  5.1  ^m.  The  n-type  epi  layer  (0.8  pai  thick)  is  also  de¬ 
pleted  in  normal  operation  of  the  CCD,  so  the  total  depletion  width  is  -5.9  p,m. 
For  the  Fairchild  CCD311,  the  range  given  fer  the  substrate  resistivity  is  15  to 
60  fi-cm.  Upon  using  Eq.  (10)  with  Vr  =  10V  and  also  adding  the  thickness  of 
the  n-type  epi  layer  (0.5  n,m),  we  obtain  a  range  of  -4.1  to  —6.5  p-m  for  the 
total  depletion  width,  with  the  mean  value  being  -5.3  p,m.  Thus,  the  ratio  of 
calculated  x^  (in.  (10))  to  apparent  x^  (Figure  5-2)  is  —2.4  for  CCDs  from 
both  suppliers.  Stated  another  way,  there  is  an  apparent  discrepancy  by  a 
factor  of  -2.4  between  CIO  dark  current  measurements  and  calculations  of  the 
neutron-induced  change  in  this  quantity  using  Eq.  (9).  Possible  reasons  for 
this  discrepancy  are  discussed  below  in  Section  5.4. 

★ 

The  present  calculations  of  x^  are  considered  more  accurate  than  the  estimated 
range  given  in  our  earlier  work  (reported  in  Section  4.1). 
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Kalma  and  Hardwick®*  measured  the  radiation  response  of  various  com¬ 
mercially  available  PIN  photodiodes.  Included  in  their  study  were  measurements 
of  neutron- induced  increases  in  dark  current  for  such  devices.  Figure  8  in 
their  pacer  Is  a  plot  of  dark  current  versus  neutron  fluence  for  diodes  from 
seven  manufacturers.  Table  5-1  In  the  present  work  lists  the  device  types 
studied  by  Kalma  and  gives  additional  Information.  The  reverse  bias  at  which 
dark  current  measurements  were  made  and  the  measured  pre-irradiation  capaci- 

CO 

tance  at  that  bias  are  listed  for  each  device.  These  values  of  bias  and 

capacitance  are  those  appropriate  for  the  data  shown  in  Figure  8  of  Kalma 's 
61  63 

paper.  ’  It  is  assumed  here  that  the  measured  capacitance  is  attributable 
solely  to  the  depletion  layer  of  width  in  a  given  PIN  diode.  This  width 
was  then  calculated  using  the  relation  -  eA/C,  where  €  is  dielectric  con¬ 
stant,  A  is  active  area,  and  C  is  capacitance.  Results  of  such  calculations 
are  listed  in  Table  5-1  for  each  device  type. 

The  next  step  was  to  replot  Raima's  data,  for  convenience  of  analysis, 
in  terms  of  the  change  in  dark  current  density  (J^  -  JdQ)  versus  neutron  flu¬ 
ence  and  Figures  5-3  to  5-6  present  the  results  of  this  replotting.  Also 
shown  in  these  figures  (solid  and  dashed  lines)  are  results  of  calculations 
performed  using  Eq.  (8).  Table  5-2  gives  the  ratio  of  calculated  to  experi¬ 
mental  values  of  the  quantity  Jd  -  0dQ  as  a  function  of  neutron  fluence  based 
on  information  in  Figures  5-3  to  5-6.  We  now  discuss  findings  for  each  device 
type. 


Figure  5-3  compares  calculation  with  experiment  for  Spectronics  devices. 

For  the  SPX-3290,  agreement  is  quite  good  at  the  two  highest  fluences.  At 
12  2 

10  n/cm  ,  the  calculated  change  in  dark  current  density  is  63%  of  the  measured 
11  2 

value.  At  10  n/cm  ,  calculation  yields  a  value  that  Is  42%  of  the  measured 
change  in  current.  For  the  SPX-2232,  agreement  is  good  at  the  highest  fluence. 
At  lower  fluences,  calculated  values  are  32  to  38%  of  those  measured.  The  in¬ 
dication  is  that  some  additional  current  component  (i.e.,  other  than  that  due 
to  carriers  thermally  generated  in  the  depletion  region)  contributes  at  lower 
fi  uences. 

Figure  5-4  presents  experimental  data  and  calculations  for  the  Hewlett- 

Packard  PIN  diode.  Agreement  is  good  at  the  highest  fluence.  (It  should  be 

6 1 

noted  that  in  Raima's  paper  there  is  a  plotting  error  at  the  highest  fluence 
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DATA  POINTS  AFTER  KALMA 
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-  44pm 

a  SPX-3290 

- CALCULATION 

-8.3  pm 


A  / 

/ 

/ 

/ 

/ 


A  / 

/ 

/ 

/ 

/ 


0!1  1012  1013  io'4 

FLUEMCE  (n/cm2) 

Figure  5-3.  Change  in  dark  current  density  versus  neutron  fluence  for  Spectronics 
PIN  diodes.  Data  points  were  derived  from  the  work  of  Kalma  and  Hardwick. 61 
Results  of  calculations  are  also  shown. 
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Figure  5-4.  Change  in  dark  current  density  versus  neutron  fluence  for  a  Hewlett 
Packard  PIN  diode.  Data  points  were  derived  from  the  work  of  Kalma  and 
Hardwick. Results  of  calculations  are  also  shown. 
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Figure  5-5.  Change  in  dark  current  density  versus  neutron  fluence  for  United 
Detector  Technology  PIN  diodes.  Data  points  were  derived  from  the  work 
of  Kalma  and  Hardwick. 61  Results  of  calculations  are  also  shown. 
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Table  5-2.  The  ratio  of  calculated  to  experimental  values  cf  the  quantity 
Jd  ’  Jdo  at  various  neutron  fluences  based  on  Information  in  Figures 
5-3  to  5-6. 


DEVICE 

1011  n/cm2 

1012  n/cm2 

1013  n/cm2 

1014  n/c*2 

SPX-2232 

0.32 

0.36 

0.38 

1.18 

SPX-3290 

0.42 

0.63 

0.86 

1.2 

5082-4207 

0.23 

0.30 

0.54 

0.86 

PIN-050 

0.39 

0.97 

2.15 

14.5 

PIN-020A 

— 

1.00 

1.98 

0.78 

YAG-100 

0.25 

0.43 

1.16 

2.54 

SHS-100 

0.35 

0.45 

1.28 

2.32 

AVERAGE 
(excludes  P 

.  — ■ 

0.31 

IN-05D) 

_ 

0.53 

1.03 

1.48 

in  Figure  8  for  the  Hewlett-Packard  device.  The  data  point  shown  here  in 

CO 

Figure  5-4  is  the  correct  one.  )  Disagreement  is  noted  at  lower  fluences, 
with  the  discrepancy  decreasing  as  the  fluences  increase.  Qualitatively,  one 
can  account  for  this  discrepancy  in  terms  of  an  additional  current  component 
that  saturates  at  some  moderate  fluence  value.  An  example  of  such  a  saturating 
current  would  be  that  attributable  to  gamma-ray- induced  charge  buildup  in  a 
surface  oxide  layer.  This  charge  buildup  would  lead  to  enhanced  thermal  genera¬ 
tion  of  carriers  due  to  interface  state  production  and/or  depletion  of  the  under¬ 
lying  silicon.  Charge  buildup  in  SK^  does  saturate,  but  further  study  would  be 
required  to  determine  the  mechanism  responsible  for  the  enhanced  current  at  low 
fluences  evident  in  Figure  5-4. 

Comparisons  between  calculations  and  experimental  data  for  United  Detector 
Technology  devices  are  shown  in  Figure  5-5.  For  the  PIN-05D,  the  data  tend  toward 
saturation  at  high  fluences.  The  calculated  line  (dashed)  describes  the  data 
fairly  well  at  the  three  lowest  fluences  but  not  at  10^  n/cm2.  At  that  fluence, 
Kalma  noted  a  significant  (i.e.,  factor  of  ~4)  drop  in  capacitance  for  this  de- 
vice,  apparently  due  to  carrier  removal.  Since  the  PIN-05D  was  not  operated  in 
full  depletion,  carrier  removal  should  result  in  a  wider  depletion  region  and  thus 
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14  2 

an  Increase  in  dark  current.  However,  the  measured  current  at  10  n/cm  is  sub¬ 
stantially  less  than  would  be  expected  if  no  change  in  depletion  width  had  occur¬ 
red.  (Note  that  a  constant,  pre- irradiation,  value  of  xd  was  used  to  perform 

calculations.  This  is  true  in  all  cases  shown  in  Figures  5-3  to  5-6.)  For  the 

12 

PIN-020A,  agreement  between  calculation  and  experiment  is  good  at  10  and 
14  2 

10  n/cm  ,  with  a  factor  of  ~2  discrepancy  occurring  at  the  intermediate  fluence. 
Kalma  noted  a  slight  (25%)  decrease  in  measured  capacitance  for  this  device  after 
1014  n/cmS  presumably  due  to  carrier  removal.®3  If  the  lower  value  of  capacitance 
is  used  to  calculate  a  new  value  of  xH.  then  the  agreement  between  experiment  and 
calculation  is  within  5%  at  101H  n/cm  for  the  PIN-020A.  In  view  of  the  scatter 
in  the  data  for  United  Detector  devices,  such  close  agreement  is  probably  fortuitous. 

Comparison  of  calculations  with  experimental  data  for  EG&G  devices  are 

shown  in  Figure  5-6.  Calculated  results  underestimate  the  data  at  the  two  lowest 

13  2 

fluences  and  overestimate  the  data  at  the  two  highest  fluences.  At  10  n/cm  , 
agreement  between  theory  and  experiment  is  within  ~25%,  with  the  discrepancy 
being  considerably  larger  at  other  fluences. 

Kern  and  McKenzie®^  measured  leakage  current  in  JFETs  before  and  after 

14.8-MeV  neutron  bombardment.  This  current  was  measured  between  source  and 

drain  while  a  gate  bias  of  sufficient  magnitude  to  deplete  the  active  volume 

was  applied.  The  primary  source  of  leakage  current  in  their  devices  should  be 

thermal  generation  of  carriers  in  the  depletion  region.  At  300°K,  the  increase 

12  2  -12 

in  leakage  current  after  1.1  x  10  n/cm  was  3  x  10  A.  The  volume  of  the 

o  3 

devices  studied,  assumed  to  be  fully  depleted,  was  2.04  x  10  cm  .  Equation  (8) 
can  be  expressed  ir.  terms  of  generation  current  instead  of  current  density: 

!g  '  !go  =  q  ni  0  *d  A/2  Kg  •  (11) 

(A  is  area.)  We  equate  the  x^A  product  to  the  volume  of  the  depletion  region 

given  by  Kern  and  McKenzie.  In  order  to  calculate  I  -  I  a  value  of  K„  ap- 

g  go  g 

propriate  for  14.8  MeV  neutrons  is  needed.  For  recombination  center  production 

64 

in  neutral  bulk  silicon,  fusion  neutrons  are  a  factor  of  2  to  3  move  effective. 
Assuming  that  14.8-MeV  neutrons  are  a  factor  of  2  more  effective  than  fission 
neutrons  in  producing  generation  centers,  K_  is  then  reduced  to  3.5  x 

GO  9  _  1  o 

10  n-sec/cm  .  Equation  (11)  then  yields  a  value  of  7.2  x  10"1  A  for  change 
in  generation  current,  which  is  a  factor  of  2.4  larger  than  that  measured  by 
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Kern  and  McKenzie.  It  is  interesting  to  note  that  this  discrepancy  factor  is 
equal  to  that  noted  above  for  CCDs  when  experimental  results  were  compared 
with  calculations.  Further  discussion  of  this  point  is  given  below.  One  ad¬ 
ditional  comment:  Kern  and  McKenzie's  data  indicate  that  neutron  bombardment 
produces  a  generation  level  very  close  to  the  level  which  dominates  before  ir¬ 
radiation.  The  same  conclusion  was  reached  in  our  work  (reported  in  Section  4.1). 

5.4  DISCUSSION  AND  CONCLUSIONS 

We  first  summarize  the  above  comparisons  of  calculations  with  experi¬ 
mental  results.  For  CCDs  fabricated  by  two  sources,  the  calculated  neutron- 

induced  increase  in  dark  current  density  was  larger  than  the  experimental  value 

62 

by  a  factor  of  2.4.  For  the  JFET  data  of  Kern  and  McKenzie,  our  calculated 

increase  in  leakage  current  is  larger  than  their  experimental  value  by  a  factor 

of  2.4,  assuming  a  fusion-fission  damage  ratio  of  2.  For  the  PIN  photodiode 

fil 

data  of  Kalma  and  Hardwick,  results  are  varied.  As  shown  in  Table  5-2,  at 

11  12  2 

low  fluences  (i.e.,  10  and  10  n/cm  )  the  calculated  increases  in  dark  cur¬ 
rent  underestimate  the  experimental  values.  At  the  two  highest  fluences,  the 
calculated  value  is  larger  than  the  experimental  value  in  most  cases.  We  have 
listed  average  ratios  at  each  fluence  in  Table  5-2  to  provide  a  rough  indica¬ 
tion  of  the  degree  of  agreement  between  experiment  and  calculation  for  the  PIN 

diodes.  (Because  of  the  anomalously  low  measured  dark  current  for  the  PIN-05D 
14  2 

after  10  n/cm  ,  ratios  for  this  device  type  were  omitted  in  calculating  average 

values.)  The  average  discrepancy  is  a  factor  of  3  at  the  lowest  fluence  and  a 

13 

factor  of  1.5  at  the  highest  fluence.  Overall,  agreement  is  better  at  10  and 

10  n/cm  than  at  the  two  lowest  fluences.  Note  that  in  the  most  successful 

cases  (SPX-3290  and  PIN-020A),  agreement  between  theory  and  experiment  is  within 

a  factor  of  «~2  when  all  f  luences  examined  are  considered.  Also  note  that  in  Eq.  (8), 

which  was  used  to  perform  dark  current  calculations,  the  scaling  factor  on  the 

right-hand  side  is  x,,  Upon  normalizing  the  change  in  dark  current  density  to 

o  61 

this  quantity,  the  data  of  Kalma  and  Hardwick  should  become  independent  of  de¬ 
vice  type  if  it  is  assumed  that  the  present  model  fully  describes  their  results. 
Figure  5-7  shows  the  results  of  such  a  scaling.  Also  shown  is  the  calculated 
value  of  (0d  -  ^o^d  versus  fluence.  The  discrepancies  evident  between  cal¬ 
culation  and  experiment  in  figure  5-7  were  discussed  above  in  connection  with  re¬ 
sults  shown  in  Figures  5-3  through  5-6  and  in  Table  5-2. 
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rigure  5-7.  The  quantity  -  J^o^d  versus  neutron  fluence  for  various  PIN 

diodes.  Data  were  derived  from  the  work  of  Kalma  and  Hardwick6*  as  de¬ 
scribed  in  the  text.  Values  of  x^  employed  are  given  in  Table  5-1. 
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(The  approach  used  tc  obtain  the  calculated  fit  in  Figure  5-7  should 
be  of  general  applicability  for  devices  In  which  measured  dark  current  is 
dominated  by  generation  current  (l.e.,  originates  In  the  depleted  region). 
We  can  express  the  change  in  dark  current,  Ald,  per  unit  volume  V  and  per 
unit  neutron  fluence  as  (using  Eq.  (11)) 


V0  2  Kg 


(12) 


At  300°K,  we  obtain  a  value  for  Al(j/V0  of  1.6  x  10"^  amps/neutron-cm  for  fis¬ 
sion  neutrons.  Multiplying  this  value  by  the  depletion-region  volume  and  by 
the  neutron  fluence  will  give  the  predicted  change  in  dark  current  for  a  given 
device. ) 


In  the  process  of  attempting  to  determine  possible  explanations,  for  the 
discrepancies  noted  above  for  various  types  of  silicon  devices,  we  noted  differ¬ 
ences  in  analytical  relationships  ror  generation  rate  and  generation  lifetime  in 

51 

the  literature.  Grove  gives  the  expression  for  G  stated  above  (Eq.  (2)),  and 

then  equates  this  expression  to  n1f'2Tg  (Eq.  (4)),  which  leads  to  the  Tg  relation 

given  in  Eq.  (3).  That  is,  Grove  states  that  G  =  n./2Tg  in  general  for  a  single 

generation  level.  (Note  that  he  does  not  assume  this  level  to  be  at  midg3p.)  On 

65 

the  other  hand,  Sze,  who  also  employs  the  same  G  expression  as  Grove  (Eq.  (2) 
above),  states  that  G  =  n./Tg  general.  He  then  obtains  an  expression  for  Tg 
which  differs  by  a  factor  of  2  from  that  civen  by  Grove  (Eq.  (3)  above).  It  is 
not  our  intent  to  resolve  this  conflict  here,  but  it  does  have  an  impact  on  the 
Tg  and  Kg  values  obtained  in  the  Present  work.  In  Eq.  (1)  of  Section  4.1,  a 
differential  equation  is  given  which  is  based  on  the  modification  of  Zerbst 
analysis  performed  by  Schroder  and  Nathanson.  That  equation  actually 
differs  slightly  from  the  one  devE loped  by  Schroder  and  Nathanson  because 
they  assume  G  =  n./-rg  whereas  Eq.  (1)  in  Section  4.1  was  developed41  as¬ 
suming  G  =  ri^/2rg.  Thus,  we  have  been  consistent  in  our  analytical  approach 
here  in  that  the  same  G  expression  was  used  in  the  Zerbst  analysis  used  to 
obtain  t  (and  hence  K  )  as  in  the  analysis  used  to  obtain  an  expression  for 
Jg  (Eq-  h)  above).  If  instead  we  had  assumed  G  =  n.j/Tg,  then  all  Tg  values 
and  also  the  value  for  K  would  be  a  factor  of  2  larger.  That  is,  K  would 

7  2 "  1 

equal  1.4 x 10  n-sec/cm  .  Thus,  care  must  be  employed  when  using  K  -  7  x 

fi  ?  ^ 

10  n-sec/cm  ;  one  must  assume  Gsn./2Tg  when  doing  so.  It  is  important  to 


132 


note  that  calculations  performed  using  Fqs.  (8)  and  (9)  (e.g.,  the  calculations 
shown  In  Figure  5-2)  are  independent  of  th*»  6  expression  assumed  as  Iona  as 
the  same  expression  is  used  consistently  (i.e.,  in  Zerbst  analysis  and  in  sub¬ 
sequent  calculations).  We  therefore  conclude  that  the  generation  rate  expres¬ 
sion  assumed  here  does  not  account  for  the  discrepancies  noted  above  between 
calculations  and  experimental  data  for  various  devices. 

66  67 

Grimmeiss  ’  has  performed  experimental  work  and  analysis  which  in¬ 
dicates  that  for  levels  deep  in  the  forbidden  gap  the  effective  width  for 
carrier  generation  in  a  reverse-biased  p-n  junction  is  less  than  the  full  de¬ 
pletion  width.  (This  topic  has  also  been  addressed  by  Calzolari  and  Graff i 
If  such  a  situation  applies  for  the  experimental  cases  considered  above,  then 
a  portion  of  the  discrepancy  between  calculation  and  experiment  can  be  accounted 
for.  In  the  case  of  CCDs,  the  maximum  reverse  voltage  appearing  across  the 
junction  between  the  n-type  epi  layer  and  the  p-type  substrate  is  10V.  Braun 
and  Grimmeiss^  calculate  for  a  linearly  graded  phn  junction  (not  the  actual 
case  here,  but  close  enough  for  the  present  purposes)  with  10V  applied  that  the 
percentage  of  the  depletion  region  width  effective  for  carrier  generation  is 
80%.  This  factor  reduces  the  discrepancy  between  experiment  and  calculation 
for  CCDs  to  a  factor  of  ~1.9.  (If  we  further  assume  that  some  portion,  say 
15**,  of  the  applied  gate  voltage  is  dropped  across  the  oxide,  the  discrepancy 
factor  is  reduced  to  ~1.7.) 

For  the  JFET  experiments  of  Kern  and  McKenzie,62  a  DC  reverse  bias  of 
5V  was  applied  to  the  gate  p-n  junction  to  deplete  the  device.  Once  again 
using  the  linearly  graded  p  n  junction  calculations  of  Braun  and  Grimmeiss,6^ 
we  obtain  an  effective  carrier  generation  width  of  75"  of  the  full  depletion 
region  width.  Applying  this  correction  reduces  the  discrepancy  factor  between 
theory  and  experiment  to  ~1.0  . 

For  the  PIN  diode  data  of  Kalma  and  Hardwick,6^  it  is  not  presently 
clear  how  to  make  an  appropriate  correction  for  the  effective  generation 
volume  (However,  see  Footnote  10  in  Ref.  69.)  As  seen  in  Table  5-2,  at  the 
two  lowest  fluences  such  a  correction  would  serve  to  increase  the  discrepancy 
between  calculation  and  experiment.  It  appears  that  an  additional  current 
component,  i.e.,  other  than  generation  current,  is  important  at  such  fluences. 

At  the  two  highest  fluences,  a  generation  volume  correction  would  decrease  the 
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observed  discrepancy  in  most  cases.  One  further  note:  Leakage  current  meas¬ 
urements  on  PIN  diodes  were  performed  at  room  temperature*  which  may  have 

6  3 

varied  a  few  degrees.  This  could  account  for  some  of  the  discrepancy  be¬ 
tween  calculations  ana  experimental  data  since  depends  strongly  on  tempera¬ 
ture  through  n^.  We  also  note  that  the  active  areas  listed  In  Table  5-1,  which 
were  used  in  all  calculations,  were  not  measured  but  were  obtained  from  manu¬ 
facturers'  data  sheets.  Actual  areas  could  differ  somewhat  from  these  typical 
values  for  specific  devices  and  therefore  increase  the  discrepancy  between 
calculated  and  experimental  results  for  PIN  diodes. 

To  conclude,  we  have  found  that  leakage  current  data  obtained  for  various 
nrut.ron-irradleted  silicon  devices,  including  CCDs,  PIN  diodes,  and  JFETs,  can 
be  described  reasonably  well  using  a  model  which  incorporates  K^,  the  generation 
lifetime  damage  coefficient,  determined  using  MOS  capacitors.  Agreement  between 
calculations  and  experimental  results  Is  within  a  factor  of  <£,2  for  most  cases, 
with  the  exception  of  PIN  diode  data  at  low  fluences  (10^  n/cm^).  Such  agree¬ 
ment  suggests  that  the  present,  calculational  approach  Is  suitable  for  making 
reasonably  accurate  estimates  of  neutron-induced  increases  in  dark  current  for 
silicon  devices.  From  a  basic  mechanisms  viewpoint,  further  work  will  be  re¬ 
quired  to  resolve  the  remaining  di screp&:«cy  between  calculated  and  experimental 
results. 
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George  C.  Messenger 
Consulting  Engineer 

ATTN:  G.  Messenger 

Georgia  Institute  of  Technology 
Georgia  Tech.  Research  Institute 
ATTN:  R.  Curry 

Georgia  Institute  of  Technology 
Office  of  Contract  Administration 

ATTN:  Res.  &  Sec.  Coord,  for  H.  Denny 

Goodyear  Aerospace  Corporation 
Arizona  Division 

ATTN:  Security  Control  Station 

Grumnan  Aerospace  Corporation 
ATTN:  J.  Rogers 

GTE  Sylvania,  Inc. 

Electronics  Systems  Group-Eastern  Division 
ATTN:  C.  Thornhill 
ATTN:  L.  Blaisdell 
ATTN:  L.  Pauples 


GTE  Sylvania,  Inc. 

ATTN:  H  &  V  Group 
ATTN:  H.  Ullman 
ATTN:  P.  Fredrickson 
ATTN:  J.  Waldron 

Harris  Corporation 

Harris  Semiconductor  Division 
ATTN:  J.  Cornell 
ATTN:  C.  Anderson 

Honeywell ,  Inc. 

Avionics  Division 
ATTN:  R.  Gunn 

Honeywell,  Inc. 

Avionics  Division 

ATTN:  C.  Cerulli 

Honeywell,  Inc. 

Radiation  Center 

ATTN:  Technical  Library 


Honeywell,  Inc. 
Defense  Systems 

Division 

ATTN: 

K. 

Gaspard 

Hughes  Aircraft 

Company 

ATTN: 

R. 

McGowan 

ATTN: 

J. 

Singletary 

Hughes  Aircraft 

Company 

ATTN: 

E. 

Smith 

ATTN:' 

W. 

Scott 

IBM  Corporation 

ATTN: 

T. 

Martin 

ATTN: 

F. 

Fietse 

ATTN: 

H. 

Mathers 

I IT  Research 

Institute 

ATTN: 

I. 

Mindel 

Institute  for  Defense  Analyses 

ATTN:  Technical  Information  Services 

Intel  Corporation 

ATTN:  M.  Jordan 

International  Business  Machine  Corporation 
Thomas  Watson  Research  Center 
ATTN:  J.  Ziegler 

International  Telephone  &  Telegraph  Corporation 
ATTN:  Department  608 
ATTN:  A.  Richardson 

Intersil,  Inc. 

ATTN:  D.  MacDonald 

IRT  Corporation 

ATTN:  J.  Harrity 

JAYCOR 

ATTN:  T.  Flanagan 
ATTH:  R.  Stahl 
ATTN:  L.  Scott 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  P.  Partridge 


DEPARTMENT  OF  0EFE.1SE  CONTRACTORS  (Continued) 

National  Academy  of  Sciences 
National  Materials  Advisory  Board 
ATTN:  R.  Shane 


daman  Sciences  Corporation 
ATTN:  J.  Lutoell 

Litton  Systems,  Inc. 

Guidance  A  Control  Systems  Division 
ATTN:  G.  Maddox 

Lockheed  Missiles  A  Space  Company,  Inc. 
ATTN:  C.  Thompson 
ATTN:  H.  Phillips 
ATTN:  E.  Smith 
ATTN:  P.  Bene 
ATTN:  M.  Smith 

Lockheed  Missiles  &  Space  Company,  Inc. 
ATTN:  J.  Smith 
ATTN:  J.  Crowley 

M.I.T.  Lincoln  Laboratory 
ATTN:  P.  McKenzie 

Martin  Marietta  Corporation 
ATTN:  W.  Brockett 
ATTN:  H.  Cates 
ATTN:  W.  Janocko 
ATTN:  R.  Gaynor 

Martin  Marietta  Corporation 
ATTN:  E.  Carter 


National  Semiconductor  Corporation 
ATTN:  R.  Hang 
ATTN:  A.  London 

University  of  New  Mexico 
Electrical  Engineering  A  Computer  Science  Dept. 
ATTN:  H.  Southward 


Northrop  Corporation 
Northrop^ Research  A  Technology  Center 


ATTN: 
ATTN: 
ATTN: 
ATTN: 
ATTN: 
2  cy  ATTN: 
2  cy  ATTN: 
15  cy  ATTN: 


Jackson 
Elsenberg 
talma 
Hartmann 
Chen 
Othmer 
Library 
J.  Srour 


Northrop  Corporation 
Electronic  Division 

ATTN:  0.  Strobel 
ATTN:  L.  Apodaca 
ATTN:  P.  Gardner 


Physics  International  Company 
ATTN:  J.  Shea 
ATTN:  Division  6000 
ATTN:  J.  Huntington 


Magnavox  Govt.  A  Indus.  Electronics  Company 
ATTN:  W.  Rlcheson 

McDonnell  Douglas  Corporation 
ATTN:  Technical  Library 

McDonnell  Douglas  Corporation 
ATTN:  M.  Stitch 
ATTN:  D.  Dohm 
ATTN:  Library 


RAD  Associates 

ATTN:  C.  MacDonald 
ATTN:  R.  Poll 
ATTN:  S.  Rogers 

Rand  Corporation 

ATTN:  C.  Crain 

Raytheon  Company 

ATTN:  J.  Clccio 


McDonnell  Douglas  Corporation 
ATTN:  J.  Holmgrem 
ATTN:  D.  Fitzgerald 


Raytheon  Company 

ATTN:  A.  Van  Doren 
ATTN:  H.  Flescher 


Mission  Research  Corporation 
ATTN:  C.  Longmire 

Mission  Research  Corporation- San  Diego 
ATTN:  J.  Azarewlcz 
ATTN:  V.  Van  Lint 
ATTN:  J.  Raymond 
ATTN:  R.  Berger 

Mitre  Corporation 

ATTN:  M.  Fitzgerald 

Motorola,  Inc. 

Government  Electronics  Division 
ATTN:  A.  Christensen 


RCA  Corporation 
Government  Systems  Division 
Astro  Flectronlcs 

ATTN:  V.  Mancl no 
ATTN:  G.  Brucker 

RCA  Corporation 
David  Samoff  Research  Center 
ATTN:  D.  O'Connor 
ATTN:  Office  N103 

RCA  Corporation 
Government  Systems  Division 
Missile  &  Surface  Radar 
ATTN:  R.  Kill  Ion 


Motorola,  Inc. 
Semiconductor  Group 
ATTN:  L.  Clark 


RCA  Corporation 
Camden  Complex 

ATTN:  J.  Saul tz 
ATTN:  E.  Van  Keuren 
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DEPARTMENT  Of  DEFENSE  CONTRACTORS  (Continued) 
RCA  Corporation 

Somerville  Plant,  Solid  State  Division 
ATTN:  W.  Allen 

Rensselaer  Polytechnic  institute 
ATTN:  R.  Gutmann 

Research  Triangle  Institute 

ATTN:  Security  Office  for  M.  Simons,  Jr. 

Rockwell  International  Corporation 


ATTN: 

V. 

Strahan 

ATTN: 

J. 

Bell 

ATTN: 

T. 

Oki 

ATTN: 

G. 

Messenger 

ATTN: 

V. 

De  Martino 

Rockwell  International  Corporation 
Space  Division 

ATTN:  0.  Stevens 

Rockwell  International  Corporation 
ATTN:  T.  Yates 
ATTN:  TIC  BA08 

Sanders  associates,  Inc. 

ATTN:  L.  Brodeur 

Science  Applications,  Inc. 

ATTN:  J.  Naber 
ATTN:  V.  Ophan 
ATTN:  V.  Verblnski 

Science  Applications,  Inc. 

ATTN:  W.  Chadsey 

Science  Applications,  Inc. 

ATTN:  D.  Stribling 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Sperry  Rand  Corporation 

Sperry  Flight  Systems 
ATTN:  D.  Schow 

Sperry  Uni  vac 

ATTN:  J.  Inda 

Spire  Corporation 

ATTN:  R.  Little 

SRI  International 

ATTN:  P.  Dolan 
ATTl:  A.  Whitson 

Systron  Donne*-  Corporation 
ATTN:  .  Indelicate 

Teledyne  Ryan  Aeronautical 
ATTN:  J.  Rawlings 

Texas  Instruments,  Inc. 

AiTN:  R.  Stehlin 
ATTN:  A.  Peletler 

TRW  Defense  &  Space  Systems  Group 
ATTN:  A.  Pavel ko 
ATTN:  A.  Wltteles 
ATTN :  0 .  Adams 
ATTN:  R.  Kings land 
ATTN:  R.  Schnieder 
ATTN:  H.  Holloway 
ATTN:  H.  Hald 
ATTN;  P.  Guilfoyle 

TRW  Defense  &  Space  Systems  Group 
ATTN:  M.  Gorman 
ATTN:  F.  Fay 
ATTN:  R.  Kltter 


Singer  Company 

ATTN:  J.  Brinkman 

Singer  Company 

Data  Systems 

ATTN:  R.  Spiegel 

Solar  Energy  Research  Institute 
ATTN:  A.  Stanley 

Sperry  Rand  Corporation 

Sperry  Microwave  Electronics 

ATTN:  Engineering  Laboratory 

Sperry  Rand  Corooration 

Sperry  Division 

ATTN:  P.  Maraffino 
ATTN:  R.  Viola 
ATTN:  F.  Scaravagl ione 
ATTN:  C.  Craig 


TRW  Systems  and  Energy 
ATTN:  D.  Millward 
ATTN:  G.  Spehar 

Vought  Corporation 
ATTN:  Library 
AITN:  Technical  Data  Center 
ATTN:  R.  Tomme 

Westinghouse  Electric  Corporation 

Aerospace  S  Electronic  Systems  Division 
ATTN:  L.  McPherson 

Westinghouse  Electric  Corporation 

Defense  and  Electronic  Systems  Center 
ATTN:  D.  Crichi 
ATTN:  H.  Kalapaca 


